Transcriptional signature of an adult brain tumor in Drosophila

Background

Mutations and gene expression alterations in brain tumors have been extensively investigated, however the causes of brain tumorigenesis are largely unknown. Animal models are necessary to correlate altered transcriptional activity and tumor phenotype and to better understand how these alterations cause malignant growth. In order to gain insights into the in vivo transcriptional activity associated with a brain tumor, we carried out genome-wide microarray expression analyses of an adult brain tumor in Drosophila caused by homozygous mutation in the tumor suppressor gene brain tumor (brat).

Results

Two independent genome-wide gene expression studies using two different oligonucleotide microarray platforms were used to compare the transcriptome of adult wildtype flies with mutants displaying the adult brat^k06028mutant brain tumor. Cross-validation and stringent statistical criteria identified a core transcriptional signature of brat^k06028neoplastic tissue. We find significant expression level changes for 321 annotated genes associated with the adult neoplastic brat^k06028tissue indicating elevated and aberrant metabolic and cell cycle activity, upregulation of the basal transcriptional machinery, as well as elevated and aberrant activity of ribosome synthesis and translation control. One fifth of these genes show homology to known mammalian genes involved in cancer formation.

Conclusion

Our results identify for the first time the genome-wide transcriptional alterations associated with an adult brain tumor in Drosophila and reveal insights into the possible mechanisms of tumor formation caused by homozygous mutation of the translational repressor brat.

Cancer is a multistep process, which involves loss of cell proliferation control, resistance to cell death, and invasion as well as metastasis. This process is often associated with multiple and progressive genetic alterations. Proto-oncogenes are activated in a dominant fashion by mutation, chromosome translocation or gene amplification, whereas tumor suppressor genes are inactivated or lost by mutation, chromosome loss, mitotic recombination or gene conversion [1–3].

Mutations and gene expression alterations in brain tumors have been extensively investigated. For example, data has been accumulated for the cerebellar tumor medulloblastoma which is the most common malignant brain tumor in children [4]. In addition to chromosome loss, medulloblastomas overexpress certain genes, including c-myc, pax5, and zic, all of which encode transcription factors [5–7]. Thus, animal models of brain tumors are required to gain insights into the correlation of genetic alterations and tumor phenotype and to better understand how genetic alterations cause malignant neoplasm [4, 8]. To better understand the genetic alterations associated with brain tumors in vivo, the fruitfly Drosophila melanogaster is an excellent model system for which unprecedented genetic and transgenic technologies as well as entire genome sequence information are available [9, 10]. Indeed, more than two thirds of the known human cancer genes are thought to have Drosophila homologs [11, 12], and various studies have shown that genetic inactivation or misregulation of Drosophila oncogenes or tumor suppressor genes lead to neoplasms that display characteristic features of malignant growth like in human cancer [13–15]. The level of gene and pathway conservation, the similarity of cellular processes, and the emerging evidence of functional conservation of Drosophila genes that are homologous to mammalian oncogenes or tumor suppressor genes, indicate that studies of tumorigenesis in flies can directly contribute to the understanding of human brain tumors [16, 17].

One of the earlier studies on mutations causing malignant neoplasms in Drosophila examined 12 recessive lethal mutations affecting the optic neuroblasts and ganglion mother cells in the larval brain, the imaginal discs, and the hematopoietic organs [13]. Among others, it could be shown, that a mutation in the lethal (2) giant larvae (lgl) gene leads to a failure of the optic neuroblasts of the inner and outer formation centers and the ganglion mother cells (GMCs) to generate adult neurons. The neuroblasts and GMCs proliferated extensively, causing a dramatic enlargement of the larval brain hemispheres. Upon transplantation of neoplastic tissue into wildtype hosts, the transplanted tissues invaded the body cavities, the abdomen, the thorax, the ovaries, the gut, and the thoracic muscles [13].

Further investigation of mutations in the dlg1 and lgl genes revealed that the resulting phenotypes exhibited most of the characteristics defining neoplastic growth in vertebrates, including rapid growth both in situ and in culture, defective intercellular junctions, low adhesiveness in between cells, defective intercellular communication, absence of terminal differentiation, lack of response to ecdysone, invasiveness, and lethality to the host ([13], reviewed in [18]).

Furthermore, it is known that expression of an activated form of RAS1 (RAS1V12) in Drosophila imaginal discs is leading to ectopic cell proliferation and hyperplastic tissue growth [19]. The mutational activation of the human homolog of Ras is associated with a wide variety of human tumors (reviewed in [20]). When RASV12 expression is combined with inactivation of scribbled. (scrib), lethal (2) giant larvae (lgl), discs large 1 (dlg1), bazooka (baz), stardust(sdt) and Cdc42 in Drosophila larvae, migration, invasion of tumor tissue and secondary tumor foundation can be found [15]. Another Drosophila gene possessing tumor suppressor properties is warts (wts). Flies homozygous for strong wts alleles show overproliferation in the imaginal discs and the central nervous system leading to pupal lethality. The observed cellular clones where irregular in shape, unpatterned and formed folded lobes. The cells to be found in wts mutant clones were larger than neighboring wildtype cells [21]. The mammalian homolog of wts, LATS1, is also known for causing soft-tissue sarcomas and ovarian stromal cell tumors when mutated in mice [22]. The current model of wts / LATS1 function suggests interaction with CyclinA / Cdc2, linking it directly to the cell cycle ([23], reviewed in [24]).

Mutations in the Drosophila gene salvador (sav) have been associated with increased cell growth and an increase in cell number. Also, Sav has been suggested to interact directly with Wts and thereby contributing to its function [25].

In this report, we have determined at the full genome level the transcriptional signature of an adult brain tumor in Drosophila caused by homozygous mutation in the tumor suppressor gene brain tumor (brat). This gene encodes a member of the conserved NHL family of proteins, which appear to regulate differentiation and growth. These factors have been directly implicated in human tumor formation [26, 27]. Based on the analysis of loss- and gain-of-function experiments, Brat acts in translation repression, ribosomal RNA synthesis, and negative regulation of cell growth [28–30]. Structural analysis indicates that the NHL beta-propeller is the essential domain for both the translational repression and cell growth inhibitory activities of Brat [31]. Inactivation of brat results in neoplastic overgrowth and tumor formation in the larval brain and generally causes lethality in the larval third instar and pupal stages [32]. It has been hypothesized that in brat loss-of-function larvae the optic neuroblasts and ganglion mother cells are incapable of generating neurons and repesent the source of the brain tumor [33]. Tumor cells derived from homozygous brat larval brain tissue can grow not only in situ, but also after transplantation into adult host flies. The transplanted cells grow rapidly, forming metastases and secondary malignant tumors, and finally kill their host [34].

Here we focus on the strong neoplastic adult brain phenotype caused by homozygous mutation of the brat^k06028allele. To analyze the transcriptional activity associated with this brain tumor, we carried out two independent genome-wide microarray expression studies using two different oligonucleotide array platforms to compare adult wildtype flies with flies displaying the adult brat^k06028mutant brain tumor. Our experiments identified 321 annotated genes that showed highly significant (p<=0.0001) changes in expression levels due to loss of the tumor suppressor activity of brat. One fifth of these genes show homology to known mammalian genes involved in cancer formation. We observed significant upregulation of genes involved in the control of asymmetric neuroblast division. We also found significant upregulation of genes involved in ribosome biogenesis, translation, and RNA processing, suggesting a link between the tumor suppressor activity of brat and its role as a translational repressor during Drosophila development.

Adult homozygous brat^k06028mutants show a strong neoplastic brain phenotype

The brat mutant allele brat^k06028was generated by a Plac W transposon insertion in the non-coding Exon 4 of the brat locus lying immediately upstream of the transcriptional start site [28, 32]. Larvae homozygous for this insertion manifest a strong neoplastic brain phenotype. We found that 15% (n = 979) of homozygous brat^k06028mutants eclose. These homozygous adult brat^k06028mutants show limited motility compared to Oregon R wildtype flies and their lifespan is strongly reduced; most of the homozygous brat^k06028mutants die after the first 10 days. Histological analysis showed that homozygous adult brat^k06028mutants display a pronounced overproliferation in their brains as compared to Oregon R wildtype brains (Figure 1). The penetrance of this phenotype is 100% (n = 300). Most of the newly hatched brat^k06028mutants showed a head capsule completely filled with neoplastic tissue. In addition, we observed ectopic neuropil-like structures; these structures appeared bilaterally in ventro-medial positions dorsal to the central brain complex and more ventro-laterally to the optic lobes (Figure 1B, arrowheads). In paraffin sections the size of the central brain structures and the optic lobes of brat^k06028mutants appeared reduced as compared to the wildtype, however, the overall morphology of neuropil subdivisions in these brains was retained.

Adult flies homozygous for the brat^k06028allele show a strong neoplastic brain phenotype. Horizontal 7 μm paraffin sections at the level of the oesophagus (oe) of the heads of two days old adult Oregon R wildtype (A) and brat^k06028(B) flies. Homozygous adult brat^k06028mutants display a pronounced neoplastic phenotype in their heads (B, arrow) and ectopic neuropil-like structures can be observed (B, arrowheads). The size of several neuropil structures like lamina (la), medulla (me), lobula (lo), lobula plate (lp), antennal lobe (al), and ventrolateral protocerebrum (vp) appears reduced (B, compare to wildtype in A). However, their overall morphology appears to be unaffected, suggesting that the cells giving rise to these neuropil structures differentiated normally in homozygous adult brat^k06028mutants.

Full size image

Transcriptional signature of the adult brat^k06028mutant brain tumor

In order to gain insights into the genetic alterations associated with a fully developed brain tumor in vivo, we determined at the full genome level the transcriptome of the neoplastic brain tumors caused by homozygous brat^k06028mutation as compared to wildtype brain. To ensure cross-validation and increased significance of our results, we conducted two independent genome-wide gene expression studies to compare wildtype flies with homozygous brat^k06028flies using two different oligonucleotide array platforms (Table 1). In the first experiment, RNA was extracted from homozygous adult brat^k06028fly heads and from isolated adult fly heads of Oregon R as control. The labeled cRNA of this experiment was hybridized to custom made full-genome GeneChips (roDROMEGa, experiment A). In the second experiment, RNA was extracted from dissected adult brains, and hybridization of labeled cRNA to commercially available full-genome Gene Chips involved a signal amplification step (DrosGenome1, experiment B; see materials and methods). Moreover, in this second experiment, the transcriptome of dissected adult brains of homozygous brat^k06028mutants was compared to dissected adult brains of flies generated by transposon excision of the brat^k06028P-element (termed brat^k06028 jumpout, see materials and methods). Precise transposon excision resulted in reversion of the brat^k06028neoplastic brain phenotype to wildtype-like brain (data not shown).

Stringent quality control and filtering for the two resulting data sets was done independently (see materials and methods). After analysis, 1778 transcripts in the wildtype condition and 2955 transcripts in the brat^k06028condition of experiment A were statistically judged as present. This represents 12.7% and 21% of all transcripts on the roDROMEGa array, respectively. In experiment B, 5063 transcripts were judged as present in the wildtype-like brat^k06028 jumpout condition and 4981 in the brat^k06028condition. This represents 36.3% and 35.7% of all transcripts on the DrosGenome1 array, respectively. Transcripts in both datasets were considered as differentially expressed when the change between conditions was larger than 2-fold, their signal strength was above 10 in the condition with the higher expression level, and the change had a significance value of p ≤ 0.0001 (t-test). In experiment A 725 transcripts were judged as differentially expressed, whereas in experiment B 1888 transcripts were judged as differentially expressed in the brat^k06028tumor condition (Table 1).

The core set of genes differentially expressed in brat^k06028neoplasms

To strengthen the significance of our data, we considered those genes for further analysis that passed the filter criteria in both experiments. This was the case for 321 transcripts, representing a core transcriptional signature that characterizes the neoplastic tissue of adult homozygous brat^k06028mutants in a highly reproducible manner (Table 1). In both experiments, the majority of the 321 differentially regulated genes (279) showed an increased expression level in neoplastic adult brat^k06028mutants. 42 transcripts out of 321 differentially regulated genes showed a decreased expression level in the brat^k06028condition. (5 transcripts out of the 321 genes showed an opposite differential expression change in experiment A and B). As shown in Figure 2, the majority of the 321 genes differentially regulated in brat^k06028neoplasms displayed comparable expression level changes in both experiments. Only 88 genes showed significantly higher expression level changes in experiment B as compared to experiment A. This may be due to the fact that in experiment B cRNA hybridization to the array involved a signal amplification step (see materials and methods).

Expression level changes observed for the 321 genes identified in both oligonucleotide array experiments as differentially regulated in the brat tumor condition. X-coordinate shows relative fold-changes for each individual gene (y-coordinate) depicted in red for experiment A and depicted in orange for experiment B, respectively. Note that the two experiments A and B disagree on the differential expression level changes only for five transcripts. Please note that the highest expression levels are cut at a value of 80 due to the scaling of the y-axis.

Full size image

The 321 genes differentially expressed in both experiments were grouped into functional classes according to their annotation in Flybase. The function of most of the differentially expressed transcripts were unknown (n = 67). A list of these genes is shown in Figure 8. The remaining genes of known function were grouped into the functional classes metabolism (n = 56), cytoskeleton/structural proteins (n = 31), RNA binding/processing (n = 27), transcription/replication/repair (n = 20), translation (n = 19), cell cycle (n = 18), transcriptional regulation (n = 17) chromatin structure (n = 13), signal transduction (n = 12), chaperones (n = 11), transport (n = 10), proteolytic systems (n = 8), stress response (n = 7), cell surface receptors/CAMs/ion channels (n = 4), and apoptosis (n = 1). You can find a list of these 321 genes along with the fold change between experiments in Supplementary Table 3 [see Additional file 3].

Genes involved in metabolism, cell cycle and apoptosis

Neoplastic overgrowth involves essential alterations in cell physiology, loss of cell proliferation control and resistance to cell death. Accordingly, we expected to observe the differential expression of genes linked to metabolism, cell cycle regulation and apoptosis in our gene expression profile studies. Indeed, transcript profiling of brat^k06028mutants identified transcripts of this type of which misregulation appears to be associated with malignant growth.

Among the genes involved in metabolism, we identified 41 transcripts with elevated expression levels (Figure 3). These genes included eyes absent (eya), whose activity is required for the survival of progenitor cells at a critical stage in eye morphogenesis [35] as well as uninitiated (und), the Drosophila methionine aminopeptidase 2, which was shown to be a putative regulator of translation initiation required for cell growth and patterning [36]. We also detected selenide, water dikinase (SelD) as differentially upregulated in brat^k06028neoplasms. SelD encodes a product involved in selenocysteine biosynthesis, and mutational inactivation in SelD results in a reduction in cell proliferation in the imaginal discs and the larval brain [37]. Within the functional class metabolism, we also found 15 transcripts as significantly down-regulated in the brat^k06028tumor condition (Figure 3). Among these were bubblegum (bgm), encoding a product with long-chain-fatty-acid-CoA-ligase activity. Mutant analysis of bgm revealed that in young homozygous mutant flies the optic lobes appear normal however as the flies grow older the optic lobes show signs of regional degeneration [38].

Differential expression of genes belonging to the functional classes metabolism, cell cycle and apoptosis. Fold-changes for each gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

As expected, all of the transcripts positively regulating cell cycle showed increased expression levels in the tumor condition as compared to the wildtype situation. Among these we identified cyclin A, cyclin B, cyclin E, string (stg), cdc2, Cdk7, Cks, Set, abnormal spindle (asp), fizzy (fzy), polo, and Myb oncogene-like (Myb). These genes are involved in the regulation of highly conserved aspects of cell proliferation such as progression through G1, S and G2/M phases, spindle orientation and the maintenance of genomic integrity throughout mitosis.

Of particular interest is the highly significant upregulation of the genes asp, fzy, rough deal (rod), and Myb. asp encodes a microtubule-associated protein that associates with the polar regions of the mitotic spindle. Mutants of asp are larval lethal, with a high frequency of aberrant (e.g. polyploid) cells arrested in metaphase in the larval brain, suggesting that asp may play a role in spindle pole organization during mitosis [39]. fzy encodes a product involved in cyclin catabolism and fzy mutants also show metaphase arrest with compact condensed chromosomes like asp mutants [40]. rod encodes a product involved in mitotic chromosome segregation which is localised to the kinetochore, and mutations in rod result in mitotic segregational failure due to delayed or incomplete release of sister chromatids [41]. Finally Myb encodes a proto-oncogene with transcriptional activator activity involved in centrosome cycle which is required to sustain the appropriate rate of proliferation, to suppress formation of supernumerary centrosomes, and to maintain genomic integrity [42].

It is noteworthy that the activity of all of these genes involved in cell cycle regulation has been linked to several aspects of neurogenesis [43, 44]. Thus, all of these genes are known to be expressed during embryonic and larval development in mitotically active cells of the nervous system, the neuroblasts and ganglion mother cells (GMCs), suggesting that their elevated levels of transcriptional activity correlate with an elevated and aberrant activity of the cell cycle machinery in neoplasms of homozygous brat^k06028mutants.

We did not detect any significant alteration of gene expression concerning transcripts belonging to the functional class apoptosis, except elevated expression levels for thioredoxin peroxidase 2 (Jafrac2), which is involved in the induction of apoptosis [45]. The fact that we only identified one transcript involved in programmed cell death suggests that the apoptotic program is either impaired below detectability or is unaltered in neoplastic tissue of adult brat^k06028mutant brain tumors.

Genes involved in translation and RNA binding/processing

Recent studies have shown that brat acts as a translational repressor [29]. Moreover, mutant analyses demonstrated that brat-mutant cells are larger than wildtype cells and have enlarged nucleoli, which are associated with an increase in total rRNA production. Also, the C. elegans homolog of brat, ncl-1, is known for being a repressor of RNA polymerase I and III transcription and for being an inhibitor of cell growth. Loss of function mutations in ncl-1, result in enlarged nucleoli. The rates of rRNA and 5S RNA transcription are increased and cells are enlarged. [46]. These data suggest that excess ribosomal synthesis and cell growth may be important aspects of the tumorous phenotype of brat^k06028mutants [30]. Indeed, in addition to an elevated and aberrant activity of the cell cycle machinery, our gene expression profile of brat^k06028adult brain tumors detected genes involved in all aspects of ribosome synthesis and translation control. Thus, transcript profiling of brat^k06028mutants identified 19 transcripts as differentially expressed belonging to the functional class translation and 27 transcripts as differentially expressed in the functional class RNA binding/processing. All of these transcripts showed increased expression level in the tumor condition as compared to the wildtype situation (Figure 4).

Differential expression of genes belonging to the functional classes translation and RNA binding / processing. Fold-changes for each individual gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

Among the transcripts that belong to the functional class translation we identified Eukaryotic initiation factor 1A (eIF-1A), Eukaryotic initiation factor 3 p40 subunit (eIF-3p40), and Int6 homologue (Int6) that encode proteins with translation initiation factor activity involved in protein biosynthesis [47]. Moreover, genes involved in translation elongation factor activity like Elongation factor 1α48D (Ef1α48), Ef1 γ, and eEF1 δ were detected as significantly upregulated in the tumor condition. In addition, genes involved in 35S primary transcript processing (Fibrillarin, Fib) and in rRNA processing (Nop5) were identified as upregulated. Interestingly, we also detected upregulation of Ribosomal protein L1 (RpL) and Ribosomal protein L3 (RpL3), both of which encode structural constituent of the ribosome involved in protein biosynthesis. RpL3 is of particular interest as it has been mapped to the chromosomal region harbouring the Minute(3) 86D mutation [48].

Among the transcripts that belong to the functional class RNA binding/processing we identified Rrp4 and Csl4, both encoding gene products with 3'-5' exoribonuclease activity involved in mRNA processing [49]. Also Developmental embryonic B (DebB) and SC35 with pre-mRNA splicing factor activity involved in mRNA splicing were detected as significantly upregulated [50]. In addition, Rm62 encoding a product with ATP dependent RNA helicase activity involved in RNA interference was upregulated in the tumor condition, as were U2 small nuclear riboprotein auxiliary factor 50 (U2af50) and U2 small nuclear riboprotein auxiliary factor 38 (U2af38), both encoding products with pre-mRNA splicing factor activity involved in mRNA splicing (reviewed in [47]). Within this functional class we also identified upregulated transcripts which are involved in various aspects of cell growth and proliferation, such as polyA-binding protein (pAbp) which encodes a product with poly(A) binding involved in positive regulation of translation. pAbp mutants display cytokinesis defects, and mutant analyses demonstrated that pAbp is required to connect the centrosome to the minus-ends of free microtubules (reviewed in [47]). Intriguingly, we also detected Small ribonucleoprotein Sm D3 (SmD3), which encodes a protein involved in mRNA splicing, as upregulated in the brain tumor condition. This is of particular interest since homozygous 12-day old SmD3 mutant larvae show overgrowth of the brain, hematopoietic organs and imaginal discs and die during third instar larval or pupal stages [51]. Finally, we also identified staufen (stau), which has been shown to be involved in the asymmetric localization of cell fate determinants during neuroblast divisions [52], as upregulated in neoplastic tissue.

Genes involved in chromatin structure, transcription/replication/repair, and transcriptional regulation

In addition to the expected upregulation of genes involved in cell cycle regulation as well as ribosome synthesis and RNA processing, we also observed that the basal transcriptional machinery was significantly upregulated in the brat^k06028tumor condition (Figure 5). Within the functional class chromatin structure we found 13 genes with elevated expression levels. Among those are genes involved in various aspects of chromatin modeling such as Nucleoplasmin (Nlp), involved in nucleosome spacing, and Nucleosome assembly protein 1 (Nap1) as well as Chromatin assembly factor 1 subunit (Caf1), both encoding proteins with histone-specific chaperone activity [53]. We also detected two suppressors of variegation, namely Suppressor of variegation 3–9 (Su(var)3-9) with histone methyltransferase activity [54] and Suppressor of variegation 205 (Su(var)205) involved in the establishment of chromatin silencing. Interestingly, mutant analysis of Su(var)205 revealed that larval brain neuroblasts show cytokinesis defects during mitosis [54]; this is also the case for gluon (glu), which we detected as significantly upregulated in the tumor condition. In glu mutants, larval brain neuroblasts show distinct abnormalities during chromosome segregation and lethality occurs during the late larval stage [55].

Differential expression of genes belonging to the functional classes chromatin structure, transcription / replication / repair, and transcriptional regulation. Fold-changes for each gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

In the functional class transcription/replication/repair we identified 20 genes as differentially regulated in the brat^k06028tumor condition, all of them having elevated expression levels as compared to wildtype. Among them are the minichromosome maintenance (MCM) genes, which encode an evolutionary conserved family of molecules that form an important part of the pre-replicative complex, required for DNA replication [56]. We found two members of this complex (mcm2 and mcm3) as upregulated in brat^k06028neoplasms. Interestingly, Mcm2 mutants persist as third instar larvae for several days, and dissection of wandering larvae demonstrates that their CNS is smaller than that of wildtype and that they lack identifiable imaginal discs [57]. We also detected mutagen-sensitive 209 (mus209), involved in mitotic spindle assembly as upregulated. Mutants of mus209 show abnormal chromosome condensation [58], whereas mutations of rad50, a gene involved in DNA repair, result in elevated levels of anaphase bridges in dividing cells of third instar larval brain and imaginal discs. Moreover, RNA polymerase II 15 kD subunit (RpII15), involved in transcription from Pol II promoter [59] and absent, small or homeotic disc 2 (ash2), a member of the trithorax group were detected as upregulated. Interestingly, mutant alleles for ash2 are larval pupal lethals and display imaginal disc and brain abnormalities [60].

Among the transcripts belonging to the functional class transcriptional regulation, we identified two members of the snail family of zinc-finger transcription factors, namely worniu (wor) and snail (sna) as significantly upregulated in brat^k06028neoplasms. During wildtype embryogenesis, both genes are expressed during neuroblast delamination, and mutant analyses suggest that Snail and Worniu function in neuroblasts, around the time of division to give rise to ganglion mother cells (GMCs) [61]. Moreover, we also detected upregulation of several other transcription factors involved in the control of proliferative activity of embryonic and postembryonic neuroblasts. Thus, we identified tailless (tll), deadpan (dpn), and castor (cas) as upregulated. tll is expressed in procephalic neuroblasts and required for their specification as tll mutants lack the anterior brain [62]. dpn encodes a pan-neural gene shown to be expressed in all neuronal lineages during embryogenesis as well as in neuroblasts of the larval CNS and in precursors of sensory neurons in imaginal discs. Loss- and gain-of-function experiments suggest that dpn activity is critical for the proper regulation of cell proliferation in the larval brain [63]. cas is involved in a cascade of sequentially expressed transcription factors, which translates information on timing of GMC formation [64]. In addition, cas is known to be involved in post-embryonic brain development as small homozygous somatic mutant cas clones in the adult brain lead to ellipsoid body and mushroom body defects [65]. Finally, we also detected diminutive (dm), the Drosophila homolog of the proto-oncogene c-myc, as significantly upregulated in the brat^k06028tumor condition. Myc genes link patterning signals to cell division by regulating primary targets involved in cellular growth and metabolism [66].

Genes involved in cytoskeleton/structural and signal transduction

In previous studies, histological analysis of brat^k06028mutants indicated that the brain tumor phenotype is primarily due to uncontrolled divisions of optic lobe neuroblasts and GMCs [33]. Our microarray expression analysis supports the assumption that neuroblasts and/or GMCs are at the cellular origin of these brain tumors. This is not only evident by the significant upregulation of genes involved in neuroblast delamination and proliferation such as wor, sna, tll, dpn, and cas, but also by the elevated transcription levels of specific genes belonging to the functional classes cytoskeleton/structural and signal transduction (Figure 6).

Differential expression of genes belonging to the functional classes cytoskeletal / structural and signal transduction. Fold-changes for each gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

In the functional class cytoskeleton/structural, we found high levels of transcription in brat^k06028tumorous tissue for the genes miranda (mira) and inscuteable (insc). This is of particular interest, since both proteins play essential roles in the asymmetric localization of cell fate determinants during neuroblast divisions. Miranda creates intrinsic differences between sibling cells by mediating the asymmetric segregation of transcription factor Prospero into only one daughter cell during neural stem-cell division. In insc mutants, mitotic spindles in neuroblasts fail to rotate into proper position and neuroblasts divide in random orientation resulting in general disorganization of the neuroblast array and defective neuroblast morphology (reviewed in [67]). Moreover, we also identified as significantly upregulated in the tumor condition ciboulot (cib), an actin binding protein, and the Drosophila profilin protein chickadee (chic), both of which cooperate in central brain metamorphosis [68]. Intriguingly, expression of big brain (bib) appears to be downregulated in the brat^k06028tumor condition. bib encodes a product with connexon channel activity, which is required autonomously in epidermal precursors to prevent neural development (reviewed in [69]).

Among those transcripts belonging to the functional class signal transduction, we identified several members of GTPase signal transduction pathways as significantly upregulated in brat^k06028neoplasms. Thus, we identified RacGAP50C, which has GTPase activator activity and is involved in Rho protein signal transduction, and also three members of the Ran GTPase signal transduction pathway, namely Segregation distorter (Sd), ran, and moleskin (msk). Moleskin, for example, appears to be involved in EGF receptor signalling pathway initiation of gene expression in response to Drosophila receptor tyrosine kinase signalling [70]. We also detected MAP kinase activated protein-kinase-2 (MAPk-Ak2), which has protein serine/threonine kinase activity, as significantly upregulated. In addition, we identified Notch (N) as significantly upregulated in the brat^k06028tumor condition. Notch encodes a transmembrane receptor mediating cell-cell communication, and Notch signalling has been implicated in a wide variety of cellular processes, including the maintenance of stem cells, specification of cell fate, differentiation, and proliferation [71]. Surprisingly, however, the Insulin-related peptide (Ilp2), a component of the insulin signalling pathway, was downregulated in the tumor condition. Overexpression of Ilp2 alters growth control in a Insulin receptor-dependent manner, suggesting a role for Ilp2 in controlling organismal size by augmenting both cell number and cell size of different organs (reviewed in [66]).

Genes involved in cell surface receptors/CAMs, chaperones, proteolytic systems, stress response and transport

Previous studies have shown that tumor cells derived from homozygous brat larval brain tissue can grow not only in situ, but also after transplantation into adult host flies. The transplanted cells grow rapidly, forming metastases and secondary malignant tumors that finally kill their host [34]. This suggests that loss of brat leads to the acquisition of invasive and metastatic properties of the resulting tumor cells. This acquisition of invasive and metastatic abilities is associated with altered binding specificities of cell adhesion molecules (CAMs), the activation of proteases, as well as chaperones and heat shock proteins mediating stress response [1, 72]. In accordance with this, our transcriptional analysis of brat^k06028tumor cells identified genes as differentially regulated in the tumor condition belonging to the functional classes cell surface receptors/CAMs, chaperones, proteolytic systems, and stress response (Figure 7). Within the class cell surface receptors/CAMs we identified Neurotactin (Nrt), a transmembrane protein that localizes to cell-cell junctions where it mediates cell adhesion and cell signalling, as upregulated [73]. Three other genes belonging to the class cell surface receptors/CAMs were differentially downregulated in the brat^k06028tumor condition. Among the classes proteolytic systems, chaperones and stress response, we identified genes upregulated like Ubiquitin carboxy-terminal hydrolase (Uch), involved in protein deubiquitination, and granny smith with leucyl aminopeptidase activity. Also the chaperones T-complex Chaperonin 5 (Cct5) and Cctgamma with chaperonin ATPase activity involved in protein folding and the heat shock proteins Hsp27, Hsp26, and Hsc70-4 were identified as significantly upregulated in the brat^k06028tumor condition. Finally, we detected 10 transcripts belonging to the functional class transport as differentially regulated in neoplastic brat^k06028tissue. Among these is Pendulin (Pen), an adaptor molecule which is expressed in embryonic neuroblasts and the proliferating regions of the larval brain lobes and may be required for the normal transmission and function of proliferative signals in the cells [74].

Differential expression of genes belonging to the functional classes cell surface receptors /CAMs, chaperones, proteolytic systems, stress response, and transport. Fold-changes for each gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

Differential expression of genes belonging to the functional class unknown function. Fold-changes for each gene are shown in red for experiment A and in orange for experiment B, respectively.

Full size image

Vertebrate homologs of genes differentially regulated in brat^k06028neoplasms are involved in mammalian cancer formation

Our analysis of the 321 genes differentially regulated in brat^k06028neoplasms revealed genes involved in general aspects of tumorigenesis such as elevated and aberrant metabolic activity, elevated and aberrant activity of the cell cycle and basal transcriptional machinery, as well as elevated and aberrant activity of ribosome synthesis and translation control. Considering the level of gene and pathway conservation, the similarity of cellular processes, and the emerging evidence of functional conservation of Drosophila genes that are homologous to mammalian oncogenes or tumor suppressor genes [16, 17], we wondered whether our data might reveal any parallels to the genetic alterations associated with mammalian cancer formation. This is of particular interest since one of the two human homologs of brat, TRIM3, has been mapped to chromosome 11p15 [75, 76], a region that has been termed the "multiple tumor-associated chromosomal region 1" due to the presence of numerous cancer-related genes in this region [77, 78].

We therefore searched available databases in order to identify vertebrate homologs of the 321 genes differentially regulated in brat^k06028neoplasms and to determine their possible involvement in mammalian cancer formation. We identified 62 out of 321 Drosophila genes as having mammalian homologs, of which misregulation has been associated with various types of tumors (see Table 2). Mammalian homologs of genes involved in the regulation of conserved aspects of cell proliferation were found to be implicated in mammalian cancer formation. Among these are mammalian homologs of cyclins, Cdks, and genes involved in various other aspects of cytokinesis and the maintenance of genomic integrity. Thus, the string homolog Cdc25B, encoding a dual-specific phosphatase that mediates cell cycle progression by activating cyclin-dependent kinases, has been shown to possess oncogenic potential and is frequently overexpressed in human prostate cancer tissues [79]. Drosophila CG9993 encodes for a long-chain fatty acid transporter and its rat homolog (FATP) shows elevated transcript levels in hepatoma cell lines [80]. Homologs of Drosophila MCM family of proteins are known as proliferation markers in vertebrates, and evidence exists that HsMCM2 is upregulated in primary human tumor tissue [81]. Furthermore, it has been shown that the rate and level of Mcm3 and MCM4 expression appears to be higher in cancer cells than in normal proliferating cells of the uterine cervix and dysplastic cells.

Prominent examples of genes involved in mammalian cancer formation are the proto-oncogenes Myb and Myc. The Myc gene is a central regulator of proliferation, differentiation, cell survival, and neoplastic transformation. It is found mutated or overexpressed in up to 30% of human cancers and the Myc network appears to regulate a large number of genes, approaching 10% of human genes [82]. Similarly, the human homolog of Drosophila CG9772, called SCF^Skp2, is an ubiquitin ligase encoding onco-protein that is one of the components of the machinery used to control Myc levels through the ubiquitin pathway. Interestingly, this ubiquitin ligase is also required for induction of Myc-responsive genes, suggesting that ubiquitination not only promotes Myc turnover but also "licenses" its transcriptional activity [83]. Examples of homologous genes involved in human brain tumor formation are the tailless ortholog Nr2f1 involved in retinoblastoma, and the homologs of Myb (MYB), dm (MYC), and Notch (NOTCH2), whose misregulation has been associated with human gliomas [71, 84].

A core transcriptional signature of neoplastic brain tissue in Drosophila

In order to gain insights into the transcriptional alterations associated with brain tumors in vivo, we have determined the transcriptional signature of an adult brain tumor in Drosophila caused by homozygous mutation in the tumor suppressor gene brat. To this end, we have conducted two independent genome-wide gene expression studies using two different oligonucleotide microarray platforms to compare adult wildtype flies with flies displaying the adult brat^k06028mutant brain tumor. Cross-validation of our comparative analysis of adult wildtype heads compared to brat^k06028heads versus dissected brat^k06028brains compared to dissected adult brains of flies generated by transposon excision of the brat^k06028P-element identified a core transcriptional signature of brat^k06028neoplastic tissue. The obtained transcriptional signature revealed genetic alterations associated with malignant growth in a highly reproducible manner, despite the fact that both the biological material and the array platforms used differed considerably between the two experiments. This involved the use of stringent filter criteria such as a significance value of p ≤ 0.0001 (t-test). Due to these very stringent filter criteria, the 321 transcripts detected as differentially regulated in the brat^k06028tumor condition probably represent only a subset of genes that could be involved in brain tumorigenesis. For example, it has been shown that PTEN plays a crucial role in the highly conserved insulin/PI3-kinase (PI3K) signaling pathway involved in cell growth and proliferation (for review see [66] and the human homolog of PTEN is frequently lost in gliomas (reviewed in [84]). Although we did detect PTEN as differentially down-regulated in the brat tumor condition in both experiments, it was not incorporated into the core set of genes differentially regulated in brat^k06028neoplasms due to p-values below our cutoff value (see supplementary material).

Nevertheless, our genome-wide microarray analysis in Drosophila identified characteristic features of malignant growth such as those found in human cancer. Thus, we detected gene expression changes indicative of elevated and aberrant metabolic activity, elevated and aberrant activity of the cell cycle and basal transcriptional machinery, and elevated and aberrant activity of ribosome synthesis and translation control associated with the adult brat^k06028mutant brain tumor. Our data are in accordance with assumptions that the vast catalog of cancer cell genotypes is a manifestation of essential alterations in cell physiology that collectively dictate malignant growth: self-sufficiency in growth signals, insensitivity to growth-inhibitory (antigrowth) signals, evasion of programmed cell death (apoptosis), limitless replicative potential, and tissue invasion and metastasis [1]. This is further exemplified by the fact that for 20% of the genes differentially regulated in the Drosophila brat^k06028tumor condition it was possible to identify mammalian homologs associated with various types of tumors. These include homologs of Myb, Myc, cyclins as well as other genes involved in cell growth and proliferation control whose misregulation is frequently associated with mammalian cancer formation and malignant growth [1].

Brat acts as a tumor suppressor during Drosophila brain development

Brat was originally identified as a growth suppressor of the larval brain of Drosophila, with mutant brains growing up to eightfold larger than normal [28, 32]. Loss of brat also conferred metastatic potential upon tissues transplanted into wildtype host flies [34]. Moreover, our analysis of the brat^k06028allele shows that a substantial amount (15%) of homozygous mutant flies eclose and show a strong neoplastic adult brain phenotype with 100% penetrance. These data demonstrate that Brat functions as a tumor suppressor during brain development of Drosophila.

In general, tumor suppressors act as molecular antagonists to the formation of malignancy, a physiological state that is characterized by indefinite cellular growth and division [1]. This in turn requires elevated levels of transcriptional activity of genes involved in all aspects of aberrant cell cycle activity. Indeed, our genome-wide microarray expression analysis of brat^k06028neoplastic tissue revealed genetic alterations that can be attributed to malignant growth. We identified essential components of the cell cycle machinery like cyclin A, B, E, string and Myb as significantly upregulated in the tumor condition. This elevated cell cycle activity is accompanied by significant upregulation of genes involved in DNA replication, chromatin structure and cytokinesis, such as gluon, Mcm2 and Mcm3, and the Drosophila Myc homolog diminutive. Gain of function of these genes is frequently associated with malignant overgrowth and cancer in Drosophila [16]. Moreover, deregulation of the proto-oncogenes Myb and Myc characterizes one of the hallmarks of cancer, as the oncogenic activity of these genes in turn impairs differentiation and promotes growth [66]. In addition, we also observed significant upregulation of the Ran GTPase pathway since three members of it, namely Sd, ran, and msk, show elevated expression levels in the brat^k06028tumor condition. Recent data suggest that Ran GTPase signalling serves as a positional marker for the eukaryotic genome throughout the cell cycle by regulating microtubule nucleation and nuclear envelope formation around chromatin (reviewed in [85]). This together with the elevated and aberrant activity of several other structural components required for proliferation and cytokinesis emphasises the self-sufficiency of the cellular machinery during unrestrained growth.

Proliferation control via ribosome synthesis regulation by brat?

Recent results demonstrate that brat functions to repress ribosomal RNA synthesis and cell growth. brat mutant cells are larger than control cells, have enlarged nucleoli and contain excess rRNA. Furthermore, brat overexpression inhibits clone and organ growth, and leads to a decreased level of rRNA per cell [30]. Based on these observations it has been suggested, that in contrast to the cell growth regulation pathway of the insulin receptor and its effectors, brat affects cell growth not through the activity of ribosomes, but rather through the regulation of their synthesis [30]. Indeed, disruption of one or more of the steps that control protein synthesis has been associated with alterations in the cell cycle and regulation of cell growth. For example mutations that inactivate the tumor-suppressor activity of RB or p53, or both, result in aberrant upregulation of essential components in the protein synthesis machinery and increase ribosome biogenesis, leading to enhanced mRNA translation rates (reviewed in [86]). Among the 321 genes differentially regulated due to loss of the tumor suppressor activity of brat, we identified several genes involved in ribosome biogenesis and RNA processing as significantly upregulated in neoplastic tissue. One of the most prominent examples is ribosomal protein RpL3 which maps to the chromosomal region harbouring the Minute(3) 86D mutation [48]. However, although we observe several genes involved in protein synthesis and ribosome biogenesis as differentially upregulated in the brat^k06028tumor condition, it remains to be determined whether these phenomena represent the cause or consequence of tumor formation.

Neuroblast polarity and asymmetric division

Interestingly, we find that inscuteable (insc), staufen (stau), snail (sna), miranda (mira), and worniu (wor) have elevated transcript levels in brat^k06028brain tissue. Both insc and stau are involved in asymmetric segregation of cell fate determinants in dividing neuroblasts and insc expression is thought to be regulated by sna [52, 61]. Miranda is an adaptor protein involved in the distrinution of the Prospero gene product. Prospero is required to distinguish the proliferative capabilities of NBs and GMCs [67]. As it is known that the Discs large (Dlg), Scribble (Scrib) and Lethal giant larvae (Lgl) tumor suppressor proteins regulate multiple aspects of neuroblast asymmetric cell division these findings strengthen the link between tumorigenesis and asymmetry in cell division [87].

Taken together, our genome-wide microarray expression analysis of an adult brain tumor in Drosophila caused by mutation in the tumor suppressor gene brat identifies a core set of genes whose deregulation can be attributed to the genetic alterations associated with tumor formation and malignant growth. We detect aberrant and elevated activity in the transcription of genes involved in metabolic activity, cell cycle regulation, cytokinesis, as well as in basal transcriptional regulation. 20% of the genes we identified as differentially regulated in the brat^k06028tumor condition have mammalian homologs whose misregulation has been associated with various types of tumors. Moreover, we observed significant alterations in transcriptional activity of genes involved in ribosome biogenesis, translation and RNA processing due to loss of the translational repressor encoded by brat. This may indicate that loss of brat leads to overexpression of ribosomal proteins and increased ribosome function, which in turn may cause malignant growth. Considering the evolutionary conservation of gene structure and function, we propose that our results obtained by the first genome-wide expression analysis of an adult brain tumor in Drosophila will also be valid for tumor formation in mammals. It will now be important to elucidate at the cellular and molecular level the mechanisms involved in the genetic alterations caused by mutation of the tumor suppressor brat, which ultimately lead to tumor formation, metastasis and invasiveness.

Fly strains and genetics

The wildtype was Oregon R. The brat allele used in this study is brat^k06028, balanced over CyO, P{w^+mc = ActGFP}JMR1. Flies used for histology and for transcript profiling were 1 to 2 days old homozygous P{w^+mc = lacW}, brat^k06028males and females, identified by the absence of GFP. For the control condition of the oligonucleotide array experiment, either wildtype Oregon R were used or flies generated by transposon excision of the k brat^k06028P-element. Excision of the brat^k06028transposon resulted in several viable strains, which were further characterized by genomic PCR. Using the primer pair 5'AACAACCAAAACAACGGCAACC 3' and 5' AAACGGAGATAAGCCGACTTAC 3', flanking the insertion site of the brat^k06028transposon in the brat genomic DNA, a fragment of 211 bp was PCR amplified and sequenced from both strands. For all characterized strains, this sequence was indistinguishable from the genomic brat sequence of Oregon R control flies; strain "brat^k06028 jumpout" was chosen for further work. For assessing the percentage of adult escapers in brat^k06028, homozygous third instar larvae were identified by the absence of GFP, counted (n = 979) and collected into bottles containing standard food. Bottles were then subsequently checked on a daily basis for hatched flies during the following 4 weeks. All flies were kept on standard cornmeal/ yeast/ agar medium at 25°C.

Histology

For sectioning, adult flies were fixed, dehydrated, embedded in paraffin, and cut into 7 μm sections as described by Heisenberg and Böhl [88]. Sections were mounted on coated glass slides with DePeX (Fluka), and neural structures were visualized with a Zeiss Axioskop microscope by autofluorescence (wavelength 488 nm) and recorded using a Prog/Res/3008 digital camera (Kontron, Zürich).

Experimental paradigms and oligonucleotide arrays

Two oligonucleotide array experiments were carried out. In the first experiment, a custom-designed Drosophila oligonucleotide array (roDROMEGAa, Affymetrix) was used [89]. This array contains 14,090 sequences representing 13,369 single transcripts encoding Drosophila proteins deposited in SWISS-PROT/TrEMBL databases (Celera Genome/BDGP Release no.1; [10] as well as prokaryotic and custom chosen control sequences. This array was used to compare the transcriptome from Oregon R control heads with heads of homozygous brat^k06028flies; 5 replicates were done for the wildtype condition and 6 were done for homozygous brat^k06028(experiment A). In the second experiment, a commercially available Affymetrix Drosophila full genome array was used (Drosophila Genome Array, Affymetrix) which contains probe sets interrogating more than 13,500 genes from Drosophila melanogaster (experiment B). Sequences used to design this array can be found in Flybase (version 1). Over 8,000 of the genes represented have at least one EST/cDNA match. Additionally, approximately 5,500 genes were identified using prediction algorithms. This array was used to compare the transcriptome of dissected Drosophila brains derived from brat^k06028 jumpout with that of homozygous brat^k06028brains. 6 replicates were done for each of the two conditions.

Preparation of biotinylated cRNA

Total RNA was isolated from 100 fly heads or 150 fly brains, respectively, using guanidinium isothiocyanate in combination with acidic phenol (pH 4.3) (fast RNA tube green kit from BIO101) in a fast prep homogenizer FP120 (Bio 101). After precipitation, RNA was dissolved in DEPC-treated water (Ambion) and spectrophotometrically quantified using a GeneQuant RNA/DNA calculator (Pharmacia Biotech). The quality of the RNA was checked on a 0.5 × TBE agarose gel and the samples were stored at -80°C. cDNA was synthesized upon total RNA as a template, using the SuperScript Choice System for cDNA synthesis (Gibco/BRL) with a T7-(T)24 DNA primer. This primer (5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGG-(T)24VN-3') was PAGE-gel purified. For first strand cDNA synthesis, a typical 40 μl reaction contained 10–15 μg RNA, 200 pmoles T7-(T)24 primer, 500 μM of each dNTPs and 800 units reverse transcriptase (AMV Superscript II). The reaction was incubated for one hour at 42°C. Second strand cDNA synthesis was carried out at 18°C for two hours in a total volume of 340 μl, using 20 units E. coli DNA ligase, 80 units E. coli DNA polymerase I and 4 units RNase H in the presence of 250 μM of each dNTP. After 2nd strand cDNA synthesis, 0.5 μl RNase A (100 mg/ml) (Qiagen) was added and the samples were incubated at 37°C for half an hour. Thereafter 7.5 μl proteinase K (10 mg/ml) (Sigma) was added and the samples were further incubated at 37°C another half hour. After cDNA synthesis was completed, samples were phenol-chloroform extracted (3 times) using Phase Lock Gel (5 Prime-3 Prime, Inc.) and precipitated overnight at -20°C with 2.5 volumes 100 % ethanol. After percipitation, the samples were stored at -20°C. Biotinylated antisense cRNA was synthesized from the dsDNA template, using T7 RNA polymerase (MEGAscript T7 Kit, Ambion, Inc.).

A 20 μl reaction volume contained between 0.3–1.5 μg cDNA, 7.5 mM of both ATP and GTP, 5.6 mM of both UTP and CTP and 1.8 mM of both biotinylated Bio-16-UTP and Bio-11-CTP (ENZO diagnostics) and 2 μl 10 × T7 enzyme mix. The reaction was incubated at 37°C for 8 hours. Thereafter the unincorporated NTPs were removed by putting the sample over an RNeasy spin column (Qiagen). Aliquots of the reaction before and after cRNA synthesis were analyzed by agarose gel electrophoresis. Samples were precipitated overnight at -20°C, taken up in 20 μl DEPC treated water and spectrophotometrically quantified. Susequently, the biotinylated antisense cRNA was fragmented by heating the sample to 95°C for 35 min in a volume of 25 μl, containing 40 mM tris-acetate (pH 8.1), 100 mM KOAc, 30 mM MgOAc. After the fragmentation, the samples were placed on ice.

Hybridization and Scanning

Gene chips were equilibrated at RT and prehybridized with 280 μl hybridization buffer (1 × MES (pH 6.7)/ 1 M NaCl/ 20 mM EDTA/ 0.01 % Tween 20) for 15 min at 45°C with rotation at 60 rpm. Hybridization was done for 16 h at 45°C (60 rpm) in a final volume of 220 μl hybridization buffer, 0.1 mg/ml herring sperm DNA, 0.5 mg/ml acetylated BSA containing 50 pM control oligo B2, spiked bacterial control RNAs (BioB, BioC, BioD, cre) and 15 μg fragmented biotinylated cRNA. After hybridization, the arrays were briefly rinsed with 6 × SSPE-T buffer (0.9 M NaCl/ 0.06 M NaH₂PO₄/ 6 mM EDTA/ 0.01 %Tween 20). Washing and staining were carried out on a Fluidics station 400 (Affymetrix, MAS 5.0) using 100 mM MES, 0.1 M NaCl, 0,01 % Tween 20 as stringent wash buffer. Two different staining protocols were used. When comparing the transcriptome of Oregon R fly heads and homozygous brat^k06028fly heads, staining was carried out in 220 μl 1 × MES buffer, 2.0 mg/ml acetylated BSA, 10 ng/ ml StreptavidinR-phycoerythrin conjugate (Molecular Probes) at 40°C for 15 min. In the second experiment comparing dissected brains of brat^k06028 jumpout flies with homozygous brat^k06028brains, an additional antibody amplification step was included. Prior to staining, 220 μl 1 × MES buffer containing 2 mg/ml acetylated BSA and 2 μg/ml biotinylated anti-streptavidin antibody (Vector) was applied onto the array for 30 min at 40°C. Following a washing step, staining was then performed using streptavidinR-phycoerythrin as described above. Following a final washing step, arrays were scanned with a commercial confocal laser scanner (Agilent).

Data analysis

The single arrays were analysed using Microarray Suite 5.0 (Affymetrix). Detailed analysis was carried out using Race-A (Roche), Excel2000 and Access2000 software (Microsoft), and GeneSpring (Silicon Genetics). In this process, the mean of all signal values of each array were set to 1, and all other signal values were expressed relative to it (normalization). A Nalimov outlier test was performed [90]. The data was grouped into control- and experimental conditions. For every gene in all 4 conditions a single mean signal value over the replicates was calculated. Also an unpaired t-test was performed for every gene per experiment to assess significance of change between conditions. The annotation of the two different arrays was linked (see below) using information from NetAffx [91] and Flybase [92]. Genes were judged as present when their presence attribute in the RACE-A software was above 0.75. Both datasets were filtered using a cutoff on the mean signal (≥ 10 in the condition the gene had to be present in), the changefold (≥ 2) and the t-test (≤ 0.0001). The various control sequences present on the arrays were excluded from analysis. Subsequently the overlap of the two datasets was calculated using the database mentioned above. When more than one probe set was representing a gene the probe set with the highest statistical significance as judged by the t-test was chosen. In those cases were the p-value was identical, the probe set with the higher expression levels was chosen. The raw data of the two experiments can be found in Supplementary Table 1 and Supplementary Table 2 [see Additional files 1 and 2].

Linking the two different oligonucleotide arrays

Using the probe set annotation information provided for the roDROMEGa array, we determined which probe sets on the Affymetrix "DrosGenome1" GeneChip were analogous to those on the Roche custom chip. To this end, a database was built with the following information: 1) Affymetrix "DrosGenome1" probe set annotations, 2) roDROMEGa annotations, and 3) FlyBase gene annotations.

The Affymetrix "DrosGenome1" annotations were extracted from the June 2003 DrosGenome1_annot_csv.zip file obtained from the NetAffx [91] Affymetrix web site http://www.affymetrix.com/Auth/analysis/downloads/taf/DrosGenome1_annot_csv.zip. The database tables contained probe set ID, gene symbol, FlyBase FBgns, and synonyms as annotated by Affymetrix for each of the 14,010 probe sets on DrosGenome1. The roDROMEGa array annotations were extracted from text dump files obtained from Roche. The database table contained probe set ID and gene symbol information for each of the 14,090 probe sets on the Roche custom chip. FlyBase gene annotations were extracted from the 29-08-2003 FlyBase genes dump file located at ftp://flybase.org/flybase/genes/genes.txt. The resulting database tables contained gene symbol, FBgns, and synonyms for all 43,177 annotations in FlyBase. During extraction of Affymetrix DrosGenome1 and FlyBase annotations from their source files, special cleanup was required to correct from HTML escape characters commonly used in Drosophila melanogaster annotation dump files. This was necessary because the annotations for the roDROMEGa array contained only fully spelled out words for gene symbols (i.e. 'alphaCop' instead of '&agr;Cop').

The 14,090 probe sets on the Roche custom Drosophila microarray represented 13,343 distinct genes which were mapped to Affymetrix's DrosGenome1. Utilizing the database described above, we first mapped gene symbols annotated by Roche on their custom chip to gene symbols and synonyms annotated by Affymetrix on their DrosGenome1 chip. This search yielded 13,377 probe sets (12,784 genes) on the Roche Drosophila custom chip. The remaining 713 probe sets on the Roche custom chip which had no direct gene symbol or synonym match to DrosGenome1 were used in a three-way search using FlyBase genes data. These 713 Roche gene symbols were searched against all FlyBase gene symbols and synonyms and in case of a match all the corresponding FlyBase FBgns for these genes were used to search the Affymetrix FBgn annotations on DrosGenome1. This search resulted in additional 449 probe sets (401 genes) on the Roche Drosophila custom chip which mapped to analogous probe sets on DrosGenome1. Of the remaining 172 unmapped Roche probe sets, 92 were found to be either Affymetrix AFFX controls, Roche internal controls, or rRNA controls. Thus, 13,826 Roche probe sets (13,185 genes) have at least one analogous Affymetrix DrosGenome1 probe set. 172 probe sets (157 genes) on the Roche Drosophila custom chip dolfed our parameters and remained orphaned with no Affymetrix DrosGenome1 match. A list with the orphaned probe sets can be found in Supplementay Table 4 [see Additional file 4]

Functional classification and links to mammalian homologs involved in cancer

The 321 genes differentially regulated in both experiments were grouped into functional classes by reviewing their annotation in Flybase. When there were vertebrate homologs annotated in Flybase we searched the PubMed bibliographic database http://www.ncbi.nlm.nih.gov/PubMed/medline.html for literature that linked the gene in question to cancer in vertebrates.

We thank Martin Müller and Bruno Bello for comments on the manuscript. We also thank Clemens Broger, Martin Strahm, and Martin Neeb (F. Hoffman-La Roche) for support with the RACE-A software. We are grateful to the students of the Neurogenomics course 2002 for actively participating in this project, to Philippe Demougin for help with the array hybridizations and to Margit Böhler for assistance with the histology. This research was supported by grants from the SNSF and Neurex/Eltem

Affiliations

1. Institute of Zoology, Biocenter/Pharmacenter, University of Basel, Klingelbergstr. 50, CH-4056, Basel, Switzerland

   • Thomas Loop
   • , Ronny Leemans
   • , Urs Stiefel
   • , Boris Egger
   • , Fukang Xie
   • , Heinrich Reichert
   •  & Frank Hirth

2. Biocenter, University of Basel, Klingelbergstr. 70, CH-4056, Basel, Switzerland

   • Leandro Hermida
   •  & Michael Primig

3. Roche Genetics Pharmaceuticals Division, F. Hoffmann-La Roche Ltd, CH-4070, Basel, Switzerland

   • Ulrich Certa

4. Institute of Biology III, University of Freiburg, Schaenzlestrasse 1, D-79104, Freiburg, Germany

   • Karl-Friedrich Fischbach

Corresponding author

Correspondence to Frank Hirth.

Authors' contributions

TL carried out the the array study of experiment B, performed the histology of flyheads, did the bioinformatics and statistics for both chip experiments and drafted the manuscript. RL carried out the array study of experiment A and performed the jumpout of the brat^k06028P-element. US dissected the flies that were used as a tissue source for the RNA extractions. LH carried out the bioinformatics for linking the annotations of the two different array platforms used in our study. BE balanced the fly strains. FX was initially involved in setting up experimental conditions and carried out pilot experiments. MP, UC, KFF and HR participated in the design and coordination of the study. FH conceived of the study, and participated in its design and coordination. All authors read and approved the final manuscript.

Reprints and Permissions