BACKGROUND
in drosophila and other arthropods, pulses of the steroid hormone 20-hydroxyecdysone  are responsible for the temporal coordination of larval molts and metamorphosis.
physiological responses during these events can be diverse; for example during metamorphosis obsolete larval tissues are destroyed and adult structures arise from imaginal disc cells.
remarkably, these actions are carried out in a coordinated, tissue specific manner.
at the site of target tissues, 20e binds to its cognate nuclear receptor triggering a cascade of gene activation.
primary 20e-inducible genes, which are directly induced by the steroid-receptor complex, are the earliest genes in the cascade to be transcribed and are insensitive to protein synthesis inhibitors.
in contrast, secondary 20e-inducible genes are expressed later and are dependent on the synthesis of primary-response genes.
early studies of the 20e cascade looking at the puffing patterns of polytene chromosomes of late 3rd instar larvae in response to 20e, predicted that primary-response genes would code for proteins that are responsible both for the induction of secondary response genes as well as for the inhibition of their own transcription  <cit> .
three of the most well characterized 20e-primary response genes, eip74ef, eip75b, and br fit this description perfectly.
furthermore, these three genes reside at chromosome cytolocations 74ef, 75b, and 2b <dig> respectively, which, along with approximately three other loci, exhibit rapid and dramatic puffing after exposure to 20e either naturally or artificially  <cit> .
likewise, the earliest characterized secondary response genes are found in a region  that forms a distinct yet delayed puff upon 20e exposure  <cit> .
although the precise function of these genes has yet to be determined, based on their genomic sequence they are thought to encode effector proteins which is consistent with early predictions of secondary-response gene function  <cit> .
since the original characterization of these 20e-response genes, however, many examples of primary and secondary response genes with diverse functions have emerged, and many of these 20e-inducible genes do not appear to be associated with any identifiable 20e-induced puffs  <cit>  underling the importance identifying individual components of the cascade for a clearer picture of 20e action.
the transcription factor complex coordinating the entire 20e-hierarchy, the 20e nuclear receptor heterodimer is composed of the vertebrate retinoid x receptor homolog, ultraspiracle  and the ecdysone receptor .
ecr isoform distribution differs between tissues and typically those tissues expressing different isoforms show different responses to ecdysone at metamorphosis  <cit> .
for example, immunolocalization experiments examining ecr isoform expression pattern at the onset of metamorphosis show that larval tissues destined to die during metamorphosis, such as the salivary gland, stain strongly for isoform b <dig> and weakly for isoform a whereas proliferating tissues, such as imaginal discs, show a reciprocal pattern staining strongly for isoform a and weakly for isoform b <dig>  <cit> .
the most recent annotation of the drosophila genome documents five ecr isoforms that differ in sequence at the amino terminus but share common dna- and ligand- binding domains.
the ligand-binding pocket of the 20e receptor shows remarkable flexibility enabling its activation by a variety of steroidal and non-steroidal 20e analogs/agonists  <cit> .
for example, the ligand-binding domain of the drosophila ecr homolog in heliothis virescens adopts different structures in the presence of different 20e analogs  <cit> .
of the insect and plant derived analogs characterized, ponsterone a  is the most potent agonist of the ecr receptor  <cit>  with an affinity approximately eight times greater than 20e  <cit> .
as such, poa, which differs from 20e by the abscence of a single hydroxyl group, has been used to work out numerous kinetic and physical parameters of the ecr  <cit> .
interestingly, the natural role of this and other phytoecdysteroids is still under debate, although most evidence disfavors a hormonal role in plants  <cit> .
instead, it is predicted that phytoecdysteroids induce precocious molting and subsequently death in insects providing an effective defense against insect feeding  <cit> .
given the increased binding affinity of poa over 20e as well as its potential use in insect control, it is of importance to understand how this phytoecdysteroid affects transcription on a global scale.
the drosophila kc cell line is one of the most well documented 20e-responsive cell lines currently available and has been used in numerous studies examining the effects of 20e.
in this study we use cdna-based microarrays representing approximately 80% of the drosophila genome to identify 20e-responsive genes in kc <dig> cells.
by exploiting the sensitivity of secondary-response gene transcription to protein synthesis, we are able to determine which of these genes are primary and which are secondary 20e-inducible genes.
furthermore, the transcriptional response to physiological levels of 20e is compared to that elicited under two other conditions: a 20-fold higher concentration of 20e or its plant derived structural analog, poa.
these analyses lead to the identification of  <dig> genes that reacted similarly to all three treatments.
in addition to examining the transcriptional response to 20e in kc cells, we were also interested in examining the response in a natural target tissue of 20e signaling, salivary glands.
salivary glands have been critical in the elucidation of the 20e signaling hierarchy through the examination of puff patterns on polytene chromosomes.
by focusing on the response of the salivary gland in the isolation from other larval tissues we were able to identify many 20e-responsive genes that were not detected when whole larvae were examined  <cit> .
taken together, this work provides a detailed picture of the genes involved in the ecdysone hierarchy of gene transcription in kc <dig> cells and salivary glands.
RESULTS
identification of 20-hydroxyecdysone -responsive genes in kc <dig> cells
the 20e-hierarchy of gene transcription serves as a good model for examining hormonal control of development.
to gain insight into those genes that are part of the hierarchy, we assayed for 20e-induced changes in gene transcription across the genome.
in the first part of this study, the transcriptional response to  <dig>  μm 20e was examined in the drosophila cell line, kc <dig>  because  <dig>  μm 20e causes the polytene chromosome puffs at the 74ef and 75b early 20e-inducible loci to reach their maximum size after  <dig> hours of organ culture  <cit>  and in situ hybridization using labeled rna produced after 20e treatment shows puff specific rnas also increase along with puff size  <cit>  the transcriptional response was examined after a  <dig> hour treatment with  <dig>  μm 20e.
the transcriptional profile of cells exposed to  <dig>  μm 20e for  <dig> hours was also examined to identify any genes that may be induced early on but turned off after  <dig> hours.
to allow identification of 20e-responsive genes in a controlled extracellular environment, in the absence of endogenous hormones, these experiments were conducted in vitro.
kc <dig> cells were used because the kc cell line is responsive to 20e both morphologically  <cit>  and transcriptionally  <cit>  and is commonly used for the study of 20e.
to examine the transcriptional response to 20e on genome-wide scale, rna was isolated from both 20e-treated cells and untreated control cells, reverse-transcribed and labelled with cy <dig> or cy <dig> dye-coupled nucleotides respectively, and co-hybridized to cdna microarrays.
following data acquisition and analysis , transcripts with at least a  <dig>  fold change in abundance in response to 20e were identified.
after a  <dig> hour treatment with 20e,  <dig> genes were induced  .
two-thirds of these genes  including the known 20e-inducible genes eip <dig>  br, eip28/ <dig>  and eip <dig>  remained significantly up-regulated after  <dig> hours in the presence of 20e while  <dig> additional genes were induced and  <dig> genes were repressed  .
among the 20e-induced genes indentified here in kc cells, the most strongly induced after  <dig> hours, eip28/ <dig>  was originally identified owing to its quick reaction to 20e in kc cells  <cit>  and has since been used to study aspects of the 20e response such as its tissue specificity  <cit> .
taken together, the response observed here after a  <dig> and  <dig> hour 20e-treatment are consistent with two prior observations on the 20e-induced puffing patterns of larval salivary glands: first, that primary 20e-response genes are likely to still be detectable after a four-hour 20e-treatment and second, that the 20e-response of kc <dig> is more robust after  <dig> hours of 20e exposure  <cit>  .
of the 20e-inducible genes identified here, 30% were also identified in a recent study by gauhar and colleagues that examined the transcriptional response of kc cells following either a  <dig>   <dig>  or  <dig> hr 20e treatment  <cit> .
however, if we apply more stringent fold-change criteria to the list of 20e indentified here the degree of overlap increases such that at a 2-fold cut-off, there is 50% overlap, and there is 100% overlap between ours and the gauhar et al.
lists if we select a 3-fold change in expression.
remaining differences in the identification of 20e responsive genes indentified is likely due in part to the use of different culturing medium as we have previously found that media composition affects the transcriptional response of kc cells to 20e .
the use of a different microarray platform in the measurement of the transcriptional profile could also be a contributing factor.
identification of primary 20e-response genes
having identified genes involved in the 20e-response, we were interested in determining where these genes fit into the 20e-hierarchy of gene activation.
in particular, we wanted to identify those genes comprising the primary response to the hormone and are thus likely to be directly induced by 20e.
maximum induction of primary-response genes by 20e is expected to occur after  <dig> hours of exposure to the hormone, however, some secondary-response gene transcription may have already begun at this point  <cit> .
one distinguishing characteristic of primary and secondary-response genes is their dependence on protein synthesis for transcription; only primary response genes are transcribed in the absence of protein synthesis  <cit> .
thus, to identify primary response genes, cells were treated with 20e in conjunction with an inhibitor of protein synthesis, cycloheximide .
rna isolated from treated cells was labelled and co-hybridized to cdna arrays along with a differentially labelled untreated control.
from this analysis we identified  <dig> genes that were induced and  <dig> genes that were repressed by the 20e-cycloheximide combination .
to determine which of these genes were induced/repressed in response to 20e, and were not responding to exposure to cycloheximide alone, additional microarrays were run using rna from kc <dig> cells treated only with cycloheximide.
from this experiment,  <dig> genes were indentified whose transcription is altered due solely to the inhibition of protein synthesis .
a two-class sam analysis identified  <dig> genes whose transcripts levels in the presence of cycloheximide differs significantly when 20e is present including the well characterized 20e-primary-response genes, eip75b and br .
although eip74ef was not among the genes on this list, we confirmed its induction by 20e by qrt-pcr  suggesting that the cdna probe on the batch of arrays used for this study was likely of poor quality and unable to detect transcripts of this gene.
of the  <dig> primary-response genes identified,  <dig> were not induced by 20e alone.
the induction of these  <dig> genes by 20e is likely masked in the absence of cycloheximide due to a negatively regulated auto-feedback loop where primary-response gene transcription is repressed by primary-response gene protein products  <cit> .
interestingly, the only gene transcript repressed as part of the primary response to 20e, cg <dig>  also evades detection unless protein synthesis is blocked .
the majority of 20e-responsive genes identified after a  <dig> hour treatment with  <dig>  μm 20e behaved as secondary-response genes as their response to 20e was blocked in the presence of cycloheximide;  <dig> genes were induced and  <dig> were repressed as part of the secondary response .
a two- class sam analysis of triplicate data identified primary-response genes that responded to 20e in the presence of cycloheximide.
the top  <dig> up-regulated primary-response genes and single down-regulated primary-response gene are listed .
secondary-response genes regulated by 20e only in the absence of cycloheximide were also identified.
the top  <dig> up- and down-regulated secondary-response genes are listed .
the fold change of each gene following a  <dig> hour treatment with  <dig>  μm 20e  and, for primary-response genes, the fold change following a  <dig> hour treatment with  <dig>  μm 20e and  <dig> μm cycloheximide  is given.
q-values determined by sam represent the lowest false discovery rate at which that gene is considered significant.
for a gene to be included as differentially expressed, it must have a q-value < 5% in that experimental condition.
q-values > 5% are not shown.
david  was used to identify enriched annotation categories among groups of genes that respond similarly to 20e.
some of the most enriched functional terms associated with each group  are given in the right most column.
a complete list of all genes in these groups and their functional annotation is in additional file  <dig>
in addition to differences in transcriptional response in the absence of protein synthesis, it is also expected that most primary-response genes encode regulators while secondary-response genes encode effectors  <cit> .
to determine if this is the case for the genes identified here as part of the primary and secondary response, we examined each group of genes for enrichment in functional annotation terms/keywords using the online resource david.
the results of this analysis are summarized in table  <dig>  several of the primary-response gene transcripts include those that code for the well known dna binding transcription factors eip75b and br as well as vrille which is consistent with the model that primary-response proteins are required for the transcriptional induction of the secondary-response genes.
"induction of programmed cell death by hormones" is among the most highly enriched terms associated with primary-response genes and comes as no surprise due to the well established involvement of 20e-regualted genes such as br and nc in cell death pathways during drosophila development  <cit> .
however, we should point out that that kc cells do not undergo apoptotic cell death in response to ecdysone treatment and therefore other factors and/or conditions must need to be present in order to carry out this particular process.
"catalytic activity" and "hydrolase activity" were both also among the most significantly enriched terms supporting the idea that primary-response genes are involved in a broad range of regulatory roles  <cit> .
for example, one of the genes belonging to the hydrolase class includes imaginal disk growth factor  <dig>  a protein that when bound to its receptor initiates signal transduction cascades important to imaginal disk development.
some primary- and secondary-response genes include genes that code for proteins involved in cell movement and organization and/or are associated with the cytoskeleton.
these include up-regulated transcripts for roadblock, a gene coding for a dynein-associated protein and down-regulated transcripts for spastin, and a gene coding for a microtubule severing protein.
transcript changes in genes of this sort are consistent with the cell movement and morphogenetic changes that occur during 20e-dependent developmental changes.
also of interest is that several of the secondary-response genes include genes that are involved in metabolic processes.
more specifically, there seems to be an increase in transcripts for genes involved in mitochondrial respiration and a decrease in transcripts for genes involved in glycolysis perhaps suggesting that in addition to the cellular organization and tissue changes that are induced by 20e, there is also a shift in how energy is being produced.
transcriptional response of kc <dig> cells to an increased concentration of 20e and to its structural analog, ponasterone a
to examine the specificity of the 20e-response of kc <dig> cells to both hormone concentration and ligand structure, two other conditions were tested: a 20-fold higher concentration of 20e,  <dig> μm, and the replacement of 20e with its structural analog, poa.
20e concentration is known to affect at least three aspects of the 20e-response: the size of both early and late puffs  <cit> , the rate of early puff regression  <cit> , and the transcription of some early gene isoforms  <cit> .
poa, often used as substitute for 20e, is a more potent activator of the 20e-response.
at least four different measures indicate that poa activity is approximately eight fold higher than 20e activity: affinity for the ecdysone receptor  <cit> , morphological changes associated with the 20e-response  <cit> , the level of induction of known 20e-inducible proteins  <cit> , and the level of transcription from an ecdysone response element   <cit> .
thus, to minimize effects caused by the difference in receptor affinity and make the comparison between analogs as similar as possible, our working concentration of poa was eight fold less than that of 20e .
genes that responded to either of the two treatments were identified by microarray analysis as described above.
compared to the number of genes induced by a  <dig> hour treatment with  <dig>  μm 20e, both  <dig> μm 20e and  <dig>  μm poa stimulated the transcription of a greater number of genes:  <dig> and  <dig> respectively versus  <dig> induced by  <dig>  μm 20e .
furthermore, poa down-regulated substantially more genes than either 20e treatment.
one-hundred and fifty genes were repressed after a  <dig> hour treatment with poa - almost double the number repressed by  <dig>  μm 20e  and nearly six fold more than the number repressed by  <dig> μm 20e  .
it is possible that the differences in the response to poa may be due to the activation of additional nuclear receptors; it is only assumed to be the 20e receptor heterodimer because it has been shown to be the case in vitro  <cit> .
overall,  <dig> genes were identified that responded to at least one of the three hormone treatments examined .
to determine which of these  <dig> hormone-responsive genes were sensitive to all three conditions tested and which are differentially regulated by either  <dig> μm 20e or poa relative to  <dig>  μm 20e, a multi-class sam analysis was performed where the transcriptional response to  <dig> μm 20e and poa was compared to the transcriptional response  <dig>  μm 20e for each gene.
this analysis revealed two things about the transcriptional response to each of the hormone conditions tested.
first, of the  <dig> genes that responded to at least one of the three treatments,  <dig> showed no significant difference in their transcriptional response to all three treatments and eight of these are primary 20e-inducible genes .
second, poa produced a more distinct transcriptional profile than did an elevated concentration of 20e with respect to the  <dig>  μm 20e treatment , suggesting that poa may not be acting solely through the ecdysone hierarchy.
we confirmed the transcriptional response of three key 20e-inducible genes, the ecdysone-receptor complex, and a novel 20e-incudible gene by qrt-pcr.
in all cases, the transcriptional response for these genes measured by qrt-pcr, although differing in magnitude when compared to the microarray analysis , had the same direction of change .
cg5346
pect
prl-1
cg4825
eater
cyp4e2
cg14523
granny-smith
a multi-class sam analysis revealed genes with similar transcriptional responses to the hormones and hormone concentrations studied.
genes identified as primary 20e-response genes  are indicated in bold.
as with table  <dig>  the fold change and q-value of each gene following the treatment indicated is given.
a missing fold change value means that gene was either not represented by at least 2/ <dig> spots and/or did not exhibit a coefficient of variation of at least  <dig> for that treatment across replicate arrays.
given the observed similarities and differences in the ecdysteroid-response to both an increased 20e-concentration and the use of a structural analog, we wanted to identify functions associated with genes that responded similarly to all treatments tested.
using david we found two categories enriched in similarly responding genes that were also enriched among primary-response genes: "catalytic activity" and "hydrolase activity" .
conversely, genes that reacted significantly different to either  <dig> μm 20e, poa or  <dig>  μm 20e were enriched in many unique categories; the most strongly enriched functional terms are associated with genes down-regulated by poa relative to  <dig>  μm 20e .
furthermore, other categories related to development are also strongly enriched in the list of genes that are down-regulated by poa relative to  <dig>  μm 20e such as "anatomical structure development", "induction of programmed cell death", and "instar larval or pupal development" .
the apparent lack of activation of developmental genes by poa, as indicated by this analysis, is suggestive of an alternate mode of action of poa induced death in non-adapted insects more so than the simple induction of premature molting via the canonical ecdysone hierarchy.
in fact, the results potentially suggest that poa might interfere with normal 20e induced transcriptional changes and prevent the induction of some of the genes that are required for 20e dependent developmental processes.
three conditions
both poa
expression differences are relative to the  <dig>  μm 20e dataset .
david  was used to identify enriched annotation categories among groups of genes that responded similarly to one or more of the hormone treatments examined in this study .
the modified fisher exact p-value represents the probability that the number of genes associated with an annotation term in the sample list is random compared to the occurrence of the term in the all of the genes on the array used for the analysis ).
biological process;  molecular function;  cellular component;  swiss prot keywords;  interpro name;  uniprot sequence feature.
transcriptional response of salivary glands to 20e
we wanted to extend our analysis to a 20e-responsive larval tissue to see how it compares to the 20e transcriptional response of kc <dig> cells.
salivary glands were chosen in this respect due to the relative ease of dissection and because the foundations of our understanding of the 20e response comes from the study of salivary glands.
cultured glands from  <dig> late third instar larvae staged by the blue gut method  were treated with  <dig>  μm 20e for four hours.
we chose to treat salivary glands prior to the natural 20e pulse rather than assay glands after the natural pulse to rule out potential confounding effects caused by other endogenous hormones and/or the presence of protein products of secondary-response genes within the 20e hierarchy.
rna was isolated from 20e treated glands and untreated control glands, amplified, and analyzed by microarray experiments as described in the methods section.
we identified  <dig> up-regulated genes including several known 20e-inducible genes such as, eip75b, eip74ef, ecr, nc, and impl <dig>  and  <dig> down-regulated genes .
in a study by beckstead and colleagues looking at 20e responsive genes in cultured larval organs  <dig> genes were identified as 20e responsive  <cit> .
their treatments, however, utilized twice the concentration of 20e and lasted for an additional two hours.
nevertheless, we were still able to identify  <dig> genes in common, the majority of which showed similar changes in expression level .
by focusing only on salivary glands in this study we were able to identify  <dig> salivary gland specific 20e-inducible genes that were not identified in total organ preparations  <cit>  .
a david analysis reveals that these salivary gland specific genes are enriched for protein binding .
genes induced by 20e in salivary glands but not in kc <dig> cells were enriched in go biological processes such as "salivary gland histolysis", "programmed cell death" and various "developmental" processes as might be expected.
"cellularization" and "transporter activity" were also among the highest scoring hits.
conversely, genes that were up regulated only in kc cells showed enrichment for "metabolic process" and "oxidoreductase activity"- categories that were also enriched among secondary response genes in kc cells suggesting that the secondary response genes differ between kc cells and salivary glands .
comparison of 20e responsive genes from both systems allowed the identification of  <dig> genes whose response is conserved across both systems suggesting that they might occupy key positions within the hierarchy .
among the most highly induced of these genes are two genes that have not been traditionally associated with the 20e response in drosophila, cg <dig> and rrp <dig> .
cg <dig> is currently un-annotated and rrp <dig> is a component of the exosome that has been shown to relocate to developmental loci during periods of active transcription and is involved in mrna processing  <cit> .
rrp <dig> showed the greatest correlation in expression profile to eip75b highlighting the importance of gene transcription in the ecdysone response.
david  was used to identify enriched annotation categories among groups of genes that respond similarly to 20e.
the modified fisher exact p-value represents the probability that the number of genes associated with an annotation term in the sample list is random compared to the occurrence of the term in the all of the genes on the array used for the analysis ).
biological process;  molecular function;  cellular component;  swiss prot keywords;  interpro name;  uniprot sequence feature.
tissue specific differences in ecdysone receptor  isoform induction by 20e
some of the differences in the transcriptional response of kc cells and salivary glands may be explained by differential chromatin states of inactive genes in the two systems.
for example, salivary gland secretion protein  genes known be under the control of the hierarchy were expressed in glands but not in cells most likely because of the difference in the two cell types.
differential gene regulation by 20e may also result from the expression of different ecr isoforms or combinations thereof  <cit> .
typically, tissues with dissimilar developmental fates express different ecr isoforms  <cit> .
our microarray results agree with previous findings that kc cells express hemocyte and plasmatocyte marker genes pxn, ush, and hml  <cit>   and that 20e treatment of kc cells causes the induction of eip28/ <dig> and eip <dig> , both of which are 20e-regulated in the hematopoietic lymph gland  <cit>  providing further support for an embryonic hemocyte origin of kc cells  <cit> .
in contrast to larval salivary glands, embryonic hemocytes persist in the adult fly  <cit>  suggesting that kc cells are likely more similar to imaginal disc or adult tissues than larval tissues and placing them in a different metamorphic class than the salivary gland  <cit> .
because ecr is a 20e-inducible gene, we first examined the overall level of induction of ecr by qrt-pcr using a primer designed against the 3' region of the gene common to all isoforms .
interestingly, unlike many of the 20e-inducible genes we identified by microarray analysis, ecr was induced to a similar extent in poa-treated kc cells, 20e-treated kc cells and 20e-treated salivary glands .
we next examined the transcriptional response of each ecr isoform to poa and 20e using primers that recognize unique regions of each isoform and found that there was no difference in isoform regulation in kc cells treated with either hormone .
comparison of isoform induction in 20e-treated kc cells to 20e-treated glands, however, revealed a tissue specific pattern in isoform regulation .
at least two ecr isoforms  were induced in kc cells by both poa and 20e but not induced at all in salivary glands treated with 20e .
it has previously been shown that the proteins encoded by the ecr-ra, ecr-rd and ecr-re isoforms show strong expression in both embryos and imaginal discs  <cit>  and ecr-a is thought to be responsible for adult differentiation  <cit> ecr-rc, on the other hand, although induced by both poa and 20e in kc cells, was more strongly induced by 20e in salivary glands .
although differences in ecr isoform regulation in kc cells and salivary glands by 20e may be linked to the variation in the observed transcriptional response to 20e, it is unlikely that the same explanation can be applied to account for the difference in response observed in kc cells treated with 20e and poa since qrt-pcr analysis reveals at best, subtle differences in ecr isoform regulation by 20e and poa .
it is tempting to speculate that the different ecr isoforms form 20e receptors with different target genes and thus contribute to differences in the 20e induced transcript profiles.
alternatively, the differences in transcript profiles may have very little to with which ecr isoforms are present and more to do with cell type differences of other proteins.
for example, differences in the concentrations of particular chromatin proteins could affect the accessibility of the various target genes differently in different cell types.
also, the relative amount of proteins that partner with ecr to make the 20e heterodimer nuclear receptor could affect the sites the 20e receptor binds to.
while the classical partner for ecr has been thought to be ultraspiracle , it is now known that biologically active 20e receptors can be made without usp  <cit> .
in mid-third instar larvae, the 20e dependent induction of the glue proteins requires ecr but not usp suggesting that ecr makes 20e receptors that consist of ecr and yet to be characterized partner protein  <cit> .
one approach to assess the relative contributions of each of the ecr isoforms would be to employ rnai knockdown strategies of the different ecr isoforms in both kc cells and larval tissues and examine what affect this has on the 20e induced transcript profiles.
CONCLUSIONS
we have identified  <dig> primary 20e-response genes that are induced by 20e in the absence of protein synthesis in drosophila kc <dig> cells, one of the most widely used cell lines employed to study the ecdysone response.
the primary 20e responsive genes were enriched for gene ontology  terms such as induction of programmed cell death by hormones, catalytic activity, hydrolase activity, and cytochrome p <dig>  we have also identified  <dig> secondary 20e-response genes .
go terms that are enriched in the secondary-response genes include an increase in transcripts for genes involved in mitochondrial respiration and a decrease in transcripts for genes involved in glycolysis suggesting that 20e induces shifts in cellular metabolism.
some primary- and secondary-response genes include genes that code for proteins involved in cell movement and organization and/or are cytoskeletal associated which is consistent with the cell movement and morphogenetic changes that occur during 20e-dependent developmental changes.
comparison of the genome-wide transcriptional response to 20e to its plant derived structural analog ponasterone a  revealed a large difference in the transcriptional targets of these molecules.
while these two compounds are structurally very similar, many more genes related to various aspects of development appear to be significantly induced by 20e than by poa.
more specifically, the most strongly enriched functional terms are associated with genes down-regulated  by poa relative to  <dig>  μm 20e.
the genes not induced by poa include go terms such as developmental process, cell differentiation, gamete generation, anatomical structure development, organelle organization and biogenesis, neuron development, and induction of programmed cell death.
the apparent lack of activation of developmental genes by poa suggests that poa induced death in non-adapted insects may be due to more than the just the simple induction of premature molting via the canonical ecdysone hierarchy.
we also compared the 20e response in kc cells to that of a natural 20e target tissue where the function of 20e has been well described, the salivary glands of wandering 3rd instar larvae, and found little overlap in 20e-responsive genes.
genes induced by 20e in salivary glands but not in kc <dig> cells were enriched in go biological processes such as "salivary gland histolysis", "programmed cell death" and various "developmental" processes as well as "cellularization" and "transporter activity".
to help identify a potential mechanism to explain the difference in the transcriptional responses in kc cells and salivary glands, we analyzed the 20e-induced transcription of the various ecr isoforms, a known 20e-inducible gene and a component of the nuclear receptor complex that binds 20e.
we did find differences in the induction of ecr isoforms ecr-ra, er-rc, and ecr-rd/e between kc cells and salivary glands in response to 20e.
this suggests that the relative amount of various ecr isoforms present in a cell is different in different cell types and possibly contributes to the transcriptional response a given tissue has to 20e.

 in drosophila and other arthropods, pulses of the steroid hormone 20-hydroxyecdysone are responsible for the temporal coordination of larval molts and metamorphosis. physiological responses during these events can be diverse; for example during metamorphosis obsolete larval tissues are destroyed and adult structures arise from imaginal disc cells. remarkably, these actions are carried out in a coordinated, tissue specific manner. at the site of target tissues, 20e binds to its cognate nuclear receptor triggering a cascade of gene activation. primary 20e-inducible genes, which are directly induced by the steroid-receptor complex, are the earliest genes in the cascade to be transcribed and are insensitive to protein synthesis inhibitors. in contrast, secondary 20e-inducible genes are expressed later and are dependent on the synthesis of primary-response genes. early studies of the 20e cascade looking at the puffing patterns of polytene chromosomes of late 3rd instar larvae in response to 20e, predicted that primary-response genes would code for proteins that are responsible both for the induction of secondary response genes as well as for the inhibition of their own transcription. three of the most well characterized 20e-primary response genes, eip74ef, eip75b, and br fit this description perfectly. furthermore, these three genes reside at chromosome cytolocations 74ef, 75b, and 2b respectively, which, along with approximately three other loci, exhibit rapid and dramatic puffing after exposure to 20e either naturally or artificially. likewise, the earliest characterized secondary response genes are found in a region that forms a distinct yet delayed puff upon 20e exposure. although the precise function of these genes has yet to be determined, based on their genomic sequence they are thought to encode effector proteins which is consistent with early predictions of secondary-response gene function. since the original characterization of these 20e-response genes, however, many examples of primary and secondary response genes with diverse functions have emerged, and many of these 20e-inducible genes do not appear to be associated with any identifiable 20e-induced puffs underling the importance identifying individual components of the cascade for a clearer picture of 20e action. the transcription factor complex coordinating the entire 20e-hierarchy, the 20e nuclear receptor heterodimer is composed of the vertebrate retinoid x receptor homolog, ultraspiracle and the ecdysone receptor. ecr isoform distribution differs between tissues and typically those tissues expressing different isoforms show different responses to ecdysone at metamorphosis. for example, immunolocalization experiments examining ecr isoform expression pattern at the onset of metamorphosis show that larval tissues destined to die during metamorphosis, such as the salivary gland, stain strongly for isoform b and weakly for isoform a whereas proliferating tissues, such as imaginal discs, show a reciprocal pattern staining strongly for isoform a and weakly for isoform b. the most recent annotation of the drosophila genome documents five ecr isoforms that differ in sequence at the amino terminus but share common dna- and ligand- binding domains. the ligand-binding pocket of the 20e receptor shows remarkable flexibility enabling its activation by a variety of steroidal and non-steroidal 20e analogs/agonists. for example, the ligand-binding domain of the drosophila ecr homolog in heliothis virescens adopts different structures in the presence of different 20e analogs. of the insect and plant derived analogs characterized, ponsterone a is the most potent agonist of the ecr receptor with an affinity approximately eight times greater than 20e. as such, poa, which differs from 20e by the abscence of a single hydroxyl group, has been used to work out numerous kinetic and physical parameters of the ecr. interestingly, the natural role of this and other phytoecdysteroids is still under debate, although most evidence disfavors a hormonal role in plants. instead, it is predicted that phytoecdysteroids induce precocious molting and subsequently death in insects providing an effective defense against insect feeding. given the increased binding affinity of poa over 20e as well as its potential use in insect control, it is of importance to understand how this phytoecdysteroid affects transcription on a global scale. the drosophila kc cell line is one of the most well documented 20e-responsive cell lines currently available and has been used in numerous studies examining the effects of 20e. in this study we use cdna-based microarrays representing approximately 80% of the drosophila genome to identify 20e-responsive genes in kc cells. by exploiting the sensitivity of secondary-response gene transcription to protein synthesis, we are able to determine which of these genes are primary and which are secondary 20e-inducible genes. furthermore, the transcriptional response to physiological levels of 20e is compared to that elicited under two other conditions: a 20-fold higher concentration of 20e or its plant derived structural analog, poa. these analyses lead to the identification of genes that reacted similarly to all three treatments. in addition to examining the transcriptional response to 20e in kc cells, we were also interested in examining the response in a natural target tissue of 20e signaling, salivary glands. salivary glands have been critical in the elucidation of the 20e signaling hierarchy through the examination of puff patterns on polytene chromosomes. by focusing on the response of the salivary gland in the isolation from other larval tissues we were able to identify many 20e-responsive genes that were not detected when whole larvae were examined. taken together, this work provides a detailed picture of the genes involved in the ecdysone hierarchy of gene transcription in kc cells and salivary glands.  identification of 20-hydroxyecdysone -responsive genes in kc cells the 20e-hierarchy of gene transcription serves as a good model for examining hormonal control of development. to gain insight into those genes that are part of the hierarchy, we assayed for 20e-induced changes in gene transcription across the genome. in the first part of this study, the transcriptional response to μm 20e was examined in the drosophila cell line, kc because μm 20e causes the polytene chromosome puffs at the 74ef and 75b early 20e-inducible loci to reach their maximum size after hours of organ culture and in situ hybridization using labeled rna produced after 20e treatment shows puff specific rnas also increase along with puff size the transcriptional response was examined after a hour treatment with μm 20e. the transcriptional profile of cells exposed to μm 20e for hours was also examined to identify any genes that may be induced early on but turned off after hours. to allow identification of 20e-responsive genes in a controlled extracellular environment, in the absence of endogenous hormones, these experiments were conducted in vitro. kc cells were used because the kc cell line is responsive to 20e both morphologically and transcriptionally and is commonly used for the study of 20e. to examine the transcriptional response to 20e on genome-wide scale, rna was isolated from both 20e-treated cells and untreated control cells, reverse-transcribed and labelled with cy or cy dye-coupled nucleotides respectively, and co-hybridized to cdna microarrays. following data acquisition and analysis, transcripts with at least a fold change in abundance in response to 20e were identified. after a hour treatment with 20e, genes were induced. two-thirds of these genes including the known 20e-inducible genes eip br, eip28/ and eip remained significantly up-regulated after hours in the presence of 20e while additional genes were induced and genes were repressed. among the 20e-induced genes indentified here in kc cells, the most strongly induced after hours, eip28/ was originally identified owing to its quick reaction to 20e in kc cells and has since been used to study aspects of the 20e response such as its tissue specificity. taken together, the response observed here after a and hour 20e-treatment are consistent with two prior observations on the 20e-induced puffing patterns of larval salivary glands: first, that primary 20e-response genes are likely to still be detectable after a four-hour 20e-treatment and second, that the 20e-response of kc is more robust after hours of 20e exposure. of the 20e-inducible genes identified here, 30% were also identified in a recent study by gauhar and colleagues that examined the transcriptional response of kc cells following either a or hr 20e treatment. however, if we apply more stringent fold-change criteria to the list of 20e indentified here the degree of overlap increases such that at a 2-fold cut-off, there is 50% overlap, and there is 100% overlap between ours and the gauhar et al. lists if we select a 3-fold change in expression. remaining differences in the identification of 20e responsive genes indentified is likely due in part to the use of different culturing medium as we have previously found that media composition affects the transcriptional response of kc cells to 20e. the use of a different microarray platform in the measurement of the transcriptional profile could also be a contributing factor. identification of primary 20e-response genes having identified genes involved in the 20e-response, we were interested in determining where these genes fit into the 20e-hierarchy of gene activation. in particular, we wanted to identify those genes comprising the primary response to the hormone and are thus likely to be directly induced by 20e. maximum induction of primary-response genes by 20e is expected to occur after hours of exposure to the hormone, however, some secondary-response gene transcription may have already begun at this point. one distinguishing characteristic of primary and secondary-response genes is their dependence on protein synthesis for transcription; only primary response genes are transcribed in the absence of protein synthesis. thus, to identify primary response genes, cells were treated with 20e in conjunction with an inhibitor of protein synthesis, cycloheximide. rna isolated from treated cells was labelled and co-hybridized to cdna arrays along with a differentially labelled untreated control. from this analysis we identified genes that were induced and genes that were repressed by the 20e-cycloheximide combination. to determine which of these genes were induced/repressed in response to 20e, and were not responding to exposure to cycloheximide alone, additional microarrays were run using rna from kc cells treated only with cycloheximide. from this experiment, genes were indentified whose transcription is altered due solely to the inhibition of protein synthesis. a two-class sam analysis identified genes whose transcripts levels in the presence of cycloheximide differs significantly when 20e is present including the well characterized 20e-primary-response genes, eip75b and br. although eip74ef was not among the genes on this list, we confirmed its induction by 20e by qrt-pcr suggesting that the cdna probe on the batch of arrays used for this study was likely of poor quality and unable to detect transcripts of this gene. of the primary-response genes identified, were not induced by 20e alone. the induction of these genes by 20e is likely masked in the absence of cycloheximide due to a negatively regulated auto-feedback loop where primary-response gene transcription is repressed by primary-response gene protein products. interestingly, the only gene transcript repressed as part of the primary response to 20e, cg also evades detection unless protein synthesis is blocked. the majority of 20e-responsive genes identified after a hour treatment with μm 20e behaved as secondary-response genes as their response to 20e was blocked in the presence of cycloheximide; genes were induced and were repressed as part of the secondary response. a two- class sam analysis of triplicate data identified primary-response genes that responded to 20e in the presence of cycloheximide. the top up-regulated primary-response genes and single down-regulated primary-response gene are listed. secondary-response genes regulated by 20e only in the absence of cycloheximide were also identified. the top up- and down-regulated secondary-response genes are listed. the fold change of each gene following a hour treatment with μm 20e and, for primary-response genes, the fold change following a hour treatment with μm 20e and μm cycloheximide is given. q-values determined by sam represent the lowest false discovery rate at which that gene is considered significant. for a gene to be included as differentially expressed, it must have a q-value < 5% in that experimental condition. q-values > 5% are not shown. david was used to identify enriched annotation categories among groups of genes that respond similarly to 20e. some of the most enriched functional terms associated with each group are given in the right most column. a complete list of all genes in these groups and their functional annotation is in additional file in addition to differences in transcriptional response in the absence of protein synthesis, it is also expected that most primary-response genes encode regulators while secondary-response genes encode effectors. to determine if this is the case for the genes identified here as part of the primary and secondary response, we examined each group of genes for enrichment in functional annotation terms/keywords using the online resource david. the results of this analysis are summarized in table several of the primary-response gene transcripts include those that code for the well known dna binding transcription factors eip75b and br as well as vrille which is consistent with the model that primary-response proteins are required for the transcriptional induction of the secondary-response genes. "induction of programmed cell death by hormones" is among the most highly enriched terms associated with primary-response genes and comes as no surprise due to the well established involvement of 20e-regualted genes such as br and nc in cell death pathways during drosophila development. however, we should point out that that kc cells do not undergo apoptotic cell death in response to ecdysone treatment and therefore other factors and/or conditions must need to be present in order to carry out this particular process. "catalytic activity" and "hydrolase activity" were both also among the most significantly enriched terms supporting the idea that primary-response genes are involved in a broad range of regulatory roles. for example, one of the genes belonging to the hydrolase class includes imaginal disk growth factor a protein that when bound to its receptor initiates signal transduction cascades important to imaginal disk development. some primary- and secondary-response genes include genes that code for proteins involved in cell movement and organization and/or are associated with the cytoskeleton. these include up-regulated transcripts for roadblock, a gene coding for a dynein-associated protein and down-regulated transcripts for spastin, and a gene coding for a microtubule severing protein. transcript changes in genes of this sort are consistent with the cell movement and morphogenetic changes that occur during 20e-dependent developmental changes. also of interest is that several of the secondary-response genes include genes that are involved in metabolic processes. more specifically, there seems to be an increase in transcripts for genes involved in mitochondrial respiration and a decrease in transcripts for genes involved in glycolysis perhaps suggesting that in addition to the cellular organization and tissue changes that are induced by 20e, there is also a shift in how energy is being produced. transcriptional response of kc cells to an increased concentration of 20e and to its structural analog, ponasterone a to examine the specificity of the 20e-response of kc cells to both hormone concentration and ligand structure, two other conditions were tested: a 20-fold higher concentration of 20e, μm, and the replacement of 20e with its structural analog, poa. 20e concentration is known to affect at least three aspects of the 20e-response: the size of both early and late puffs, the rate of early puff regression, and the transcription of some early gene isoforms. poa, often used as substitute for 20e, is a more potent activator of the 20e-response. at least four different measures indicate that poa activity is approximately eight fold higher than 20e activity: affinity for the ecdysone receptor, morphological changes associated with the 20e-response, the level of induction of known 20e-inducible proteins, and the level of transcription from an ecdysone response element. thus, to minimize effects caused by the difference in receptor affinity and make the comparison between analogs as similar as possible, our working concentration of poa was eight fold less than that of 20e. genes that responded to either of the two treatments were identified by microarray analysis as described above. compared to the number of genes induced by a hour treatment with μm 20e, both μm 20e and μm poa stimulated the transcription of a greater number of genes: and respectively versus induced by μm 20e. furthermore, poa down-regulated substantially more genes than either 20e treatment. one-hundred and fifty genes were repressed after a hour treatment with poa - almost double the number repressed by μm 20e and nearly six fold more than the number repressed by μm 20e. it is possible that the differences in the response to poa may be due to the activation of additional nuclear receptors; it is only assumed to be the 20e receptor heterodimer because it has been shown to be the case in vitro. overall, genes were identified that responded to at least one of the three hormone treatments examined. to determine which of these hormone-responsive genes were sensitive to all three conditions tested and which are differentially regulated by either μm 20e or poa relative to μm 20e, a multi-class sam analysis was performed where the transcriptional response to μm 20e and poa was compared to the transcriptional response μm 20e for each gene. this analysis revealed two things about the transcriptional response to each of the hormone conditions tested. first, of the genes that responded to at least one of the three treatments, showed no significant difference in their transcriptional response to all three treatments and eight of these are primary 20e-inducible genes. second, poa produced a more distinct transcriptional profile than did an elevated concentration of 20e with respect to the μm 20e treatment, suggesting that poa may not be acting solely through the ecdysone hierarchy. we confirmed the transcriptional response of three key 20e-inducible genes, the ecdysone-receptor complex, and a novel 20e-incudible gene by qrt-pcr. in all cases, the transcriptional response for these genes measured by qrt-pcr, although differing in magnitude when compared to the microarray analysis, had the same direction of change. cg5346 pect prl-1 cg4825 eater cyp4e2 cg14523 granny-smith a multi-class sam analysis revealed genes with similar transcriptional responses to the hormones and hormone concentrations studied. genes identified as primary 20e-response genes are indicated in bold. as with table the fold change and q-value of each gene following the treatment indicated is given. a missing fold change value means that gene was either not represented by at least 2/ spots and/or did not exhibit a coefficient of variation of at least for that treatment across replicate arrays. given the observed similarities and differences in the ecdysteroid-response to both an increased 20e-concentration and the use of a structural analog, we wanted to identify functions associated with genes that responded similarly to all treatments tested. using david we found two categories enriched in similarly responding genes that were also enriched among primary-response genes: "catalytic activity" and "hydrolase activity". conversely, genes that reacted significantly different to either μm 20e, poa or μm 20e were enriched in many unique categories; the most strongly enriched functional terms are associated with genes down-regulated by poa relative to μm 20e. furthermore, other categories related to development are also strongly enriched in the list of genes that are down-regulated by poa relative to μm 20e such as "anatomical structure development", "induction of programmed cell death", and "instar larval or pupal development". the apparent lack of activation of developmental genes by poa, as indicated by this analysis, is suggestive of an alternate mode of action of poa induced death in non-adapted insects more so than the simple induction of premature molting via the canonical ecdysone hierarchy. in fact, the results potentially suggest that poa might interfere with normal 20e induced transcriptional changes and prevent the induction of some of the genes that are required for 20e dependent developmental processes. three conditions both poa expression differences are relative to the μm 20e dataset. david was used to identify enriched annotation categories among groups of genes that responded similarly to one or more of the hormone treatments examined in this study. the modified fisher exact p-value represents the probability that the number of genes associated with an annotation term in the sample list is random compared to the occurrence of the term in the all of the genes on the array used for the analysis ). biological process; molecular function; cellular component; swiss prot keywords; interpro name; uniprot sequence feature. transcriptional response of salivary glands to 20e we wanted to extend our analysis to a 20e-responsive larval tissue to see how it compares to the 20e transcriptional response of kc cells. salivary glands were chosen in this respect due to the relative ease of dissection and because the foundations of our understanding of the 20e response comes from the study of salivary glands. cultured glands from late third instar larvae staged by the blue gut method were treated with μm 20e for four hours. we chose to treat salivary glands prior to the natural 20e pulse rather than assay glands after the natural pulse to rule out potential confounding effects caused by other endogenous hormones and/or the presence of protein products of secondary-response genes within the 20e hierarchy. rna was isolated from 20e treated glands and untreated control glands, amplified, and analyzed by microarray experiments as described in the methods section. we identified up-regulated genes including several known 20e-inducible genes such as, eip75b, eip74ef, ecr, nc, and impl and down-regulated genes. in a study by beckstead and colleagues looking at 20e responsive genes in cultured larval organs genes were identified as 20e responsive. their treatments, however, utilized twice the concentration of 20e and lasted for an additional two hours. nevertheless, we were still able to identify genes in common, the majority of which showed similar changes in expression level. by focusing only on salivary glands in this study we were able to identify salivary gland specific 20e-inducible genes that were not identified in total organ preparations. a david analysis reveals that these salivary gland specific genes are enriched for protein binding. genes induced by 20e in salivary glands but not in kc cells were enriched in go biological processes such as "salivary gland histolysis", "programmed cell death" and various "developmental" processes as might be expected. "cellularization" and "transporter activity" were also among the highest scoring hits. conversely, genes that were up regulated only in kc cells showed enrichment for "metabolic process" and "oxidoreductase activity"- categories that were also enriched among secondary response genes in kc cells suggesting that the secondary response genes differ between kc cells and salivary glands. comparison of 20e responsive genes from both systems allowed the identification of genes whose response is conserved across both systems suggesting that they might occupy key positions within the hierarchy. among the most highly induced of these genes are two genes that have not been traditionally associated with the 20e response in drosophila, cg and rrp. cg is currently un-annotated and rrp is a component of the exosome that has been shown to relocate to developmental loci during periods of active transcription and is involved in mrna processing. rrp showed the greatest correlation in expression profile to eip75b highlighting the importance of gene transcription in the ecdysone response. david was used to identify enriched annotation categories among groups of genes that respond similarly to 20e. the modified fisher exact p-value represents the probability that the number of genes associated with an annotation term in the sample list is random compared to the occurrence of the term in the all of the genes on the array used for the analysis ). biological process; molecular function; cellular component; swiss prot keywords; interpro name; uniprot sequence feature. tissue specific differences in ecdysone receptor isoform induction by 20e some of the differences in the transcriptional response of kc cells and salivary glands may be explained by differential chromatin states of inactive genes in the two systems. for example, salivary gland secretion protein genes known be under the control of the hierarchy were expressed in glands but not in cells most likely because of the difference in the two cell types. differential gene regulation by 20e may also result from the expression of different ecr isoforms or combinations thereof. typically, tissues with dissimilar developmental fates express different ecr isoforms. our microarray results agree with previous findings that kc cells express hemocyte and plasmatocyte marker genes pxn, ush, and hml and that 20e treatment of kc cells causes the induction of eip28/ and eip, both of which are 20e-regulated in the hematopoietic lymph gland providing further support for an embryonic hemocyte origin of kc cells. in contrast to larval salivary glands, embryonic hemocytes persist in the adult fly suggesting that kc cells are likely more similar to imaginal disc or adult tissues than larval tissues and placing them in a different metamorphic class than the salivary gland. because ecr is a 20e-inducible gene, we first examined the overall level of induction of ecr by qrt-pcr using a primer designed against the 3' region of the gene common to all isoforms. interestingly, unlike many of the 20e-inducible genes we identified by microarray analysis, ecr was induced to a similar extent in poa-treated kc cells, 20e-treated kc cells and 20e-treated salivary glands. we next examined the transcriptional response of each ecr isoform to poa and 20e using primers that recognize unique regions of each isoform and found that there was no difference in isoform regulation in kc cells treated with either hormone. comparison of isoform induction in 20e-treated kc cells to 20e-treated glands, however, revealed a tissue specific pattern in isoform regulation. at least two ecr isoforms were induced in kc cells by both poa and 20e but not induced at all in salivary glands treated with 20e. it has previously been shown that the proteins encoded by the ecr-ra, ecr-rd and ecr-re isoforms show strong expression in both embryos and imaginal discs and ecr-a is thought to be responsible for adult differentiation ecr-rc, on the other hand, although induced by both poa and 20e in kc cells, was more strongly induced by 20e in salivary glands. although differences in ecr isoform regulation in kc cells and salivary glands by 20e may be linked to the variation in the observed transcriptional response to 20e, it is unlikely that the same explanation can be applied to account for the difference in response observed in kc cells treated with 20e and poa since qrt-pcr analysis reveals at best, subtle differences in ecr isoform regulation by 20e and poa. it is tempting to speculate that the different ecr isoforms form 20e receptors with different target genes and thus contribute to differences in the 20e induced transcript profiles. alternatively, the differences in transcript profiles may have very little to with which ecr isoforms are present and more to do with cell type differences of other proteins. for example, differences in the concentrations of particular chromatin proteins could affect the accessibility of the various target genes differently in different cell types. also, the relative amount of proteins that partner with ecr to make the 20e heterodimer nuclear receptor could affect the sites the 20e receptor binds to. while the classical partner for ecr has been thought to be ultraspiracle, it is now known that biologically active 20e receptors can be made without usp. in mid-third instar larvae, the 20e dependent induction of the glue proteins requires ecr but not usp suggesting that ecr makes 20e receptors that consist of ecr and yet to be characterized partner protein. one approach to assess the relative contributions of each of the ecr isoforms would be to employ rnai knockdown strategies of the different ecr isoforms in both kc cells and larval tissues and examine what affect this has on the 20e induced transcript profiles.  we have identified primary 20e-response genes that are induced by 20e in the absence of protein synthesis in drosophila kc cells, one of the most widely used cell lines employed to study the ecdysone response. the primary 20e responsive genes were enriched for gene ontology terms such as induction of programmed cell death by hormones, catalytic activity, hydrolase activity, and cytochrome p we have also identified secondary 20e-response genes. go terms that are enriched in the secondary-response genes include an increase in transcripts for genes involved in mitochondrial respiration and a decrease in transcripts for genes involved in glycolysis suggesting that 20e induces shifts in cellular metabolism. some primary- and secondary-response genes include genes that code for proteins involved in cell movement and organization and/or are cytoskeletal associated which is consistent with the cell movement and morphogenetic changes that occur during 20e-dependent developmental changes. comparison of the genome-wide transcriptional response to 20e to its plant derived structural analog ponasterone a revealed a large difference in the transcriptional targets of these molecules. while these two compounds are structurally very similar, many more genes related to various aspects of development appear to be significantly induced by 20e than by poa. more specifically, the most strongly enriched functional terms are associated with genes down-regulated by poa relative to μm 20e. the genes not induced by poa include go terms such as developmental process, cell differentiation, gamete generation, anatomical structure development, organelle organization and biogenesis, neuron development, and induction of programmed cell death. the apparent lack of activation of developmental genes by poa suggests that poa induced death in non-adapted insects may be due to more than the just the simple induction of premature molting via the canonical ecdysone hierarchy. we also compared the 20e response in kc cells to that of a natural 20e target tissue where the function of 20e has been well described, the salivary glands of wandering 3rd instar larvae, and found little overlap in 20e-responsive genes. genes induced by 20e in salivary glands but not in kc cells were enriched in go biological processes such as "salivary gland histolysis", "programmed cell death" and various "developmental" processes as well as "cellularization" and "transporter activity". to help identify a potential mechanism to explain the difference in the transcriptional responses in kc cells and salivary glands, we analyzed the 20e-induced transcription of the various ecr isoforms, a known 20e-inducible gene and a component of the nuclear receptor complex that binds 20e. we did find differences in the induction of ecr isoforms ecr-ra, er-rc, and ecr-rd/e between kc cells and salivary glands in response to 20e. this suggests that the relative amount of various ecr isoforms present in a cell is different in different cell types and possibly contributes to the transcriptional response a given tissue has to 20e. 
the 20-hydroxyecdysone hierarchy of gene activation serves as an attractive model system for studying the mode of steroid hormone regulated gene expression and development. we also find that many genes regulated by 20e in kc cells are not regulated by 20e in salivary glands of wandering 3rd instar larvae and we show that 20e-induced levels of ecr isoforms ecr-ra, er-rc, and ecr-rd/e differ between kc cells and salivary glands suggesting a possible cause for the observed differences in 20e-regulated gene transcription between the two cell types. in addition, we reveal large differences between kc cells and salivary glands with regard to their genome-wide transcriptional response to 20e and show that the level of induction of certain ecr isoforms differ between kc cells and salivary glands. we hypothesize that the differences in the transcriptional response may in part be due to differences in the ecr isoforms present in different cell types. furthermore, in drosophila different cell types elicit different morphological responses to 20e yet the cell type specificity of the 20e transcriptional response has not been examined on a genome-wide scale. we aim to characterize the transcriptional response to 20e and poa in drosophila kc cells and to 20e in salivary glands and provide a robust comparison of genes involved in each response. our genome-wide microarray analysis of kc cells treated with 20e or poa revealed that far more genes are regulated by poa than by 20e and that there is very little overlap between the transcriptional responses to each hormone. we report significant differences in the transcriptional responses of 20e and poa, two steroid hormones that differ by only a single hydroxyl group. poa has a higher affinity for the 20e nuclear receptor, composed of the ecysone receptor and ultraspiracle proteins, than 20e and a comparison of the genes regulated by these hormones has not been performed.
 the 20-hydroxyecdysone hierarchy of gene activation serves as an attractive model system for studying the mode of steroid hormone regulated gene expression and development. many structural analogs of 20e exist in nature and among them the plant-derived ponasterone a is the most potent. poa has a higher affinity for the 20e nuclear receptor, composed of the ecysone receptor and ultraspiracle proteins, than 20e and a comparison of the genes regulated by these hormones has not been performed. furthermore, in drosophila different cell types elicit different morphological responses to 20e yet the cell type specificity of the 20e transcriptional response has not been examined on a genome-wide scale. we aim to characterize the transcriptional response to 20e and poa in drosophila kc cells and to 20e in salivary glands and provide a robust comparison of genes involved in each response.  our genome-wide microarray analysis of kc cells treated with 20e or poa revealed that far more genes are regulated by poa than by 20e and that there is very little overlap between the transcriptional responses to each hormone. interestingly, genes induced by 20e relative to poa are enriched in functions related to development. we also find that many genes regulated by 20e in kc cells are not regulated by 20e in salivary glands of wandering 3rd instar larvae and we show that 20e-induced levels of ecr isoforms ecr-ra, er-rc, and ecr-rd/e differ between kc cells and salivary glands suggesting a possible cause for the observed differences in 20e-regulated gene transcription between the two cell types.  we report significant differences in the transcriptional responses of 20e and poa, two steroid hormones that differ by only a single hydroxyl group. we also provide evidence that suggests that poa induced death of non-adapted insects may be related to poa regulating different set of genes when compared to 20e. in addition, we reveal large differences between kc cells and salivary glands with regard to their genome-wide transcriptional response to 20e and show that the level of induction of certain ecr isoforms differ between kc cells and salivary glands. we hypothesize that the differences in the transcriptional response may in part be due to differences in the ecr isoforms present in different cell types. 