Effects of small Hspgenes on developmental stability and microenvironmental canalization
© Takahashi et al; licensee BioMed Central Ltd. 2010
Received: 23 June 2010
Accepted: 16 September 2010
Published: 16 September 2010
Progression of development has to be insulated from the damaging impacts of environmental and genetic perturbations to produce highly predictable phenotypes. Molecular chaperones, such as the heat shock proteins (HSPs), are known to buffer various environmental stresses, and are deeply involved in protein homeostasis. These characteristics of HSPs imply that they might affect developmental buffering and canalization.
We examined the role of nine Hsp genes using the GAL4/UAS-RNAi system on phenotypic variation of various morphological traits in Drosophila melanogaster. The stability of bristle number, wing size and wing shape was characterized through fluctuating asymmetry (FA) and the coefficient of variation (CV), or among-individual variation. Progeny of the GAL4/Hsp-RNAi crosses tended to have reduced trait means for both wing size and wing shape. Transcriptional knockdown of Hsp67Bc and Hsp22 significantly increased FA of bristle number, while knockdown of Hsp67Ba significantly increased FA and among-individual variation of wing shape but only in males. Suppression of Hsp67Bb expression significantly increased among-individual variation of bristle number. The knockdown of gene expression was confirmed for Hsp67Ba, Hsp67Bc, Hsp22, and Hsp67Bb. Correlation between FA and CV or among-individual variation of each trait is weak and not significant except for the case of male wing shape.
Four small Hsp genes (Hsp22, Hsp67Ba, Hsp67Bb and Hsp67Bc) showed involvement in the processes of morphogenesis and developmental stability. Due to possible different functions in terms of developmental buffering of these small Hsps, phenotypic stability of an organism is probably maintained by multiple mechanisms triggered by different environmental and genetic stresses on different traits. This novel finding may lead to a better understanding of non-Hsp90 molecular mechanisms controlling variability in morphological traits.
Progression of development has to be insulated from the damaging impacts of environmental and genetic perturbations to produce highly predictable phenotypes. Waddington  suggested a conceptual mechanism called canalization that buffers developmental processes from environmental and genetic perturbations and therefore helps to produce constant phenotypes. Molecular chaperones, such as the heat shock proteins (HSPs), are known to buffer various environmental stresses, and are deeply involved in protein homeostasis . Those characteristics of HSPs imply that they have potential to be candidates for developmental buffering and canalization.
Inhibition of Hsp90, one of the molecular chaperones, has been found to increase phenotypic diversity in various organisms such as Drosophila, Arabidopsis, and zebrafish [3–5]. It suggests that Hsp90 buffers developmental perturbations on morphological traits in these species. However, Milton et al.  and Debat et al.  found that a reduction of HSP90 activity did not affect phenotypic variation. Based on the mixed results, Debat et al.  suggested that Hsp90 is one of the multiple factors that participate in the developmental buffering of morphological traits rather than the only controlling factor.
Hsp70 is one of the most well-studied stress response genes, and is inducible by thermal and nutritional stresses and inbreeding in D. melanogaster [8–12]. Previous research suggests that Hsp70 contributes to stabilization of developmental processes, although results are somewhat inconsistent. Roberts and Feder  showed that increased copy number of Hsp70 significantly reduced developmental abnormalities, while Williams et al.  observed the opposite effect.
Other molecular chaperones such as Hsp22, Hsp67, Hsp68, and Hsc70 are also known to respond to environmental stresses [15–19]. It has been suggested that they contribute to thermotolerance , and some of them (Hsp22, Hsp68, and Hsp70) affect longevity . Although details of the chaperone activity and the molecular mechanism of the Hsp-mediated stress resistance are largely unknown, it is possible that these genes may also affect developmental processes.
In D. melanogaster, expression of Hsp genes is not only rapidly up-regulated under environmental stresses, but is also regulated during normal development. In the absence of environmental stress, Hsp68 and Hsp70 mRNA are expressed at very low levels in most developmental stages, but they are at higher concentration in pupae . Hsp22, Hsp67Ba, and Hsp67Bc also reach a maximum level of expression in the early pupal stage in the absence of environmental stress [21, 22]. Hsp67Bb mRNA is detected during all the larval stages to early pupal stage in a tissue-specific manner . The higher expression of the Hsps in embryos and pupae coincides with major developmental events. At the embryonic stage the body plan is being formed, while key metamorphosis processes occur during the pupal stage. Hsp activities during these periods occur at a time when meticulous spatial expression patterns develop. Such temporal regulation of expression suggests that these Hsp genes might help stabilize developmental processes at critical times.
This study aims to investigate the potential role of a subset of non-Hsp90 heat shock protein genes in phenotypic variability using RNA interference. If Hsp genes are involved in buffering phenotypic variability, suppression of their expression might result in a decrease in developmental stability and canalization. Here, we define developmental stability as a set of mechanisms buffering developmental variation among replicated or symmetrical organs within a single organism, and microenvironmental canalization as an increase in the phenotypic variance of a morphological trait, following Debat and David . We infer developmental stability by measuring fluctuating asymmetry (FA) and the coefficient of variation (CV) or among-individual variation of several bristle traits and wing traits. Using the actin-GAL4/UAS-Hsp-RNAi system, we found that transcriptional suppression of Hsp22 and Hsp67Bc affected FA of bristle traits. We also detected a significant positive correlation between FA and among-individual variation of wing shape under Hsp67Ba knockdown. Furthermore, Hsp67Bb may influence among-individual variation of bristle traits. Our findings suggest developmental roles for non-Hsp90 stress proteins in controlling the expression of morphological variability.
RNAi strains and their specificity score (s19), target and off-target genes.
Pkcδ, CG14656, CG15725, CG15803, ASPP, CG30377, Samuel, shep, CG32541, Hs3st-A, CG33988, vir, Smr, bun, CG5697, CG5794, JIL-1, shn, Rgl, brk
Hsp70Ba, Hsp70Bb, Hsp70Bbb, Hsp70Bc
Hsp68, Hsp70Aa, Hsp70Ab, Hsc70-2, Hsc70-1
Hsp70Ba, Hsp70Bb, Hsp70Bbb, Hsp70Bc
Hsp70Aa, Hsp70Ab, Hsc70-1, Hsc70-2, CG14786
One hundred eggs were collected from each cross, and placed into a 42 ml glass vial with 10 ml of the fly medium. Fly food was made according to the Bloomington stock centre, where corn syrup was replaced with dextrose medium http://flystocks.bio.indiana.edu/Fly_Work/media-recipes/bloomfood.htm. Vials were then put into 27°C cabinets. Emerging adults were collected every day and preserved in 70% ethanol for morphological measurements. Five replicate vials were set up for each cross, and we measured three females and males from each vial.
Morphological traits and shape analysis
To visualize the effect of RNAi on WS, we performed thin-plate spline analysis. Thin-plate splines indicate the differences in two configurations of landmarks as a continuous deformation using regression functions in which corresponding landmarks are matched between configurations to minimize the bending energy . Bending energy is the energy required to bend an infinitely thin metal plate over one set of landmarks so that the height over each landmark is equal to the coordinates of the corresponding landmark in the other configuration . We compared the mean wing shape of the GAL4/+ to each of the GAL4/UAS-RNAi crosses and then exaggerated their difference by 30-fold for graphical display.
For bristle traits and CS, FA was evaluated as |L-R|/(L+R)/2, where L indicates a trait value on the left side and R on the right side of the body. To evaluate FA of bristle traits collectively, we used a composite index of FA (CFA) proposed by Leung et al. (2000). To calculate this index (CFA2 in Leung et al. ), individual FA value was divided by the average FA of a given trait in the population of interest so that all traits contribute equally to CFA measure. The FA values were then summed across traits for each individual, creating a composite FA score for each individual. CFA2 is one of the best CFAs according to the simulation by Leung et al. (2000). We used all bristle traits to calculate CFA2. For WS, we used a univariate measure of FA devised by Klingenberg and Monteiro . This measure is based on the idea of one-sample standard distance [34, 35], and is equivalent to the one-sample version of the Mahalanobis distance , automatically providing a correction for directional asymmetry .
To evaluate the effect of RNAi knockdowns on microenvironmental canalization, we calculated among-individual variation of bristle traits, CS and WS. As for CS, variance and mean were calculated for each replicate vial to obtain the coefficient of variation (CV). As for bristle traits, we calculated the trace of the total covariance matrix for among-individual variation from all the components of PCA, and used it as an index of among-individual variation. As for WS, we performed the same analysis based on the relative warp, and used the trace as an index of among-individual variation. All these analyses were performed for females and males separately.
Repeatability and measurement error of wing traits
To evaluate the measurement error in landmark acquisition based on two repeated measurements for each landmark, repeatability (R) was calculated for each landmark coordinate and CS . The repeatability measure determines the proportion of variance due to variation between individuals where zero indicates that all variance is attributable to variance within individuals (100% measurement error) whilst one indicates all variance is found between individuals (0% measurement error).
To assess the relative amounts of directional asymmetry (DA), FA, and measurement error in wing shape variation, we employed Procrustes ANOVA  with degrees of freedom under the isotropic model . In this analysis, we included individual, and side and their interaction terms, and added sums of squares across landmarks and coordinates, assuming equal and isotropic variation at each landmark.
To evaluate the effect of RNAi knockdown on FA and CV or among-individual variation, we performed Dunnett tests with bristle CFA, CS FA, WS FA, bristle among-individual variation, CS CV or WS among-individual variation as a dependent variable and RNAi strain as an independent variable. We applied Bonferroni correction to account for multiple Dunnett tests performed for each trait for each sex. The effect of RNAi knockdown on mean trait values was tested with the same analysis using Bristle PC1, CS or WS as dependent variables. We checked normality of the distributions using the Kolmogorov-Smirnov test. No significant deviations from normal distribution were detected for any measure.
Quantitative real time PCR
We measured relative expression levels of Hsp genes in the early pupal stage of the GAL4/RNAi individuals and the control (GAL4/+) by performing quantitative real time PCR (RTPCR). Test larvae were produced by mating males from RNAi line with females of the GAL4 strain to drive RNAi, while control samples were obtained by crossing males of the control stain to females of the GAL4 strain. One hundred eggs were collected from each cross and placed on fly media as described earlier. Eggs were reared at 27°C, and larvae were sampled at a 70% of average pre-adult developmental time. Female and male pupae were not distinguished. We set up three biological replications and sampled 15 to 30 individuals for each replication.
Total RNA was extracted using TRIzol® Reagent (Invitrogen Cat. No. 15596-026) for the GAL4/UAS-RNAi individuals using strains 26416, 33207, 36641, 43632, 47145, 49795, and 49796. Five μg of total RNA was treated with 0.73 units of RQ1 RNase-Free DNase (Promega Cat. No. M6101) before converting to first strand cDNA using the SuperScript® III First-Strand Synthesis SuperMix system (Invitrogen Cat. No. 18080-400). Oligo-dT was used to prime the reverse transcription. As for the GAL4/UAS-RNAi individuals using strain 21608, total RNA was extracted using the SV Total RNA Isolation System (Promega Cat. No. Z3100), and then converted to first strand cDNA using SuperScript III First Strand Synthesis System for RT-PCR (Cat. No. 18080-051) with Oligo-dT primer. cDNA was diluted 20 to 50 times in water.
Real time PCR was performed on the Roche LightCycler® 480. Five replications were set up per gene per cDNA sample. For the GAL4/UAS-RNAi individuals using strains 26416, 33207, 36641, 43632, 47145, 49795, 49796, we used 10 μl reaction containing the following components: 1 μl of cDNA, 1 μl of Immolase 10× buffer, 0.4 μl of MgCl2 (25 mM), 0.8 μl of dNTP (2 mM), 4 μl of forward and reverse primer mix (1 μM each), 0.25 μl of the LightCycler® 480 High Resolution Melting Master Mix (Roche Cat. No. 04909631001), 0.01 μl of Immolase DNA polymerase (Bioline Cat. No. BIO-21047; 5 units per μl), and 2.54 μl of water. The thermocycling conditions were as follows: 95°C for 10 min, 50 cycles of 95°C for 5 sec, 58°C for 15 sec, and 72°C for 15 sec. Fluorescence signal was recorded at the end of each 72°C elongation phase. For GAL4/UAS-RNAi individuals using strain 21806, we used 20 μl reaction that consisted of 2 μl of cDNA, 4 μl of forward and reverse primer mix (5 μM), 10 μl of LightCycler® 480 SYBR Green I Master (Cat. No. 4707516), and 4 μl of water. The thermo cycling conditions were as follows: 95°C for 5 min, 45 cycles of 95°C for 10 sec, 60°C for 20 sec, and 72°C for 20 sec.
The Crossing point (Cp) estimates were acquired from the LightCycler® 480 using the Absolute Quantification Module. The conventional 2-ΔΔCt method was used to estimate relative gene expression. The mean Cp value was calculated for each cDNA sample based on five technical replicates. The relative expression of target gene was normalized with a housekeeping gene, ribosomal protein L11 (RpL11). The PCR primers for RpL11 for GAL4/21806 were described in Bogwitz et al. . Mean normalized expression of the target gene was estimated based on three biological replications in each treatment. One-way ANOVA was performed to compare between two treatments (GAL4/UAS versus GAL4/+). The primers for real-time PCR are listed in Additional file: Appendix 1.
Repeatability and measurement error
Procrustes ANOVA for the wing landmarks.
Individual × Side
Individual × Side
RNAi effect on trait means
Mean score of Bristle PC1(OR), centroid size (CS), and wing shape (WS) of RNAi driven offspring from crosses between GAL4 and RNAi strains with standard deviation in parenthesis.
21806 (Hsp67Ba )
26416 (Hsp67Bc )
33207 (Hsp70B s)
36641 (Hsp70B s)
43632 (Hsp22 )
47145 (Hsp68 )
49795 (Hsp67Bb )
49796 (Hsp67Bb )
The phenotypic effect of each cross seemed to be trait specific except for GAL4/26416 (targeting Hsp67Bc) in females and GAL4/47145 (targeting Hsp68) in males. For example, the GAL4/21086 (targeting Hsp67Ba) affected CS and WS but there was no significant effect on bristle PC1. The GAL4/36641 (targeting Hsp68) also affected WS but not CS or bristle PC1 in either sex.
RNAi effect on phenotypic variation and expression analysis
Correlation coefficient of the pair of FA and CV of centroid size (CS), and among-individual variation of bristle and wing shape (WS)
Bristle CFA2 vs. Bristle among-individual variation
CS FA vs. CS CV
WS FA vs. WS among-individual variation
In this study, we tested whether dampened Hsp expression affects developmental stability and canalization in D. melanogaster. We attempted to knock down individual Hsp genes using the GAL4/UAS-RNAi system and measure various aspects of wing and bristle morphology. Evidence of gene suppression and phenotypic consequences on developmental stability and canalization were obtained for Hsp22, Hsp67Ba, Hsp67Bb, and Hsp67Bc knockdown, whereas results from RNAi targeting of Hsp68 and Hsp70Bb were inconclusive. We observed that the phenotypic outcomes on developmental stability and canalization of the four small Hsp knockdown were sex- and trait-dependent. It is known that the random integration of P-element-UAS constructs in the genome causes variability in the expression level of double-stranded RNA , and this might affect RNAi efficiency in this study. However, given that all RNAi strains used in the current study have an isogenic background except for the P-element insertion, the phenotypic effects are likely to be caused by reduced expression of the target and/or off-target genes.
Our results indicate that Hsp22 affects bristle but not wing shape asymmetry. Drosophila bristles or macrochaetae are important peripheral sensory organs; their spatial organization is likely to be a canalized trait. In a genetic modifier screen, Pena-Rangel et al.  found that alteration of Hsp22 expression (line EP(3)3247) could suppress the pannier (pnr) mutant phenotype, suggesting a genetic interaction between Hsp22 and mutant expression. PNR is a zinc finger transcription factor which activates wingless and the achaete-scute genes, which are crucial for bristle determination [43, 44]. Together with other patterning regulators, PNR helps define the permissive (bristle mother cells) and restrictive regions for subsequent formation of the bristles . The function of Hsp22 in buffering the phenotypic effect of a pnr mutation supports the idea that Hsp22 is involved in developmental stability of bristle traits. By examining the phenotypic effect of pnr mutant phenotype under overexpression or suppression of the expression of Hsp22 in a developmental stage specific manner, a possible link between Hsp22 and pnr could be investigated further.
The development of Drosophila bristles might parallel that of the nervous system. Held  suggested a connection between bristle patterning and embryonic neuroblast development ). Other small Hsps such as Hsp26 and Hsp27 have been found through in situ hybridization to be expressed in embryonic central as well as peripheral nervous systems (BDGP database). Given the fact that RNAi knockdown of Hsp22 affected the mean bristle PC1 in this study, there may be co-localization of Hsp22 in the embryonic CNS tissues, and Hsp22 might function to buffer developmental perturbation in both tissues, but this remains to be tested.
Similar to Hsp22, knockdown of Hsp67Bc affects bristle trait asymmetry but not wing trait asymmetry. Hsp22 and Hsp67Bc are located at the chromosomal position of 67B on chromosome arm 3R, together with other small Hsps, forming the 67B Hsp cluster. RNAi knockdown of either gene alone was sufficient to alter bristle phenotype. This suggests that these two genes have unique and indispensable roles in bristle development. At the peptide level, Hsp22 and Hsp67Bc share no obvious similarity in two regions that are suggested to affect chaperone activity . At the transcriptional level, the tissue enrichment patterns of Hsp67Bc only partially overlap with those of Hsp22 (FLYATLAS http://www.flyatlas.org). Based on the dissimilarity in protein sequence, and non-overlapping spatial expression patterns, it is likely that Hsp22 and Hsp67Bc interact with different sets of protein clients, and contribute to bristle development through different mechanisms.
Hsp67Ba showed significant effects on both FA and among-individual variation of wing shape, suggesting that some developmental buffering mechanisms affect both within- and among-individual phenotypic variation. There was also a significant reduction of trait means in wing shape. However, there are at least 20 potential off-targets of this particular RNAi construct (VDRC construct ID 11237). It is therefore unclear whether the wing phenotypes we observed were due to successful suppression of Hsp67Ba, or collateral knockdown effects of these off-target genes. Future experiments should be directed to using an Hsp67Ba specific RNAi construct to assay for its causal knockdown phenotype. In fact, a new RNAi strain library, the KK library is available now at the Vienna Drosophila RNAi Center, in which each construct is targeted to the same position in the genome to reduce the number of off-target genes and positional effect. However coverage of the KK library is currently less than the RNAi strain library (GD library).
The RNAi lines 49795 and 49796 both target Hsp67Bb using the identical RNAi construct (VDRC construct ID 17696), and yet the phenotypic outcomes between these crosses were different. Bristle among-individual variation was significantly affected only in progeny derived from GAL4/49796. RTPCR results, however, indicated that Hsp67Bb was successfully suppressed in both crosses (Figure 4). One possible explanation for this inconsistent result is that there is a position effect of the construct, and the expression of Hsp67Bb may be reduced at different developmental stages in the RNAi lines. It is unclear at which pre-adult developmental stage(s) the mRNA knockdown exerts its effects on adult morphology. Future experiments could utilize stage-specific GAL4 driver lines to clarify the timing of the gene action.
The manner in which sHsps tested in the present study buffer developmental variation is unknown. Rutherford et al.  suggested a hypothetical mechanism of Hsp90 developmental buffering, based on the idea of thresholds for the expression of phenotypes in response to continuously varying strengths of signaling through Hsp90 targeted pathways. In their hypothesis, when Hsp90 levels are decreased, signal transduction clients begin to lose activity, and the strength of target pathways becomes severely reduced. In specific genetic interactions between Hsp90 and signaling pathways, reduction of the signaling to the threshold for the expression of mutant phenotype reveals cryptic variation. Whether sHsps affecting developmental buffering in this study have similar interactions with signaling pathways is unknown. The expression of sHsps is regulated by the steroid molting hormone ecdysone and other enhancer elements , suggesting a possibility of their interaction with a number of signaling pathways. Recently, Specchia et al.  suggested a novel hypothesis that Hsp90 prevents phenotypic abnormalities by suppressing the mutagenic activity of transposons. They found that functional alteration of Hsp90 resulted in transposon activation and the induction of morphological mutants which indicated that Hsp90 mutation or inhibition can generate new genetic variation by transposon-mediated mutagenesis . It is not known if sHsps affect transposon activity.
Our data provide information on a long standing argument about whether there is a single or multiple developmental buffering mechanisms controlling developmental stability and canalization . As suggested by some previous work on Hsp90 [6, 50] and this study, multiple developmental buffering mechanisms may operate with trait specific effects. We found that overall the correlation between within- and among-individual phenotypic variations was very weak for most traits. This suggests that developmental stability and canalization are not always mediated by the same molecular machinery. Given that four sHsps (Hsp22, Hsp67Ba, Hsp67Bb, and Hsp67Bc) showed gene-specific functions in terms of developmental buffering, phenotypic stability of an organism is probably maintained by multiple buffering mechanisms, activated by different environmental and/or genetic stresses for different traits.
We identified four small Hsp genes (Hsp22, Hsp67Bc, Hsp67Ba and Hsp67Bb) that may influence developmental stability and/or canalization, possibly through multiple buffering mechanisms. Our findings provide new insights into the endogenous roles of these heat shock genes, and contribute to the understanding of the molecular mechanisms controlling variability in morphological traits. They also highlight that multiple genes can influence patterns of phenotypic variability.
We wish to thank Sarah De Garis for assistance with morphological measurements. We also thank for Ryu Ueda for valuable comments on the experiments. This work was financially supported by Research Fellowships of the Japan Society for the Promotion of Science for Young scientist, and by Special Coordination Funds for Promoting Sciences and Technology of MEXT (The Ministry of Education, Sport, Culture, Science and Technology) of Japan to KHT. LR, AAH and SFL were supported by funding from the Australian Research Council through their Discovery and Fellowship schemes.
- Waddington CH: Canalization of development and the inheritance of acquired characters. Nature. 1942, 150: 563-565. 10.1038/150563a0.View ArticleGoogle Scholar
- Morimoto RI: Regulation of the heat shock transcriptional response: cross talk between a family of heat shock factors, molecular chaperones, and negative regulators. Genes Dev. 1998, 12: 3788-3796. 10.1101/gad.12.24.3788.View ArticlePubMedGoogle Scholar
- Rutherford S, Lindquist S: Hsp90 as a capacitor for morphological evolution. Nature. 1998, 396: 336-342. 10.1038/24550.View ArticlePubMedGoogle Scholar
- Queitsch C, Sangster TA, Lindquist S: Hsp90 as a capacitor of phenotypic variation. Nature. 2002, 417: 618-624. 10.1038/nature749.View ArticlePubMedGoogle Scholar
- Yeyati PL, Bancewicz RM, Maule J, van Heyningen V: Hsp90 selectively modulates phenotype in vertebrate development. PLoS Genet. 2007, 3: e43-10.1371/journal.pgen.0030043.PubMed CentralView ArticlePubMedGoogle Scholar
- Milton CC, Huynh B, Batterham P, Rutherford S, Hoffmann AA: Quantitative trait symmetry independent of Hsp90 buffering: Distinct modes of genetic canalization and developmental stability. Proc Natl Acad Sci USA. 2003, 100: 13396-13401. 10.1073/pnas.1835613100.PubMed CentralView ArticlePubMedGoogle Scholar
- Debat V, Milton CC, Rutherford S, Klingenberg CP, Hoffmann AA: HSP90 and the quantitative variation of wing shape in Drosophila melanogaster. Evolution. 2006, 60: 2529-2538.View ArticlePubMedGoogle Scholar
- Dahlgaard J, Loeschcke V, Michalak P, Justesen J: Induced thermotolerance and associated expression of the heat-shock protein Hsp70 in adult Drosophila melanogaster. Funct Ecol. 1998, 12: 786-793. 10.1046/j.1365-2435.1998.00246.x.View ArticleGoogle Scholar
- Sorensen JG, Michalak P, Justesen J, Loeschcke V: Expression of the heat-shock protein HSP70 in Drosophila buzzatii lines selected for thermal resistance. Hereditas. 1999, 131: 155-164. 10.1111/j.1601-5223.1999.00155.x.View ArticlePubMedGoogle Scholar
- Sorensen JG, Loeschcke V: Larval crowding in Drosophila melanogaster induces Hsp70 expression, and leads to increased adult longevity and adult thermal stress. J Insect Physiol. 2001, 47: 1301-1307. 10.1016/S0022-1910(01)00119-6.View ArticlePubMedGoogle Scholar
- Bettencourt BR, Kim I, Hoffmann AA, Feder ME: Response to natural and laboratory selection at the Drosophila HSP70 genes. Evolution. 2002, 56: 1796-1801.View ArticlePubMedGoogle Scholar
- Pedersen KS, Kristensen TN, Loeschcke V: Effects of inbreeding and rate of inbreeding in Drosophila melanogaster - Hsp70 expression and fitness. J Evol Biol. 2005, 18: 756-762. 10.1111/j.1420-9101.2005.00884.x.View ArticlePubMedGoogle Scholar
- Feder ME, Hofmann GE: Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999, 61: 243-282. 10.1146/annurev.physiol.61.1.243.View ArticlePubMedGoogle Scholar
- Williams KD, Helin AB, Posluszny J, Roberts SP, Feder ME: Effect of heat shock, pretreatment and hsp70 copy number on wing development in Drosophila melanogaster. Mol Ecol. 2003, 12: 1165-1177. 10.1046/j.1365-294X.2003.01771.x.View ArticlePubMedGoogle Scholar
- McColl G, Hoffmann AA, McKechnie SW: Response of two heat shock genes to selection for knockdown heat resistance in Drosophila melanogaster. Genetics. 1996, 143 (4): 1615-1627.PubMed CentralPubMedGoogle Scholar
- Bettencourt BR, Hogan CC, Nimali M, Drohan BW: Inducible and constitutive heat shock gene expression responds to modification of Hsp70 copy number in Drosophila melanogaster but does not compensate for loss of thermotolerance in Hsp70 null flies. BMC Biol. 2008, 6: 5-10.1186/1741-7007-6-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Calderwood SK: Evolving connections between molecular chaperones and neuronal function. Int J Hyperthermia. 2005, 21 (5): 375-378. 10.1080/02656730500165480.View ArticlePubMedGoogle Scholar
- Morrow G, Tanguay RM: Heat shock proteins and aging in Drosophila melanogaster. Seminars in Cell & Developmental Biology. 2003, 14: 291-299.View ArticleGoogle Scholar
- Bronk P, Wenniger JJ, Dawson-Scully K, Guo X, Hong S, Atwood HL, Zinsmaier KE: Drosophila Hsc70-4 is critical for neurotransmitter exocytosis in vivo. Neuron. 2001, 30 (2): 475-488. 10.1016/S0896-6273(01)00292-6.View ArticlePubMedGoogle Scholar
- Tower J: Hsps and aging. Trends in Endocrinology and Metabolism. 2009, 20: 216-222. 10.1016/j.tem.2008.12.005.View ArticlePubMedGoogle Scholar
- Mason PJ, Hall LMC, Gausz J: The expression of heat shock genes during normal development in Drosophila melanogaster (heat shock/abundant transcripts/developmental regulation). Mol Genet Genomics. 1984, 194: 73-78. 10.1007/BF00383500.View ArticleGoogle Scholar
- Pauli D, Arrigo AP, Vazquez J, Tonka CH, Tissieres A: Expression of the small heat shock genes during Drosophila development: comparison of the accumulation of hsp23 and hsp27 mRNAs and polypeptides. Genome. 1989, 31 (2): 671-676.View ArticlePubMedGoogle Scholar
- Pauli D, Tonka CH, Ayme-Southgate A: An unusual split Drosophila heat shock gene expressed during embryogenesis, pupation and in testis. J Mol Biol. 1988, 200 (1): 47-53. 10.1016/0022-2836(88)90332-4.View ArticlePubMedGoogle Scholar
- Debat V, David P: Mapping phenotypes: canalization, plasticity and developmental stability. Trends Ecol Evol. 2001, 16: 555-561. 10.1016/S0169-5347(01)02266-2.View ArticleGoogle Scholar
- Dietzl G, et al: A genome-wide transgenic RNAi library for conditional gene inactivation in Drosophila. Nature. 2007, 448: 151-157. 10.1038/nature05954.View ArticlePubMedGoogle Scholar
- Kellermann VM, Hoffmann AA, Sgro CM: Hsp90 inhibition and the expression of phenotypic variability in the rainforest species Drosophila birchii. Biol J Linn Soc. 2007, 92: 457-465. 10.1111/j.1095-8312.2007.00875.x.View ArticleGoogle Scholar
- Rohlf FJ, Slice D: Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990, 39: 40-59. 10.2307/2992207.View ArticleGoogle Scholar
- Bookstein FL: Morphometric tools for landmark data: geometry and biology. 1991, Cambridge: Cambridge University PressGoogle Scholar
- Rohlf FJ, Marcus LF: A revolution in morphometrics. Trends Ecol Evol. 1993, 8-Google Scholar
- Richtsmeier JT, Cheverud JM, Lele S: Advances in anthropological morphometrics. Annual Review of Anthropology. 1992, 21: 283-305. 10.1146/annurev.an.21.100192.001435.View ArticleGoogle Scholar
- Bookstein FL: Principal warps: thin-plate splines and the decomposition of deformations. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1989, 11: 567-585. 10.1109/34.24792.View ArticleGoogle Scholar
- Leung B, Forbes MR, Houle D: Fluctuating asymmetry as a bioindicator of stress: comparing efficacy of analyses involving multiple traits. Am Nat. 2000, 155: 101-115. 10.1086/303298.View ArticlePubMedGoogle Scholar
- Klingenberg CP, Monteiro LR: Distances and directions in multidimensional shape spaces: implications for morphometric applications. Syst Biol. 2005, 54: 678-688. 10.1080/10635150590947258.View ArticlePubMedGoogle Scholar
- Flury B, Riedwyl H: Standard distance in univariate and multivariate analysis. Am Stat. 1986, 40: 249-251. 10.2307/2684560.Google Scholar
- Flury B: A first course in multivariate statistics. 1997, New York: SpringerView ArticleGoogle Scholar
- Mardia KV, Kent JT, Bibby JM: Multivariate analysis. 1979, London: Academic PressGoogle Scholar
- Arnqvist G, Martensson T: Measurement error in geometric morphometrics: empirical strategies to assess and reduce its impact on measures of shape. Acta Zool Acad Sci Hung. 1998, 44: 73-96.Google Scholar
- Klingenberg CP, McIntyre GS: Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with procrustes methods. Evolution. 1998, 52: 1363-1375. 10.2307/2411306.View ArticleGoogle Scholar
- Klingenberg CP, Barluenga M, Meyer A: Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry. Evoluion. 2002, 56: 1909-1920.View ArticleGoogle Scholar
- Bogwitz MR, Chung H, Magoc L, Rigby S, Wong W, O'Keefe M, McKenzie JA, Batterham P, Daborn PJ: Cyp12a4 confers lufenuron resistance in a natural population of Drosophila melanogaster. Proc Natl Acad Sci USA. 2005, 102: 12807-12812. 10.1073/pnas.0503709102.PubMed CentralView ArticlePubMedGoogle Scholar
- Ni J, Markstein M, Binari R, Pfeiffer B, Liu L, Villalta C, Booker M, Perkins L, Perrimon N: Vector and parameters for targeted transgenic RNA interference in Drosophila melanogaster. Nat Methods. 2008, 5: 49-51. 10.1038/nmeth1146.PubMed CentralView ArticlePubMedGoogle Scholar
- Pena-Rangel MT, Rodriguez I, Riesgo-Escovar JR: A misexpression study examining dorsal thorax formation in Drosophila melanogaster. Genetics. 2002, 160 (3): 1035-1050.PubMed CentralPubMedGoogle Scholar
- Simpson P: The stars and stripes of animal bodies: evolution of regulatory elements mediating pigment and bristle patterns in Drosophila. Trends Genet. 2007, 23 (7): 350-358. 10.1016/j.tig.2007.04.006.View ArticlePubMedGoogle Scholar
- Calleja M, Herranz H, Estella C, Casal J, Lawrence P, Simpson P, Morata G: Generation of medial and lateral dorsal body domains by the pannier gene of Drosophila. Development. 2000, 127 (18): 3971-3980.PubMedGoogle Scholar
- Held LI: Bristle patterning in Drosophila. Bioessays. 1991, 13 (12): 633-640. 10.1002/bies.950131203.View ArticlePubMedGoogle Scholar
- Nakamoto H, Vigh L: The small heat shock proteins and their clients. Cell Mol Life Sci. 2007, 64 (3): 294-306. 10.1007/s00018-006-6321-2.View ArticlePubMedGoogle Scholar
- Rutherford S, Knapp JR, Csermely P: Hsp90 and developmental networks. Molecular aspect of the stress response: Chaperones, Membranes and Networks. 2007, New York: Springer New YorkGoogle Scholar
- Specchia V, Piacentini L, Tritto P, Fanti L, D'Alessandro R, Palumbo G, Pimpinelli S, Bozzetti MP: Hsp90 prevents phenotyic variation by suppressing the mutagenic activity of transposons. Nature. 2010, 463: 662-665. 10.1038/nature08739.View ArticlePubMedGoogle Scholar
- Meiklejohn CD, Hartl DL: A single mode of canalization. Trends Ecol Evol. 2002, 17: 468-473. 10.1016/S0169-5347(02)02596-X.View ArticleGoogle Scholar
- Milton CC, Batterham P, McKenzie JA, Hoffmann AA: Effect of E (sev) and Su (Raf) Hsp83 mutants and Trans-heterozygotes on bristle trait means and variation in Drosophila melanogaster. Genetics. 2005, 171: 119-130. 10.1534/genetics.104.038463.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.