The ftz upstream element drives late ftz stripes but is not required for regulation of Ftz target genes.

The regulation of gene expression in precise, rapidly changing spatial patterns is essential for embryonic development. Multiple enhancers have been identified for the evolving expression patterns of the cascade of Drosophila segmentation genes that establish the basic body plan of the fly. Classic reporter transgene experiments identified multiple cis-regulatory elements (CREs) that are sufficient to direct various aspects of the evolving expression pattern of the pair-rule gene fushi tarazu (ftz). These include enhancers that coordinately activate expression in all seven stripes and stripe-specific elements that activate expression in one or more ftz stripes. Of the two 7-stripe enhancers, analysis of reporter transgenes demonstrated that the upstream element (UPS) is autoregulatory, requiring direct binding of Ftz protein to direct striped expression. Here, we asked about the endogenous role of the UPS by precisely deleting this 7-stripe enhancer. In ftzΔUPS7S homozygotes, ftz stripes appear in the same order as wildtype, and all but stripe 4 are expressed at wildtype levels by the end of the cellular blastoderm stage. This suggests that the zebra element and UPS harbor information to direct stripe 4 expression, although previous deletion analyses failed to identify a stripe-specific CRE within these two 7-stripe enhancers. However, the UPS is necessary for late ftz stripe expression, with all 7 stripes decaying earlier than wildtype in ftzΔUPS7S homozygotes. Despite this premature loss of ftz expression, downstream target gene regulation proceeds as in wildtype, and segmentation is unperturbed in the overwhelming majority of animals. We propose that this late-acting enhancer provides a buffer against perturbations in gene expression but is not required for establishment of Ftz cell fates. Overall, our results demonstrate that multiple enhancers, each directing distinct aspects of an overall gene expression pattern, contribute to fine-tuning the complex patterns necessary for embryonic development.

A bistable autoregulatory module in the developing embryo commits cells to binary expression fates.

Bistable autoactivation has been proposed as a mechanism for cells to adopt binary fates during embryonic development. However, it is unclear whether the autoactivating modules found within developmental gene regulatory networks are bistable, unless their parameters are quantitatively determined. Here, we combine in vivo live imaging with mathematical modeling to dissect the binary cell fate dynamics of the fruit fly pair-rule gene fushi tarazu (ftz), which is regulated by two known enhancers: the early (non-autoregulating) element and the autoregulatory element. Live imaging of transcription and protein concentration in the blastoderm revealed that binary Ftz fates are achieved as Ftz expression rapidly transitions from being dictated by the early element to the autoregulatory element. Moreover, we discovered that Ftz concentration alone is insufficient to activate the autoregulatory element, and that this element only becomes responsive to Ftz at a prescribed developmental time. Based on these observations, we developed a dynamical systems model and quantitated its kinetic parameters directly from experimental measurements. Our model demonstrated that the ftz autoregulatory module is indeed bistable and that the early element transiently establishes the content of the binary cell fate decision to which the autoregulatory module then commits. Further in silico analysis revealed that the autoregulatory element locks the Ftz fate quickly, within 35 min of exposure to the transient signal of the early element. Overall, our work confirms the widely held hypothesis that autoregulation can establish developmental fates through bistability and, most importantly, provides a framework for the quantitative dissection of cellular decision-making.

Two pair-rule responsive enhancers regulate wingless transcription in the Drosophila blastoderm embryo.

BACKGROUND: While many developmentally relevant enhancers act in a modular fashion, there is growing evidence for nonadditive interactions between distinct cis-regulatory enhancers. We investigated if nonautonomous enhancer interactions underlie transcription regulation of the Drosophila segment polarity gene, wingless. RESULTS: We identified two wg enhancers active at the blastoderm stage: wg 3613u, located from -3.6 to -1.3 kb upstream of the wg transcription start site (TSS) and 3046d, located in intron two of the wg gene, from 3.0 to 4.6 kb downstream of the TSS. Genetic experiments confirm that Even Skipped (Eve), Fushi-tarazu (Ftz), Runt, Odd-paired (Opa), Odd-skipped (Odd), and Paired (Prd) contribute to spatially regulated wg expression. Interestingly, there are enhancer specific differences in response to the gain or loss of function of pair-rule gene activity. Although each element recapitulates aspects of wg expression, a composite reporter containing both enhancers more faithfully recapitulates wg regulation than would be predicted from the sum of their individual responses. CONCLUSION: These results suggest that the regulation of wg by pair-rule genes involves nonadditive interactions between distinct cis-regulatory enhancers.

Maternal Nanos inhibits Importin-α2/Pendulin-dependent nuclear import to prevent somatic gene expression in the Drosophila germline.

Repression of somatic gene expression in germline progenitors is one of the critical mechanisms involved in establishing the germ/soma dichotomy. In Drosophila, the maternal Nanos (Nos) and Polar granule component (Pgc) proteins are required for repression of somatic gene expression in the primordial germ cells, or pole cells. Pgc suppresses RNA polymerase II-dependent global transcription in pole cells, but it remains unclear how Nos represses somatic gene expression. Here, we show that Nos represses somatic gene expression by inhibiting translation of maternal importin-α2 (impα2) mRNA. Mis-expression of Impα2 caused aberrant nuclear import of a transcriptional activator, Ftz-F1, which in turn activated a somatic gene, fushi tarazu (ftz), in pole cells when Pgc-dependent transcriptional repression was impaired. Because ftz expression was not fully activated in pole cells in the absence of either Nos or Pgc, we propose that Nos-dependent repression of nuclear import of transcriptional activator(s) and Pgc-dependent suppression of global transcription act as a 'double-lock' mechanism to inhibit somatic gene expression in germline progenitors.

Identification of evolutionarily conserved downstream core promoter elements required for the transcriptional regulation of Fushi tarazu target genes.

The regulation of transcription initiation is critical for developmental and cellular processes. RNA polymerase II (Pol II) is recruited by the basal transcription machinery to the core promoter where Pol II initiates transcription. The core promoter encompasses the region from -40 to +40 bp relative to the +1 transcription start site (TSS). Core promoters may contain one or more core promoter motifs that confer specific properties to the core promoter, such as the TATA box, initiator (Inr) and motifs that are located downstream of the TSS, namely, motif 10 element (MTE), the downstream core promoter element (DPE) and the Bridge, a bipartite core promoter element. We had previously shown that Caudal, an enhancer-binding homeodomain transcription factor and a key regulator of the Hox gene network, is a DPE-specific activator. Interestingly, pair-rule proteins have been implicated in enhancer-promoter communication at the engrailed locus. Fushi tarazu (Ftz) is an enhancer-binding homeodomain transcription factor encoded by the ftz pair-rule gene. Ftz works in concert with its co-factor, Ftz-F1, to activate transcription. Here, we examined whether Ftz and Ftz-F1 activate transcription with a preference for a specific core promoter motif. Our analysis revealed that similarly to Caudal, Ftz and Ftz-F1 activate the promoter containing a TATA box mutation to significantly higher levels than the promoter containing a DPE mutation, thus demonstrating a preference for the DPE motif. We further discovered that Ftz target genes are enriched for a combination of functional downstream core promoter elements that are conserved among Drosophila species. Thus, the unique combination (Inr, Bridge and DPE) of functional downstream core promoter elements within Ftz target genes highlights the complexity of transcriptional regulation via the core promoter in the transcription of different developmental gene regulatory networks.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin.

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Timing of embryonic quiescence determines viability of embryos from the calanoid copepod, Acartia tonsa (Dana).

Like 41 other calanoid copepods, Acartia tonsa, are capable of inducing embryonic quiescence when experiencing unfavorable environmental conditions. The ecdysone-signaling cascade is known to have a key function in developmental processes like embryogenesis and molting of arthropods, including copepods. We examined the role of ecdysteroid-phosphate phosphatase (EPPase), ecdysone receptor (EcR), ß fushi tarazu transcription factor 1 (ßFTZ-F1), and the ecdysteroid-regulated early gene E74 (E74), which represent different levels of the ecdysone-signaling cascade in our calanoid model organism. Progression of embryogenesis was monitored and hatching success determined to evaluate viability. Embryos that were induced quiescence before the gastrulation stage would stay in gastrulation during the rest of quiescence and exhibited a slower pace of hatching as compared to subitaneous embryos. In contrast, embryos developed further than gastrulation would stay in gastrulation or later stages during quiescence and showed a rapid pace in hatching after quiescence termination. Expression patterns suggested two peaks of the biological active ecdysteroids, 20-hydroxyecdysone (20E). The first peak of 20E was expressed in concert with the beginning of embryogenesis originating from yolk-conjugated ecdysteroids, based on EPPase expression. The second peak is suggested to originate from de novo synthesized 20E around the limb bud stage. During quiescence, the expression patterns of EPPase, EcR, ßFTZ-F1, and E74 were either decreasing or not changing over time. This suggests that the ecdysone-signaling pathway play a key role in the subitaneous development of A. tonsa embryogenesis, but not during quiescence. The observation is of profound ecological and practical relevance for the dynamics of egg banks.

MicroRNA-34 directly targets pair-rule genes and cytoskeleton component in the honey bee.

MicroRNAs (miRNAs) are key regulators of developmental processes, such as cell fate determination and differentiation. Previous studies showed Dicer knockdown in honeybee embryos disrupt the processing of functional mature miRNAs and impairs embryo patterning. Here we investigated the expression profiles of miRNAs in honeybee embryogenesis and the role of the highly conserved miR-34-5p in the regulation of genes involved in insect segmentation. A total of 221 miRNAs were expressed in honey bee embryogenesis among which 97 mature miRNA sequences have not been observed before. Interestingly, we observed a switch in dominance between the 5-prime and 3-prime arm of some miRNAs in different embryonic stages; however, most miRNAs present one dominant arm across all stages of embryogenesis. Our genome-wide analysis of putative miRNA-target networks and functional pathways indicates miR-34-5p is one of the most conserved and connected miRNAs associated with the regulation of genes involved in embryonic patterning and development. In addition, we experimentally validated that miR-34-5p directly interacts to regulatory elements in the 3'-untranslated regions of pair-rule (even-skipped, hairy, fushi-tarazu transcription factor 1) and cytoskeleton (actin5C) genes. Our study suggests that miR-34-5p may regulate the expression of pair-rule and cytoskeleton genes during early development and control insect segmentation.

Different modes of enhancer-specific regulation by Runt and Even-skipped during Drosophila segmentation.

The initial metameric expression of the Drosophila sloppy paired 1 (slp1) gene is controlled by two distinct cis-regulatory DNA elements that interact in a nonadditive manner to integrate inputs from transcription factors encoded by the pair-rule segmentation genes. We performed chromatin immunoprecipitation on reporter genes containing these elements in different embryonic genotypes to investigate the mechanism of their regulation. The distal early stripe element (DESE) mediates both activation and repression by Runt. We find that the differential response of DESE to Runt is due to an inhibitory effect of Fushi tarazu (Ftz) on P-TEFb recruitment and the regulation of RNA polymerase II (Pol II) pausing. The proximal early stripe element (PESE) is also repressed by Runt, but in this case, Runt prevents PESE-dependent Pol II recruitment and preinitiation complex (PIC) assembly. PESE is also repressed by Even-skipped (Eve), but, of interest, this repression involves regulation of P-TEFb recruitment and promoter-proximal Pol II pausing. These results demonstrate that the mode of slp1 repression by Runt is enhancer specific, whereas the mode of repression of the slp1 PESE enhancer is transcription factor specific. We propose a model based on these differential regulatory interactions that accounts for the nonadditive interactions between the PESE and DESE enhancers during Drosophila segmentation.

Non-specificity of transcription factor function in Drosophila melanogaster.

A major problem in developmental genetics is how HOX transcription factors, like Proboscipedia (PB) and Ultrabithorax (UBX), regulate distinct programs of gene expression to result in a proboscis versus a haltere, respectively, when the DNA-binding homeodomain (HD) of HOX transcription factors recognizes similar DNA-binding sequences. Indeed, the lack of DNA-binding specificity is a problem for all transcription factors (TFs), as the DNA-binding domains generally recognize small targets of five to six bases in length. Although not the initial intent of the study, I found extensive non-specificity of TF function. Multiple TFs including HOX and HD-containing and non-HD-containing TFs induced both wingless and eyeless phenotypes. The TFs Labial (LAB), Deformed (DFD), Fushi tarazu (FTZ), and Squeeze (SQZ) induced ectopic larval thoracic (T) 1 beard formation in T2 and T3. The TF Doublesex male (DSXM) rescued the reduced maxillary palp pb phenotype. These examples of non-specificity of TF function across classes of TFs, combined with previous observations, compromise the implicit, initial assumption often made that an intrinsic mechanism of TF specificity is important for function. Interestingly, the functional complementation of the pb phenotype may suggest a larger role for regulation of expression of TFs in restriction of function as opposed to an intrinsic specificity of TF function. These observations have major ramifications for analysis of functional conservation in evolution and development.

Activation of Ftz-F1-Responsive Genes through Ftz/Ftz-F1 Dependent Enhancers.

The orphan nuclear receptor Ftz-F1 is expressed in all somatic nuclei in Drosophila embryos, but mutations result in a pair-rule phenotype. This was explained by the interaction of Ftz-F1 with the homeodomain protein Ftz that is expressed in stripes in the primordia of segments missing in either ftz-f1 or ftz mutants. Ftz-F1 and Ftz were shown to physically interact and coordinately activate the expression of ftz itself and engrailed by synergistic binding to composite Ftz-F1/Ftz binding sites. However, attempts to identify additional target genes on the basis of Ftz-F1/ Ftz binding alone has met with only limited success. To discern rules for Ftz-F1 target site selection in vivo and to identify additional target genes, a microarray analysis was performed comparing wildtype and ftz-f1 mutant embryos. Ftz-F1-responsive genes most highly regulated included engrailed and nine additional genes expressed in patterns dependent on both ftz and ftz-f1. Candidate enhancers for these genes were identified by combining BDTNP Ftz ChIP-chip data with a computational search for Ftz-F1 binding sites. Of eight enhancer reporter genes tested in transgenic embryos, six generated expression patterns similar to the corresponding endogenous gene and expression was lost in ftz mutants. These studies identified a new set of Ftz-F1 targets, all of which are co-regulated by Ftz. Comparative analysis of enhancers containing Ftz/Ftz-F1 binding sites that were or were not bona fide targets in vivo suggested that GAF negatively regulates enhancers that contain Ftz/Ftz-F1 binding sites but are not actually utilized. These targets include other regulatory factors as well as genes involved directly in morphogenesis, providing insight into how pair-rule genes establish the body pattern.

Identification of MBF2 family genes in Bombyx mori and their expression in different tissues and stages and in response to Bacillus bombysepticus infection and starvation.

The Multiprotein bridge factor 2 (MBF2) gene was first identified as a co-activator involved in BmFTZ-F1-mediated activation of the Fushi tarazu gene. Herein, nine homologous genes of MBF2 gene are identified. Evolutionary analysis showed that this gene family is insect-specific and that the family members are closely related to response to pathogens (REPAT) genes. Tissue distribution analysis revealed that these genes could be expressed in a tissue-specific manner. Developmental profiles analysis showed that the MBF2 gene family members were highly expressed in the different stages. Analysis of the expression patterns of nine MBF2 family genes showed that Bacillus bombysepticus treatment induced the up-regulation of several MBF2 family genes, including MBF2-4, -7, -9, -8. Furthermore, we found the MBF2 family genes were modulated by starvation and the expression of these genes recovered upon re-feeding, except for MBF2-5, -9. These findings suggested roles for these proteins in insect defense against pathogens and nutrient metabolism, which has an important guiding significance for designing pest control strategies.

Acute phenanthrene toxicity to juvenile diploid and triploid African catfish (Clarias gariepinus): Molecular, biochemical, and histopathological alterations.

Information on the biological responses of polyploid animals towards environmental contaminants is scarce. This study aimed to compare reproductive axis-related gene expressions in the brain, plasma biochemical responses, and the liver and gill histopathological alterations in diploid and triploid full-sibling juvenile African catfish (Clarias gariepinus). Fish were exposed for 96 h to one of the two waterborne phenanthrene (Phe) concentrations [mean measured (SD): 6.2 (2.4) and 76 (4.2) µg/L]. In triploids, exposure to 76 µg/L Phe increased mRNA level of fushi tarazu-factor 1 (ftz-f1). Expression of tryptophan hydroxylase2 (tph2) was also elevated in both ploidies following the exposure to 76 µg/L Phe compared to the solvent control. In triploids, 76 µg/L Phe increased plasma alkaline phosphatase (ALP) and lactate dehydrogenase (LDH) levels compared to the other Phe-exposed group. It also elevated lactate and glucose contents relative to the other groups. In diploids, however, biochemical biomarkers did not change. Phenanthrene exposures elevated glycogen contents and the prevalence of histopathological lesions in the liver and gills of both ploidies. This study showed substantial differences between diploids and triploids on biochemical and molecular biomarker responses, but similar histopathological alterations following acute Phe exposures.

Quantitative dynamics and increased variability of segmentation gene expression in the Drosophila Krüppel and knirps mutants.

Here we characterize the response of the Drosophila segmentation system to mutations in two gap genes, Kr and kni, in the form of single or double homozygotes and single heterozygotes. Segmentation gene expression in these genotypes was quantitatively monitored with cellular resolution in space and 6.5 to 13min resolution in time. As is the case with wild type, we found that gene expression domains in the posterior portion of the embryo shift to the anterior over time. In certain cases, such as the gt posterior domain in Kr mutants, the shifts are significantly larger than is seen in wild type embryos. We also investigated the effects of Kr and kni on the variability of gene expression. Mutations often produce variable phenotypes, and it is well known that the cuticular phenotype of Kr mutants is variable. We sought to understand the molecular basis of this effect. We find that throughout cycle 14A the relative levels of eve and ftz expression in stripes 2 and 3 are variable among individual embryos. Moreover, in Kr and kni mutants, unlike wild type, the variability in positioning of the posterior Hb domain and eve stripe 7 is not decreased or filtered with time. The posterior Gt domain in Kr mutants is highly variable at early times, but this variability decreases when this domain shifts in the anterior direction to the position of the neighboring Kni domain. In contrast to these findings, positional variability throughout the embryo does not decrease over time in double Kr;kni mutants. In heterozygotes the early expression patterns of segmentation genes resemble patterns seen in homozygous mutants but by the onset of gastrulation they become similar to the wild type patterns. Finally, we note that gene expression levels are reduced in Kr and kni mutant embryos and have a tendency to decrease over time. This is a surprising result in view of the role that mutual repression is thought to play in the gap gene system.

Ecdysone-responsive transcription factors determine the expression region of target cuticular protein genes in the epidermis of Bombyx mori.

In the present study, we found that different ecdysone-responsive transcription factors were expressed differentially in different regions of the epidermis at around pupation. ßFTZ-F1 transcripts were strongly but E74A transcripts were barely observed in the thoracic region of the epidermis, and vice versa in the abdominal region. Transcripts of all the examined transcription factors were observed in wing disc. Transcript of a cuticular protein gene, BMWCP4, which does not have a ßFTZ-F1 binding site in the 2-kb upstream region, was not observed in the thoracic region of the epidermis. Transcript of BMWCP9, which does not have an E74 binding site in the 2-kb upstream region, was not observed in the abdominal region of the epidermis. BMWCP2 has all the transcription factor binding sites examined and was expressed in the thoracic and abdominal region of the epidermis. Thus, it is suggested that ecdysone-responsive transcription factors determined the space where the cuticular protein genes were expressed, which, in turn, determined the character of the cuticle that was characterized by the combination of cuticular proteins.

Huckebein is part of a combinatorial repression code in the anterior blastoderm.

The hierarchy of the segmentation cascade responsible for establishing the Drosophila body plan is composed by gap, pair-rule and segment polarity genes. However, no pair-rule stripes are formed in the anterior regions of the embryo. This lack of stripe formation, as well as other evidence from the literature that is further investigated here, led us to the hypothesis that anterior gap genes might be involved in a combinatorial mechanism responsible for repressing the cis-regulatory modules (CRMs) of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) anterior-most stripes. In this study, we investigated huckebein (hkb), which has a gap expression domain at the anterior tip of the embryo. Using genetic methods we were able to detect deviations from the wild-type patterns of the anterior-most pair-rule stripes in different genetic backgrounds, which were consistent with Hkb-mediated repression. Moreover, we developed an image processing tool that, for the most part, confirmed our assumptions. Using an hkb misexpression system, we further detected specific repression on anterior stripes. Furthermore, bioinformatics analysis predicted an increased significance of binding site clusters in the CRMs of h 1, eve 1, run 1 and ftz 1when Hkb was incorporated in the analysis, indicating that Hkb plays a direct role in these CRMs. We further discuss that Hkb and Slp1, which is the other previously identified common repressor of anterior stripes, might participate in a combinatorial repression mechanism controlling stripe CRMs in the anterior parts of the embryo and define the borders of these anterior stripes.

How to make stripes: deciphering the transition from non-periodic to periodic patterns in Drosophila segmentation.

The generation of metameric body plans is a key process in development. In Drosophila segmentation, periodicity is established rapidly through the complex transcriptional regulation of the pair-rule genes. The 'primary' pair-rule genes generate their 7-stripe expression through stripe-specific cis-regulatory elements controlled by the preceding non-periodic maternal and gap gene patterns, whereas 'secondary' pair-rule genes are thought to rely on 7-stripe elements that read off the already periodic primary pair-rule patterns. Using a combination of computational and experimental approaches, we have conducted a comprehensive systems-level examination of the regulatory architecture underlying pair-rule stripe formation. We find that runt (run), fushi tarazu (ftz) and odd skipped (odd) establish most of their pattern through stripe-specific elements, arguing for a reclassification of ftz and odd as primary pair-rule genes. In the case of run, we observe long-range cis-regulation across multiple intervening genes. The 7-stripe elements of run, ftz and odd are active concurrently with the stripe-specific elements, indicating that maternal/gap-mediated control and pair-rule gene cross-regulation are closely integrated. Stripe-specific elements fall into three distinct classes based on their principal repressive gap factor input; stripe positions along the gap gradients correlate with the strength of predicted input. The prevalence of cis-elements that generate two stripes and their genomic organization suggest that single-stripe elements arose by splitting and subfunctionalization of ancestral dual-stripe elements. Overall, our study provides a greatly improved understanding of how periodic patterns are established in the Drosophila embryo.

On the border of the homeotic function: re-evaluating the controversial role of cofactor-recruiting motifs: the role of cofactor-recruiting motifs in conferring Hox evolutionary flexibility may critically depend on the protein environment.

In this review we present concepts that challenge a recently emerging paradigm explaining how similar Hox proteins perform different developmental functions across evolution, despite relatively limited sequence variability. This paradigm relates to the transcription factor, Fushi tarazu (Ftz), whose evolutionary plasticity has been shown to rely on the shuffling between two short protein recognition motifs. We discuss the Ftz paradigm and consider alternative interpretations to the evolutionary flexibility of this Hox protein. In particular, we propose that the protein environment might have played a critical role in the functional shuffling of Ftz during arthropod evolution.

Visual exploration of three-dimensional gene expression using physical views and linked abstract views.

During animal development, complex patterns of gene expression provide positional information within the embryo. To better understand the underlying gene regulatory networks, the Berkeley Drosophila Transcription Network Project (BDTNP) has developed methods that support quantitative computational analysis of three-dimensional (3D) gene expression in early Drosophila embryos at cellular resolution. We introduce PointCloudXplore (PCX), an interactive visualization tool that supports visual exploration of relationships between different genes' expression using a combination of established visualization techniques. Two aspects of gene expression are of particular interest: 1) gene expression patterns defined by the spatial locations of cells expressing a gene and 2) relationships between the expression levels of multiple genes. PCX provides users with two corresponding classes of data views: 1) Physical Views based on the spatial relationships of cells in the embryo and 2) Abstract Views that discard spatial information and plot expression levels of multiple genes with respect to each other. Cell Selectors highlight data associated with subsets of embryo cells within a View. Using linking, these selected cells can be viewed in multiple representations. We describe PCX as a 3D gene expression visualization tool and provide examples of how it has been used by BDTNP biologists to generate new hypotheses.