<title>Background</title> <p>Development is a variable and time-dependent process that can be thought of as mapping genotypes on to adult phenotypes [##REF##11845200##1##,##REF##17227362##2##]. The time elapsed from the embryo to the reproductive phase, commonly known as developmental time (DT), is a trait of great relevance to fitness in all organisms [##UREF##0##3##]. In particular, <italic>Drosophila </italic>species occupy ephemeral habitats, such as rotting fruits that may result in selection for rapid development. Quoting Gould's <italic>Ontogeny and Phylogeny: </italic>"The timing of maturation is a primary variable in setting life history strategies. We have a prima facie case for ascribing direct significance to the change in developmental timing itself, not only to its morphological consequences"[##UREF##1##4##]. In addition, changes in the timing of developmental processes – heterochrony- may partly account for many evolutionary changes we observe because of genetic correlations of DT with other life-history related traits. Incidentally, heterochrony can also proceed by truncation in addition to retardation/acceleration of developmental processes [##REF##14756337##5##]. The genes that control development are broadly shared among distantly related groups, and some of the pathways are surprisingly similar in superficially different organisms [##UREF##2##6##,##REF##12509752##7##]. In this sense, experimental studies of egg-to-adult heterochrony at the molecular level can provide insights into the evolutionary processes by uncovering the genetic basis of ecologically relevant adaptive traits.</p> <p>DT is a complex trait that displays considerable genetic variation in natural populations [##REF##12583594##8##,##REF##16674586##9##]. Laboratory selection for reduced DT triggered a negatively correlated response of other life-history traits, such as adult weight at eclosion [##UREF##3##10##], adult size [##UREF##3##10##, ####UREF##4##11##, ##UREF##5##12##, ##REF##11525460##13####11525460##13##], pre-adult survival [##UREF##5##12##,##REF##11525460##13##] and longevity [##UREF##6##14##]. The reduction of adult fitness as a consequence of the acceleration of DT is a pattern called <italic>fast development syndrome </italic>and it illuminates the direct connection of the pre-adult and adult stages through energetic trade-offs [##UREF##6##14##]. These studies highlight the negative genetic correlations between DT and a large number of adult life-history traits, suggesting that genes involved in early development may also have pleiotropic effects on adult traits. However, specific genes contributing to variation in DT, and its correlations with other life-history traits, remain largely unknown.</p> <p>The analysis of the genetic architecture of a trait requires not only the identification of the set of genes involved in its expression, but also their response to environmental variation (phenotypic plasticity) [##REF##15172667##15##, ####REF##16755288##16##, ##REF##17028343##17####17028343##17##]. Phenotypic plasticity refers to the ability of a genotype to produce alternative morphological, behavioral and physiological characteristics in response to environmental conditions [##UREF##7##18##, ####UREF##8##19##, ##REF##16731814##20##, ##REF##17418441##21####17418441##21##]. Temperature is an important environmental factor that affects the developing organism and surely has played a major role during the evolution of developmental traits in ectotherms [##UREF##9##22##, ####UREF##10##23##, ##REF##16942614##24####16942614##24##]. Indeed, ectotherms must adapt their developmental program to a wide range of environmental temperatures during the day and also in different seasons along the year.</p> <p><italic>Drosophila </italic>brings an impressive toolkit for dissecting multiple interacting loci with individually small and environmentally sensitive effects that affect complex traits [##REF##15172667##15##]. <italic>P</italic>-element mutagenesis is an effective strategy for studying the effects of single mutations on complex phenotypes [##REF##15172667##15##]. A comprehensive understanding of DT requires the knowledge of: (1) the identities of the genes involved in the expression of the phenotype, (2) the underlying biological processes and molecular pathways, (3) the pleiotropic effects on adult traits and (4) the genetic basis of phenotypic plasticity in response to environmental variation.</p> <p>Here we report a study of DT that reveals a large number of heterochronic mutants, most of which retard the time to reach the adult stage. The fastest developing line in our survey is a mutant of <italic>Merlin</italic>, a gene encoding a protein that regulates cell proliferation in developing imaginal discs. Interestingly, <italic>Merlin </italic>is the homolog of the human tumor suppressor gene <italic>Neurofibromatosis Type 2 </italic>(NT2) involved in deregulation of cell proliferation in pathologies of the central nervous system [##REF##16872256##25##]. In addition, we provide evidence for extensive pleiotropic effects of growth control genes on adult traits. Finally, we discovered several candidate genes involved in the plastic response to temperature variation.</p>
<title>Methods</title> <title>Drosophila stocks and Developmental time assays</title> <p>We scored 179 homozygous viable <italic>P[GT1] </italic>insertion lines, contructed in a co-isogenic <italic>Canton-S B </italic>background [##REF##11156992##63##] as part of the <italic>Berkeley Drosophila Genome Project </italic>(see Availability and requirements section for URL) for DT.</p> <p>Batches of 30–40 P-element insertion lines were assessed simultaneously. To account for environmental variation in DT between batches, 8 replicated vials of the control strain (a co-isogenic P-element insertion free line with the same genetic background) was run in parallel with each batch. For each line, 300 pairs of sexually mature flies were placed in oviposition chambers for 8 hours. Eggs were allowed to hatch and batches of 30 first-instar larvae were transferred to culture vials containing a cornmeal-agar-molasses medium (4 replicates vial per line). Emerged flies from each vial were collected every 12 hours and sorted by sex. We estimated DT as the time elapsed since the transfer of first-instar larvae to the vials until adult emergence. Each vial was kept in an incubator at 25°C ± 0.5, under a 12:12 h light:dark photoperiod and at 60–70% of humidity. Thirteen lines with low viability (those that reached less than 50% of control pre-adult survival) were excluded from our survey to avoid the confusing effect of high pre-adult mortality with DT. A subset of 42 lines selected at random was also screened at 17 ± 0.5°C.</p> <title>Statistical analysis</title> <title>Quantitative genetic analysis</title> <p>Analysis of variance (ANOVA) was utilized to assess the magnitude of mutational variance for DT at each temperature. In order to include all lines tested in different batches, individual DT scores were expressed as deviations from the mean of their contemporaneous co-isogenic controls, separately for males and females. Three-way ANOVA was computed for each thermal treatment, following the mixed model: Y = μ + L + S + L × S + E, where μ is the overall mean, L is the random effect of line, S is the fixed effect of sex and E represents the error term. An ANOVA was used to assess the magnitude of mutational variance for DT induced by <italic>P[GT1] </italic>insertions. We also estimated the relative contribution of all random sources of variation (Line, Line × Sex, and error) to the total variance.</p> <p>To identify which lines were responsible of the significant line or line-by-sex interaction, phenotypic differences between <italic>P</italic>-element insertion lines and the control were tested using <italic>Dunnett </italic>contrasts for each temperature and batch. Those lines that exhibited significant DT differences relative to the control (heterochronic mutants) were considered as lines bearing an insertion in a candidate DT gene.</p> <title>Gene-by-thermal treatment analysis</title> <p>In those lines screened at both temperatures (17 and 25°C) we studied the effect of phenotypic plasticity by thermal treatments following the mixed ANOVA model, Y = μ + L + S + T + L × S + L × T + S × T + L × S × T + E, where μ is the overall mean, L is the random effect of line, S is the fixed effect of sex, T is the fixed effect of thermal treatment and E represents the error term. Since significant G × E (genotype × environment) can arise due to: (1) differences among- lines variance in separate environments (change in scale) and/or (2) deviations from unity of the cross-environment genetic correlation (<italic>rG </italic>× <italic>E </italic>&lt; 1, see below) (changes in rank order), the contribution of these two sources of variation was analyzed by means of the equation derived by Robertson [##UREF##15##64##], <italic>V</italic><sub><italic>G </italic>× <italic>E </italic></sub>= [(σ<sub>E1 </sub>- σ<sub>E2</sub>)<sup>2 </sup>+ 2σ<sub>E1</sub>σ<sub>E2</sub>(1 - <italic>rG </italic>× <italic>E</italic>)]/2, where <italic>V</italic><sub><italic>G </italic>× <italic>E </italic></sub>is the G × E variance component, <italic>rG </italic>× <italic>E </italic>is the cross-environment genetic correlation and, σ<sub>E1 </sub>and σ<sub>E2 </sub>are the square roots of the among-line variance components of the two thermal environments studied [##REF##14504242##65##]. The first term corresponds to differences in among-line variance while the second corresponds to deviations from the perfect correlation across environments (<italic>rG </italic>× <italic>E </italic>&lt; 1). The cross-environment genetic correlation (<italic>rG </italic>× <italic>E </italic>&lt; 1) is the genetic correlation of measurements in different environments and here reflects the degree in which the same genes control the expression across environments. <italic>rG </italic>× <italic>E </italic>was estimated as: <italic>rG </italic>× <italic>E </italic>= COV<sub>E1E2</sub>/σ<sub>E1 </sub>σ<sub>E2</sub>, where COV<sub>E1E2 </sub>is the covariance of lines in different environments and it was calculated as the covariance of the lines in different environments.</p> <p>Finally, we studied the phenotypic change by the effect of the <italic>P </italic>element insertion in each line showing significant differences with respect to the control in at least one thermal assay. In this case we used the following fixed ANOVA model:</p> <p>Y = μ + T + S + T × S + E, where T stands for the effect of thermal treatment, S for the effect of sex and E represents the error term. Those lines that exhibited significant differences in DT between thermal treatments (T factor) were considered as lines bearing a mutation in a candidate gene involved in the plastic response of DT to temperature variation.</p> <p>All statistical analyses were performed using STATISTICA software packages (StatSoft, Inc. 1999, 2001).</p> <title>Gene identification and functional analysis</title> <p>In order to identify the mutated genes, nucleotide sequences flanking the <italic>P</italic>-element insertion in each candidate line were used to search for homologous regions in the <italic>D. melanogaster </italic>genome. Searches of homologous sequences were performed against Release 5 of the published <italic>D. melanogaster </italic>genomic sequence using the Pubmed server. Flybase (see Availability and requirements section for URL) was used to identify target genes. Candidate DT genes were distributed in different biological processes according to their gene-ontology (GO) terms [##REF##10802651##66##]. This analysis was performed automatically with the aid of the program FatiGO+, from the Babelomics suite of bioinformatic tools, available at (see Availability and requirements section for URL)[##REF##15980512##67##]. This program distributes the genes among the GO terms in a non-exclusive manner (i.e. a given gene might be associated to more than one GO term). Using the same program we made a comparison of the distribution of GO terms between the sample of lines that developed faster and the population of lines that developed longer in order to find GO terms that were significantly over-represented in one of the samples.</p>
<title>Results</title> <title>Pattern of variation among co-isogenic <italic>P[GT1] </italic>insertion lines</title> <p>Substantial phenotypic variation among lines was observed for DT (Figure ##FIG##0##1##), wherein the most commons effect was an increase of DT relative to the control. This is not unexpected as the most likely effect of mutations on a fitness-related trait is reasonable to be deleterious, although, we found a significant positive correlation (r: 0.3, p &lt; 3 × 10<sup>5</sup>) between DT and viability [see Additional file ##SUPPL##0##1##]. This pattern revealed that the result observed of the heterochornic effect of the <italic>P[GT1] </italic>insertion lines can not be explain as a consequence of unfit flies. The most extreme phenotypes were exhibited by lines BG01543 and BG01412. The former developed 60 and 50 hours (males and females, respectively) faster and the latter 119 and 146 hours (males and females, respectively) slower than the control. We observed a significant line effect indicating mutational variance among single <italic>P</italic>-element insertion lines (Table ##TAB##0##1##). Indeed, the line effect accounts for 84% of total variance, five times higher than the error. In the ANOVA, differences between sexes were also significant, with males developing faster than females. In addition, the contribution of the line-by-sex interaction to the total phenotypic variance was not significant.</p> <title>Identification of co-isogenic <italic>P[GT1] </italic>insertion lines affecting DT and functional analysis</title> <p><italic>Dunnett </italic>analyses revealed that 107 out of 179 lines (60%) differed significantly from the control. Mean values of all significant <italic>P[GT1] </italic>insertion lines tested at 25°C are shown in Additional file ##SUPPL##1##2##. In seventy-four lines (40%) the insertion of the <italic>P</italic>-element caused an increase in DT, while in the rest of the significant lines DT was shortened as compared to the control (20%) (Figure ##FIG##0##1##). Those lines that exhibited significant DT differences relative to the control (heterochronic mutants) were considered as lines bearing an insertion in a candidate DT gene. Among fast developing lines, BG01543 and BG01902 stand out, since DT was reduced in each by more than 2 days. The <italic>P</italic>-element in line BG01902 disrupts <italic>mastermind </italic>(<italic>mam</italic>), a gene encoding a glutamine-rich nuclear protein [##REF##1701150##26##] that is a positive transcriptional regulator of the <italic>Notch </italic>signaling pathway [##REF##10221902##27##, ####REF##10545243##28##, ##REF##10801423##29##, ##REF##12050117##30####12050117##30##]. <italic>Notch </italic>is a central element in the cell signaling mechanism that controls a broad spectrum of cell fate choices during development across metazoans [##REF##14973298##31##, ####REF##16921404##32##, ##REF##17535917##33####17535917##33##]. Interestingly, in humans, abnormalities in <italic>Notch </italic>signaling have been linked to a number of diseases, such as T cell acute lymphoblastic leukemia [##REF##17363738##34##] and aortic valve disease [##REF##16025100##35##]. In line BG01543, the <italic>P</italic>-element insertion occurs in <italic>Merlin </italic>(<italic>Mer</italic>), a well-known negative regulator required for cell proliferation in developing imaginal discs of <italic>Drosophila </italic>[##REF##16872256##25##]. It is a component of the <italic>Hippo </italic>signaling pathway, which is essential for the regulation of organ size during development [##REF##16341207##36##]. <italic>Mer </italic>might control tissue growth by regulating endocytosis of membrane receptors in imaginal epithelia [##REF##17442573##37##]. Interestingly, <italic>Mer </italic>is the homolog of the human tumor suppressor gene <italic>Neurofibromatosis Type 2 </italic>(NT2) involved in the de-regulation of cell proliferation in tumors of central nervous system pathologies [##REF##16872256##25##,##REF##16204178##38##].</p> <p>Another gene involved in a cell proliferation pathway identified in our screen is <italic>forkhead box, sub-group O (foxo)</italic>. The insertion line BG01573 bearing a <italic>P</italic>-element in <italic>foxo </italic>increased DT by 24 hours in both sexes. <italic>foxo </italic>is a key regulatory component of the insulin-signaling pathway, which in <italic>Drosophila </italic>regulates the control of growth size of cells, organs, and the entire body in response to nutrient availability [##REF##12908874##39##]. Puig <italic>et al</italic>. [##REF##12893776##40##] established that <italic>foxo </italic>activates transcription in downstream as well as upstream targets of this signaling cascade by a transcriptional feedback mechanism that regulates cell growth and proliferation. It has been suggested that the insulin-signaling pathway regulates cell proliferation in imaginal discs, though the duration of the proliferation phases are controlled by Juvenile Hormone (JH) and ecdysteroids, that are themselves unaffected by the insulin-signaling pathway [##REF##16086608##41##]. Premature ecdysone release leads to the rapid development of small adults, while the delay in ecdysone release extends developmental time and produces larger adults [##REF##14597569##42##]. Recent studies have shown that the activity of the insulin-signaling pathway in the prothoracic gland modulates ecdysone release and influence both the length and the rate of larval growth [##REF##16182526##43##, ####REF##16182527##44##, ##REF##16179433##45####16179433##45##]. It has been suggested that the neuropeptide <italic>Amnesiac (Amn) </italic>participates in regulating ecdysone synthesis in the prothoracic gland of Drosophila [##REF##16182526##43##]. Interestingly, the <italic>P</italic>-element insertion mutant at <italic>Amn </italic>(line BG02286) also extended DT by more than two days in both sexes [see Additional file ##SUPPL##1##2##]. However, ecdysone can directly inhibit insulin and growth in the fat body and other peripheral tissues, an event that triggers metamorphosis. Furthermore, the ecdysone suppression of growth rate is lost in <italic>foxo </italic>mutants [##REF##16179433##45##], indicating complex cross-talk between the ecdysone and the insulin signaling pathway during <italic>Drosophila </italic>development.</p> <p>Our survey revealed extreme heterochronic mutants, particularly in the direction of increased DT [see Additional file ##SUPPL##1##2##]. For instance, the <italic>P</italic>-element insertion in the gene <italic>Karl </italic>(BG01010) extended DT by more than 100 hours in both sexes. The molecular function inferred from protein sequence suggests that <italic>Karl </italic>is a retinol-binding protein with no clear association to a biological process. In line BG00372 we registered a delay of 93 hours and 108 hours for males and females, respectively. In this case the <italic>P</italic>-element has inserted 1498 bp downstream of the 3' end of the <italic>CG1678 </italic>gene. The molecular function and the biological processes in which this gene is involved are unknown. However, there is evidence that <italic>CG1678 </italic>interacts with other genes such as <italic>And, ewg, CG1472, rl, Bsg25D, CG11275 </italic>and <italic>msb1l </italic>[##REF##14605208##46##]. Line BG01412 exhibited the most extreme phenotype in our screen, the <italic>P</italic>-element insertion was associated with an extended developmental time of 119 and 146 hours in males and females, respectively. Unfortunately, information about the nucleotide sequences flanking the <italic>P</italic>-element insertion site is not available.</p> <p>Candidate DT genes are involved in a wide range of biological processes according to their gene-ontology (GO) terms (Table ##TAB##1##2##). Interestingly, DT genes not only are involved in processes associated to organismal development but also to biosynthetic and cellular metabolic processes. We did not find significant differences between groups of candidate genes accelerating development vs those that extended development in the distribution among GO terms, suggesting that similar ontogenetic pathways may be responsible for both types of heterochronic phenotypes.</p> <title>Candidate genes for plasticity</title> <p>Growth and development of ectotherms are determined in part by their thermal environment [##UREF##9##22##,##UREF##10##23##]. In particular, temperature during ontogeny exerts a strong influence shaping the evolution of larval traits [##REF##16942614##24##,##UREF##11##47##,##UREF##12##48##]. In fact, phenotypic responses result of adaptation to different thermal environments and/or may be an unavoidable consequence of the effect of temperature on the organism's physiology during development [##UREF##13##49##]. In this context, the pattern of phenotypic effects of <italic>P</italic>-element insertion lines reared at different developmental temperatures would provide new insights in the study of phenotypic evolution of larval traits. The ANOVA showed that differences among lines and between thermal treatments were significant (Table ##TAB##2##3##). More importantly, our screen revealed that the line-by-sex and the line-by-temperature interactions were also highly significant, indicating that the behavior of each line depended on the temperature at which it was reared and the sex. However, there was a large difference in the magnitude of these genotype-by-environment interactions. The former accounts for only 1% of the total phenotypic variance (a percentage similar to that obtained in the 25°C assay) while the line-by-temperature interaction explained 52% of the variation. Moreover, note that in the general assay, the percentage of total phenotypic variance explained by differences among lines was 84%, whereas in our assays of phenotypic plasticity this percentage dropped to 30%. It may be argued that part of the effect was obscured by the high value of the line-by-temperature interaction term. This observation opens an excellent opportunity for studying the genetic basis of phenotypic plasticity of developmental time. Thus, we decided to analyze whether the significant line-by-temperature interaction can be explained by changes in magnitude of among-line variance across thermal treatments or changes in the rank order among lines, i.e. a cross-temperature genetic correlation lower than unity. Our results showed that about half of the interaction variance can be explained by a greater among-line variance observed at 17°C than 25°C and the other half by temperature-specific effects on DT of the lines (Figure ##FIG##1##2##).</p> <p>Thirty out of 42 lines tested at 17°C and 25°C showed significant differences relative to the control at one or both temperatures (Table ##TAB##3##4##). In addition, DT measured at different temperatures were positive correlated for the set of 30 heterochronic lines in both sexes (males: r: 0.61, p &lt; 0.05; females: r: 0.55, p &lt; 0.05; Figure ##FIG##2##3##). These results imply that most of the mutations that affected DT at both temperatures did so in the same phenotypic direction (either increasing or decreasing DT). Forty seven percent of the lines showed this pattern (bottom left and top right quadrants in Figure ##FIG##2##3##). Nevertheless, we also identified a set of lines in which the effect of the mutation on DT was temperature-dependent, suggesting that the mutated loci may be possible candidates for temperature plasticity genes. Indeed, forty percent of the lines affected DT in only one of the temperatures. We observed three lines showing a heterochronic phenotype only at 25°C. The line bearing an insertion in <italic>CG14478 </italic>is an example of this pattern. In contrast, nine lines showed a heterochronic phenotype only at 17°C. For example, the insertion in <italic>Imp </italic>exhibited significant differences at 17°C in both sexes, but not at 25°C. Surprisingly, four lines (13%) affected DT in opposite phenotypic directions in the two rearing temperatures. On the one hand, line BG02159, bearing a mutation in <italic>CG32666</italic>, decreased and increased DT at 17°C and 25°C, respectively. On the other hand, mutant lines for genes <italic>Nmdar</italic>, <italic>sugarless </italic>and <italic>CG11226 </italic>exhibited the opposite pattern (Figure ##FIG##3##4##).</p>
<title>Discussion</title> <title>Genetic architecture of developmental time</title> <p>The genetic architecture of a trait determines the variational properties of the phenotype, that is, its evolutionary potential. Our study revealed 116 heterochronic independent mutants involved in the expression of DT (107 lines showed a mutant phenotype at 25°C, while 9 lines were heterochronic only at 17°C). In this sense, we were capable to identify components of different kind of pathways involved in DT expression. On the one hand, some components correspond to molecular mechanisms directly involved in timekeeping processes, such as, the ecdysone and the insulin-signaling pathway [##REF##16086608##41##,##REF##15454083##50##,##REF##12852852##51##]. It has been established that environmental inputs during ontogeny influences multiple steps of the these kind of heterochronic gene pathways, such as, the nutrition status and the light: dark regimes. A paradigmatic example of this kind of environmental input is the circadian rhythm <italic>period </italic>mutants that alter DT in <italic>Drosophila </italic>[##REF##2113515##52##]. On the other hand, we also observed heterochronic phenotypes that seem to be a by-product of the disruption of a pathway that is not explicitly involved in the control of temporal process but plays a major role in organ growth This is the case of <italic>Merlin</italic>, a negative regulator of the <italic>Hippo </italic>signaling pathway, required for cell proliferation in developing imaginal discs. Interestingly, <italic>Merlin </italic>is a component of a tumour-suppressor network associated with human tumorigenesis [##REF##16341207##36##]. All in all, our results support the hypothesis raised by Moss [##REF##17320389##53##] that timing of development is control directly and indirectly by different ontogenetic pathways.</p> <p>Since the same co-isogenic lines were also screened for abdominal and sternopleural bristle number, starvation resistance and olfactory behavior [##REF##17028343##17##,##REF##12932322##54##,##REF##15126400##55##] we also examined whether DT candidate genes have pleiotropic effects on adult fitness-related traits. Several DT heterochronic lines (genes) affect starvation resistance (4 genes), bristle number (8 genes) and olfactory behavior (13 genes) [see Additional file ##SUPPL##2##3##]. Once again, <italic>Merlin </italic>seems to plays a role not only in egg-to-adult DT but also in starvation resistance and adult olfactory behavior. In BG01543 (<italic>Mer </italic>mutant), the <italic>P[GT1] </italic>insertion in exon 1 caused a reduction in the RNA expression levels in embryos, third instar larvae and adults but not in the pupal stage compared to the control [##REF##17028343##17##] indicating that anomalies caused by the disruption of <italic>Merlin </italic>on DT took place before metamorphosis. Since insects do not grow as adults, their final size is a product of growth rate and length at larval stages [##REF##17139322##56##]. Surprisingly, the mutated <italic>Merlin </italic>line did not differ significantly from the control in several morphological traits such as face width, head width, thorax length and wing size suggesting a decoupled mechanism of timing control and size development (Carreira, Mensch &amp; Fanara, unpublished results). Regarding the pleiotropic effect of <italic>Mer </italic>on starvation resistance and adult olfactory behavior, it is not easy to discern whether aberrant phenotypes are the consequence of early developmental events or are directly related to the functional disruption of the gene product in the adult. As another example, the insertion line for <italic>tramtrack </italic>(<italic>ttk</italic>), a transcription repressor gene involved in cell-fate determination [##REF##15018932##57##], increased DT by two days in both sexes and also affected starvation resistance and olfactory behavior. These interesting examples constitute clear evidence of an intricate network of pleiotropic effects of key developmental genes throughout the life-cycle of a holometabolous organism. Such pervasive pleiotropy of key genes that control development, affecting early and late fitness-related traits, show an integrated picture of the evolution of life-history traits at the molecular level.</p> <p>Regarding the effect of sex, males developed significantly faster than females, an observation that is quite striking since it is at odds with extensive literature showing that <italic>Drosophila </italic>females usually reach adulthood before males in the species studied so far [##UREF##14##58##]. Moreover, a recent study by Paranjpe <italic>et al</italic>. [##REF##15725348##59##] showed that females developed faster than males in <italic>D. melanogaster</italic>. However, some differences between our study and Paranjpe <italic>et al.'</italic>s in terms of the experimental design are worth mentioning. Paranjpe <italic>et al</italic>. found that females developed faster than males under several photoperiod regimes during ontogeny, including 12:12 L:D (the regime performed in our study). Nevertheless, we must note that different genetic backgrounds were analyzed in these two studies. Indeed, our control line, which shares the same genetic background with all mutant lines tested, did not show any sign of sexual dimorphism in DT, stressing the need to take into account the effect of the genetic background on the expression of the analyzed trait since genetic background × gene interactions are known to be quite pervasive [##REF##10924473##60##,##REF##16532011##61##]. At the same time the null input of the line-by-sex interaction to the total phenotypic variance adds a complementary aspect to the discussion of the absence of sexual dimorphism for DT. In the last few years a similar set of coisogenic insertion lines was studied for several adult traits for which the line-by-sex interaction explained significant portions of variance. In this context, our DT survey yielded the lowest value in this category, followed by sternopleural bristle number (5%), adult odor-behavior (13%), abdominal bristle number (21%) and starvation resistance (66%) [##REF##17028343##17##,##REF##12932322##54##,##REF##15126400##55##]. Among the heterochronic mutants, nine had a male-specific effect and in seven only female DT was affected. It is important to note that none of the lines showed opposite phenotypic effects across sexes, meaning that we did not find any lines that decreased DT in one sex and increased the trait in the other. In conclusion, our results indicate that there are features of the genetic architecture of DT that are highly conserved in both sexes.</p> <title>Gene-by-environment interaction</title> <p>Although the multiple interacting genes affecting complex traits can readily be dissected, how much genotype-environment interactions contribute to variation in these traits remains elusive. The quantitative dissection of developmental time reported here reveals that almost all heterochronic lines tested at different developmental temperature presented an effect with quantitative variation between temperatures, a remarkable consequence of the large gene-by-environment interaction (Table ##TAB##3##4##). The most striking examples of these developmental reaction norms are the mutant lines for genes <italic>CG11226, CG32666, Nmdar1</italic>, and <italic>sugarless </italic>(Figure ##FIG##3##4##). CG32666 gene product is a receptor signaling protein with serine/threonine kinase activity, NMDAR1 is a ionotropic glutamate receptor and SUGARLESS is a protein with UDP-glucose 6-dehrygenase activity indicating that this routine cellular metabolic processes are not only involved in developmental time but also that are sensitive to environmental variation.</p> <p>The part of the genetic diversity that has the potential to affect the phenotype, but that is not expressed under the current genotypic or environmental conditions, is referred to as "cryptic" or "hidden" genetic variation [##REF##18079017##62##]. Molecular mechanisms responsible for this particular genetic variation include epistasis and genotype-by-environment interactions. In our case, in addition to the plastic developmental reaction norms, large input to the gene-by-environment interaction refers to the bigger variance magnitude at 17°C (Figure ##FIG##1##2##). It is possible that at 25°C genetic variation was buffered in comparison to the variation expressed at 17°C, an unusual rearing temperature for these lines. Taken together, our results highlight the large potential of the <italic>Drosophila </italic>genome to changes in this relevant environmental factor, although this genetic plasticity was exhibited by these particular inbred lines, a pattern that it is not necessarily true for outbred populations. Future efforts will focus on elucidating the molecular mechanisms controlling the temperature plastic response observed. Since in nature insects are exposed to a wide range of environmental temperatures during the day and also in different seasons along the year with local-specific thermal regimes, we are also interested in studying the evolutionary forces shaping the variation in the candidate DT genes.</p>
<title>Conclusion</title> <p>The study of heterochrony has played an important role in the intersection of evolution, genetics and developmental biology since the late nineteenth century. During the last decade the concept has been revitalized with studies at cellular and molecular level. Here we identified components of several integrated time-dependent pathways affecting egg-to-adult DT and also components of pathways that are not explicitly involved in the control of temporal processes. In this sense, we identified 116 egg-to-adult heterochronic mutants, most of them developed slower compared to the control. However, we identified a set of lines that developed faster by more than two days. This is the case of <italic>Mer </italic>mutant line, a gene involved in a cell proliferation pathway. <italic>Mer </italic>and others DT candidate genes have pleitropic effects on adult traits as well. Most of developmental time candidate genes were sensitive to the rearing temperature, a fact that stresses the need to take into account the effect of environmental variation on the genetic architecture of complex traits.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Understanding the genetic architecture of ecologically relevant adaptive traits requires the contribution of developmental and evolutionary biology. The time to reach the age of reproduction is a complex life history trait commonly known as developmental time. In particular, in holometabolous insects that occupy ephemeral habitats, like fruit flies, the impact of developmental time on fitness is further exaggerated. The present work is one of the first systematic studies of the genetic basis of developmental time, in which we also evaluate the impact of environmental variation on the expression of the trait.</p> <title>Results</title> <p>We analyzed 179 co-isogenic single <italic>P[GT1]-</italic>element insertion lines of <italic>Drosophila melanogaster </italic>to identify novel genes affecting developmental time in flies reared at 25°C. Sixty percent of the lines showed a heterochronic phenotype, suggesting that a large number of genes affect this trait. Mutant lines for the genes <italic>Merlin </italic>and <italic>Karl </italic>showed the most extreme phenotypes exhibiting a developmental time reduction and increase, respectively, of over 2 days and 4 days relative to the control (a co-isogenic <italic>P</italic>-element insertion free line). In addition, a subset of 42 lines selected at random from the initial set of 179 lines was screened at 17°C. Interestingly, the gene-by-environment interaction accounted for 52% of total phenotypic variance. Plastic reaction norms were found for a large number of developmental time candidate genes.</p> <title>Conclusion</title> <p>We identified components of several integrated time-dependent pathways affecting egg-to-adult developmental time in <italic>Drosophila</italic>. At the same time, we also show that many heterochronic phenotypes may arise from changes in genes involved in several developmental mechanisms that do not explicitly control the timing of specific events. We also demonstrate that many developmental time genes have pleiotropic effects on several adult traits and that the action of most of them is sensitive to temperature during development. Taken together, our results stress the need to take into account the effect of environmental variation and the dynamics of gene interactions on the genetic architecture of this complex life-history trait.</p>
<title>Availability and requirements</title> <p>Berkeley Drosophila Genome Project<italic>: </italic><ext-link ext-link-type="uri" xlink:href="http://flypush.imgen.bcm.tmc.edu/pscreen/"/></p> <p>Flybase: <ext-link ext-link-type="uri" xlink:href="http://flybase.bio.indiana.edu/"/></p> <p>FatiGO+: <ext-link ext-link-type="uri" xlink:href="http://www.fatigo.org/"/></p> <title>Abbreviations</title> <p>DT: developmental time; ANOVA: analysis of variance; GO: gene-ontology.</p> <title>Authors' contributions</title> <p>JM and JJF conceived the study. JM, NL, VPC and AM performed the developmental time assays. JM analyzed the data and wrote the first draft of this manuscript. JM, EH and JJF wrote the final version. NL and VPC helped to draft the final version of this manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>The authors wish to thank T.F.C. Mackay and R. R. H. Anholt for comments and suggestions that helped to improve this article. We also want to thank A. Yamamoto for technical support, two anonymous reviewers for their constructive criticisms and all colleagues at the Laboratorio de Evolución. This work was supported by grants of the Universidad de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) and Agencia Nacional de Promoción Científica y Técnica. J. Mensch, N. Lavagnino and V.P. Carreira are recipients of post-graduate scholarships of CONICET. E.H. and J.J. Fanara are members of Carrera del Investigador Científico (CONICET).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Frequency distribution of mutational effects on developmental time of 179 <italic>P[GT1] </italic>insertion lines. </bold>Grey bar represents lines showing similar effects relative to the control line (40%); green and blue bars represent lines showing a reduction (20%) and an increase (40%) in DT, respectively. Values in x-axis show the limits of each DT intervals.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Reaction norms of 45 <italic>P[GT1] </italic>insertion lines tested at 17°C and 25°C.</bold> Values are shown as the deviation of the insert line mean from contemporaneous control line.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Correlation of 30 <italic>P[GT1] </italic>heterochronic mutants tested at 17°C and 25°C. </bold>Values are shown as the deviation of the insert line mean from contemporaneous control line for each sex separately; males (solid circles) and females (open circles).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Heterochronic mutants having opposite phenotypic effects between thermal treatments. </bold>17°C assay is represent with open bars and 25°C assay with red bars. Males are represented with striped bars and females with open bars.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Analysis of variance of mutational effects on developmental time.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Source</td><td align="center">d.f.</td><td align="center"><italic>F</italic></td><td align="center"><italic>P</italic></td><td align="center">σ<sup>2</sup>(%)</td></tr></thead><tbody><tr><td align="center">Line</td><td align="center">178</td><td align="center">54.68</td><td align="center">&lt;0.0001</td><td align="center">84</td></tr><tr><td align="center">Sex</td><td align="center">1</td><td align="center">46.87</td><td align="center">&lt;0.0001</td><td align="center">Fixed</td></tr><tr><td align="center">Line × Sex</td><td align="center">178</td><td align="center">0.76</td><td align="center">0.99</td><td align="center">0</td></tr><tr><td align="center">Error</td><td align="center">1066</td><td/><td/><td align="center">16</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Distribution of candidate DT genes among biological process GO terms.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Biological process GO terms</td><td align="center">% of genes</td></tr></thead><tbody><tr><td align="left">Cellular metabolic process</td><td align="center">57.9</td></tr><tr><td align="left">Multicellular organismal development</td><td align="center">39.5</td></tr><tr><td align="left">Macromolecule metabolic process</td><td align="center">36.8</td></tr><tr><td align="left">Anatomical structure development</td><td align="center">34.2</td></tr><tr><td align="left">Regulation of biological process</td><td align="center">31.6</td></tr><tr><td align="left">Cellular developmental process</td><td align="center">31.6</td></tr><tr><td align="left">Cell communication</td><td align="center">29.0</td></tr><tr><td align="left">Biosynthetic process</td><td align="center">23.7</td></tr><tr><td align="left">Cellular component organization and biogenesis</td><td align="center">18.4</td></tr><tr><td align="left">Sexual reproduction</td><td align="center">18.4</td></tr><tr><td align="left">Establishment of localization</td><td align="center">18.4</td></tr><tr><td align="left">Nitrogen compound metabolic process</td><td align="center">10.5</td></tr><tr><td align="left">Behavior</td><td align="center">10.5</td></tr><tr><td align="left">Response to stress</td><td align="center">10.5</td></tr><tr><td align="left">Cell proliferation</td><td align="center">10.5</td></tr><tr><td align="left">Cell adhesion</td><td align="center">10.5</td></tr><tr><td align="left">Cell cycle</td><td align="center">7.9</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Analysis of variance of mutational effects on developmental time among lines tested at different thermal treatments.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Source</td><td/><td align="center">d.f.</td><td align="center"><italic>F</italic></td><td align="center"><italic>P</italic></td><td align="center">σ<sup>2</sup>(%)</td></tr></thead><tbody><tr><td align="center">Line</td><td/><td align="center">41</td><td align="center">2.11</td><td align="center">0.0085</td><td align="center">30</td></tr><tr><td align="center">Sex</td><td/><td align="center">1</td><td align="center">8.82</td><td align="center">0.0047</td><td align="center">Fixed</td></tr><tr><td align="center">Thermal Treatment</td><td/><td align="center">1</td><td align="center">4.56</td><td align="center">0.03</td><td align="center">Fixed</td></tr><tr><td align="center">Line × Sex</td><td/><td align="center">41</td><td align="center">1.73</td><td align="center">0.041</td><td align="center">1</td></tr><tr><td align="center">Line × Thermal Treatment</td><td/><td align="center">41</td><td align="center">41.4</td><td align="center">&lt;0.0001</td><td align="center">52</td></tr><tr><td align="center">Sex × Termal Treatment</td><td/><td align="center">1</td><td align="center">9.26</td><td align="center">0.0038</td><td align="center">Fixed</td></tr><tr><td align="center">Line × Sex × Thermal Treatment</td><td/><td align="center">41</td><td align="center">0.56</td><td align="center">0.9886</td><td align="center">0</td></tr><tr><td align="center">Error</td><td/><td align="center">479</td><td/><td/><td align="center">17</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Effects of significant <italic>P[GT1] </italic>insertion lines screened at different thermal treatments.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Line</td><td align="center">Gene</td><td align="center">Males 17°C</td><td align="center">Females 17°C</td><td align="center">Males 25°C</td><td align="center">Females 25°C</td><td align="center">S</td><td align="center">T</td><td align="center">S × T</td></tr></thead><tbody><tr><td align="left">BG00369</td><td align="center"><italic>CG13334</italic></td><td align="center">-34.25*</td><td align="center">-44.73**</td><td align="center">-44.50***</td><td align="center">-56.99**</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG00373</td><td align="center"><italic>CG11226</italic></td><td align="center">84.18***</td><td align="center">72.11***</td><td align="center">16.22*</td><td align="center">23.52***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG00386</td><td align="center"><italic>Nmdar1</italic></td><td align="center">NS</td><td align="center">54.81**</td><td align="center">-30.73***</td><td align="center">-36.13***</td><td align="center">N</td><td align="center">***</td><td align="center">*</td></tr><tr><td align="left">BG00737</td><td align="center"><italic>Hsp27</italic></td><td align="center">-64.24***</td><td align="center">-55.94***</td><td align="center">NS</td><td align="center">NS</td><td align="center">**</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG00930</td><td align="center">no sequence</td><td align="center">45.08***</td><td align="center">NS</td><td align="center">NS</td><td align="center">NS</td><td align="center">**</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01028</td><td align="center"><italic>CG33260</italic></td><td align="center">-35.5*</td><td align="center">-41.24**</td><td align="center">NS</td><td align="center">-21.81**</td><td align="center">N</td><td align="center">**</td><td align="center">N</td></tr><tr><td align="left">BG01214</td><td align="center"><italic>CG10064</italic></td><td align="center">47.86***</td><td align="center">NS</td><td align="center">-16.82*</td><td align="center">-28.59***</td><td align="center">***</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01339</td><td align="center"><italic>clt</italic></td><td align="center">92.63***</td><td align="center">93.64***</td><td align="center">39.14***</td><td align="center">29.46***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01354</td><td align="center"><italic>CG30492</italic></td><td align="center">NS</td><td align="center">41.53**</td><td align="center">NS</td><td align="center">NS</td><td align="center">**</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01488</td><td align="center"><italic>msi</italic></td><td align="center">NS</td><td align="center">-41.18**</td><td align="center">-19.40**</td><td align="center">-26.76***</td><td align="center">N</td><td align="center">*</td><td align="center">N</td></tr><tr><td align="left">BG01672</td><td align="center"><italic>CG14591</italic></td><td align="center">-85.53***</td><td align="center">-48.45**</td><td align="center">-27.40**</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01716</td><td align="center"><italic>paps</italic></td><td align="center">-38.8*</td><td align="center">-39.38*</td><td align="center">NS</td><td align="center">-22.74**</td><td align="center">*</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01726</td><td align="center"><italic>CG11382</italic></td><td align="center">-39.82**</td><td align="center">-38.27*</td><td align="center">NS</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01735</td><td align="center"><italic>bib</italic></td><td align="center">NS</td><td align="center">NS</td><td align="center">NS</td><td align="center">-19.91**</td><td align="center">N</td><td align="center">N</td><td align="center">*</td></tr><tr><td align="left">BG01763</td><td align="center"><italic>CG33960</italic></td><td align="center">173***</td><td align="center">166.28***</td><td align="center">34.41***</td><td align="center">NS</td><td align="center">N</td><td align="center">**</td><td align="center">N</td></tr><tr><td align="left">BG01780</td><td align="center"><italic>CG11226</italic></td><td align="center">42.8**</td><td align="center">49.41***</td><td align="center">-20.18*</td><td align="center">-32.31***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01822</td><td align="center"><italic>Imp</italic></td><td align="center">-47.95**</td><td align="center">-52.38**</td><td align="center">NS</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01902</td><td align="center"><italic>mam</italic></td><td align="center">-59.92***</td><td align="center">-62.62***</td><td align="center">-47.64***</td><td align="center">-53.05***</td><td align="center">N</td><td align="center">*</td><td align="center">N</td></tr><tr><td align="left">BG01912</td><td align="center"><italic>pxb</italic></td><td align="center">-79.85***</td><td align="center">-97.30***</td><td align="center">-33.10***</td><td align="center">-45.92***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG01990</td><td align="center"><italic>CG30492</italic></td><td align="center">-43.82**</td><td align="center">-43.86**</td><td align="center">-35.87***</td><td align="center">-46.05***</td><td align="center">N</td><td align="center">N</td><td align="center">N</td></tr><tr><td align="left">BG02042</td><td align="center"><italic>eas</italic></td><td align="center">47.78**</td><td align="center">NS</td><td align="center">NS</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02088</td><td align="center"><italic>CG15309</italic></td><td align="center">49.55**</td><td align="center">46.28**</td><td align="center">NS</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02102</td><td align="center"><italic>CG13434</italic></td><td align="center">-61.79***</td><td align="center">-70.34***</td><td align="center">-16.73*</td><td align="center">-17.89*</td><td align="center">*</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02157</td><td align="center"><italic>CG8177</italic></td><td align="center">90.67***</td><td align="center">76.60***</td><td align="center">NS</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02159</td><td align="center"><italic>CG32666</italic></td><td align="center">-85.12***</td><td align="center">-96.73***</td><td align="center">29.51***</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">*</td></tr><tr><td align="left">BG02239</td><td align="center"><italic>CG11550</italic></td><td align="center">63.36***</td><td align="center">60.57***</td><td align="center">29.40***</td><td align="center">22.72***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02462</td><td align="center"><italic>CG6301</italic></td><td align="center">50.77***</td><td align="center">38.37*</td><td align="center">NS</td><td align="center">NS</td><td align="center">*</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02690</td><td align="center"><italic>CG14478</italic></td><td align="center">NS</td><td align="center">NS</td><td align="center">-32.02***</td><td align="center">-37.68***</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02747</td><td align="center"><italic>rut</italic></td><td align="center">136.73***</td><td align="center">129.40***</td><td align="center">28.00***</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr><tr><td align="left">BG02823</td><td align="center"><italic>scyl</italic></td><td align="center">NS</td><td align="center">NS</td><td align="center">17.44*</td><td align="center">NS</td><td align="center">N</td><td align="center">***</td><td align="center">N</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Correlation between viability (VT) and developmental time (DT) for all <italic>P[GT1] </italic>insertion lines analyzed. VT is presented in angular transformation (arcsine of the square root of the proportion of survival) and DT expressed in hours. Solid line represents the correlation plot and dashed lines the 95% confidence intervals (r = 0.305).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Mean values of significant <italic>P[GT1] </italic>insertion lines on DT (25°C). Values are shown as the deviation of the insertion line mean from contemporaneous control for each sex separately.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Genetic information of candidate DT genes including gene names, <italic>P[GT1] </italic>insertion sites, cytologial positions, biological process gene-ontologies and pleiotropy.</p></caption></supplementary-material>
<table-wrap-foot><p>d.f.: degree of freedom; σ<sup>2</sup>: component variance.</p></table-wrap-foot> <table-wrap-foot><p>Genes are distributed in a non-exclusive manner (i.e. a given gene might be associated to more than one GO term, see text for more details). The percentage of genes related to each GO term is shown.</p></table-wrap-foot> <table-wrap-foot><p>d.f.: degree of freedom;: σ<sup>2 </sup>component variance.</p></table-wrap-foot> <table-wrap-foot><p>Values are shown as the deviation of the insert line mean from contemporaneous control line. Asterisks in columns S, T and S × T represent the quantitative effect of Sex, Thermal treatment and Sex-by-Thermal treatment in response to environmental variation (see materials and methods for detail statistical treatment). *p &lt; 0.05; **p &lt; 0.01; ***p &lt; 0.001; NS: no significant.</p></table-wrap-foot>
<graphic xlink:href="1471-213X-8-78-1"/> <graphic xlink:href="1471-213X-8-78-2"/> <graphic xlink:href="1471-213X-8-78-3"/> <graphic xlink:href="1471-213X-8-78-4"/>
<media xlink:href="1471-213X-8-78-S1.bmp" mimetype="image" mime-subtype="bmp"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-213X-8-78-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-213X-8-78-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Stearns"], "given-names": ["SC"], "source": ["The evolution of Life Histories"], "year": ["1992"], "publisher-name": ["Oxford University Press"]}, {"surname": ["Gould"], "given-names": ["SJ"], "source": ["Ontogeny and phylogeny"], "year": ["1977"], "publisher-name": ["Harvard University Press"]}, {"surname": ["Carrol"], "given-names": ["S"], "source": ["Endless forms most beautiful"], "year": ["2005"], "publisher-name": ["WW Norton & Company"]}, {"surname": ["Nunney"], "given-names": ["L"], "article-title": ["The response to selection for fast larval development in Drosophila melanogaster and its effect on adult weight: an example of a fitness trade-off"], "source": ["Evolution"], "year": ["1996"], "volume": ["50"], "fpage": ["1193"], "lpage": ["1204"], "pub-id": ["10.2307/2410660"]}, {"surname": ["Zwaan", "Bijsma", "Hoekstra"], "given-names": ["BJ", "R", "RF"], "article-title": ["Artificial selection for development time in "], "italic": ["Drosophila melanogaster "], "source": ["Evolution"], "year": ["1995"], "volume": ["49"], "fpage": ["635"], "lpage": ["648"], "pub-id": ["10.2307/2410317"]}, {"surname": ["Chippindale", "Alipaz", "Chen", "Rose"], "given-names": ["AK", "JA", "H", "M"], "article-title": ["Experimental evolution of accelerated development in "], "italic": ["Drosophila"], "source": ["Evolution"], "year": ["1997"], "volume": ["51"], "fpage": ["1536"], "lpage": ["1551"], "pub-id": ["10.2307/2411206"]}, {"surname": ["Chippindale", "Alipaz", "Rose", "Rose MR, Passananti HB, Matos M"], "given-names": ["AK", "JA", "MR"], "article-title": ["Experimental evolution of accelerated development in "], "italic": ["Drosophila"], "source": ["Methuselah Flies: A case study in the evolution of aging"], "year": ["2004"], "publisher-name": ["Singapore: World Scientific"], "fpage": ["413"], "lpage": ["435"]}, {"surname": ["West-Eberhard"], "given-names": ["MJ"], "source": ["Developmental plasticity and evolution"], "year": ["2003"], "publisher-name": ["Oxford University Press"]}, {"surname": ["Schlichting", "Pigliucci"], "given-names": ["C", "M"], "source": ["Phenotypic Evolution: a reaction norm perspective"], "year": ["1998"], "publisher-name": ["Sinauer Associates"]}, {"surname": ["Davidowitz", "Nijhout"], "given-names": ["G", "F"], "article-title": ["The physiological basis of reaction norms: the interaction among growth rate, the duration of growth and body size"], "source": ["Integr Comp Biol"], "year": ["2004"], "volume": ["44"], "fpage": ["443"], "lpage": ["449"], "pub-id": ["10.1093/icb/44.6.443"]}, {"surname": ["Angilletta", "Steury", "Sears"], "given-names": ["M", "T", "M"], "article-title": ["Temperature, growth rate, and body size in ectotherms: fitting pieces of a life-history puzzled"], "source": ["Integr Comp Biol"], "year": ["2004"], "volume": ["44"], "fpage": ["498"], "lpage": ["509"], "pub-id": ["10.1093/icb/44.6.498"]}, {"surname": ["James", "Partridge"], "given-names": ["AC", "L"], "article-title": ["Thermal evolution of the rate of larval development in "], "italic": ["Drosophila melanogaster "], "source": ["J Evol Biol"], "year": ["1995"], "volume": ["8"], "fpage": ["315"], "lpage": ["330"], "pub-id": ["10.1046/j.1420-9101.1995.8030315.x"]}, {"surname": ["Folguera", "Ceballos", "Spezzi", "Fanara", "Hasson"], "given-names": ["G", "S", "L", "JJ", "E"], "article-title": ["Clinal variation in developmental time and viability, and the response to thermal treatments in two species of "], "italic": ["Drososphila"], "source": ["Biol J Linn Soc"]}, {"surname": ["Stillwell", "Fox"], "given-names": ["RC", "CW"], "article-title": ["Complex patterns of phenotypic plasticity: interactive effects of temperature during rearing and oviposition"], "source": ["Ecology"], "year": ["2005"], "volume": ["86"], "fpage": ["924"], "lpage": ["934"], "pub-id": ["10.1890/04-0547"]}, {"surname": ["Markow", "O'Grady"], "given-names": ["T", "P"], "italic": ["Drosophila"], "year": ["2005"], "publisher-name": ["Academic Press"]}, {"surname": ["Robertson"], "given-names": ["A"], "article-title": ["The sampling variance of the genetic correlation coefficient"], "source": ["Biometrics"], "year": ["1959"], "volume": ["15"], "fpage": ["469"], "lpage": ["485"], "pub-id": ["10.2307/2527750"]}]
BMC Dev Biol. 2008 Aug 8; 8:78