The long non-coding RNA ET-20 mediates EMT by impairing desmosomes in breast cancer cells.

The vast majority of breast cancer-associated deaths are due to metastatic spread of cancer cells, a process aided by epithelial-to-mesenchymal transition (EMT). Mounting evidence has indicated that long non-coding RNAs (lncRNAs) also contribute to tumor progression. We report the identification of 114 novel lncRNAs that change their expression during TGFß-induced EMT in murine breast cancer cells (referred to as EMT-associated transcripts; ETs). Of these, the ET-20 gene localizes in antisense orientation within the tenascin C (Tnc) gene locus. TNC is an extracellular matrix protein that is critical for EMT and metastasis formation. Both ET-20 and Tnc are regulated by the EMT master transcription factor Sox4. Notably, ablation of ET-20 lncRNA effectively blocks Tnc expression and with it EMT. Mechanistically, ET-20 interacts with desmosomal proteins, thereby impairing epithelial desmosomes and promoting EMT. A short transcript variant of ET-20 is shown to be upregulated in invasive human breast cancer cell lines, where it also promotes EMT. Targeting ET-20 appears to be a therapeutically attractive lead to restrain EMT and breast cancer metastasis in addition to its potential utility as a biomarker for invasive breast cancer.

Homeodomain proteins: an update.

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.

Active transcriptomic and proteomic reprogramming in the C. elegans nucleotide excision repair mutant xpa-1.

Transcription-blocking oxidative DNA damage is believed to contribute to aging and to underlie activation of oxidative stress responses and down-regulation of insulin-like signaling (ILS) in Nucleotide Excision Repair (NER) deficient mice. Here, we present the first quantitative proteomic description of the Caenorhabditis elegans NER-defective xpa-1 mutant and compare the proteome and transcriptome signatures. Both methods indicated activation of oxidative stress responses, which was substantiated biochemically by a bioenergetic shift involving increased steady-state reactive oxygen species (ROS) and Adenosine triphosphate (ATP) levels. We identify the lesion-detection enzymes of Base Excision Repair (NTH-1) and global genome NER (XPC-1 and DDB-1) as upstream requirements for transcriptomic reprogramming as RNA-interference mediated depletion of these enzymes prevented up-regulation of genes over-expressed in the xpa-1 mutant. The transcription factors SKN-1 and SLR-2, but not DAF-16, were identified as effectors of reprogramming. As shown in human XPA cells, the levels of transcription-blocking 8,5'-cyclo-2'-deoxyadenosine lesions were reduced in the xpa-1 mutant compared to the wild type. Hence, accumulation of cyclopurines is unlikely to be sufficient for reprogramming. Instead, our data support a model where the lesion-detection enzymes NTH-1, XPC-1 and DDB-1 play active roles to generate a genomic stress signal sufficiently strong to result in transcriptomic reprogramming in the xpa-1 mutant.

The LIM homeobox gene ceh-14 is required for phasmid function and neurite outgrowth.

Transcription factors play key roles in cell fate specification and cell differentiation. Previously, we showed that the LIM homeodomain factor CEH-14 is expressed in the AFD neurons where it is required for thermotaxis behavior in Caenorhabditis elegans. Here, we show that ceh-14 is expressed in the phasmid sensory neurons, PHA and PHB, a number of neurons in the tail, i.e., PHC, DVC, PVC, PVN, PVQ, PVT, PVW and PVR, as well as the touch neurons. Analysis of the promoter region shows that important regulatory elements for the expression in most neurons reside from -4kb to -1.65kb upstream of the start codon. Further, within the first introns are elements for expression in the hypodermis. Phylogenetic footprinting revealed numerous conserved motifs in these regions. In addition to the existing deletion mutation ceh-14(ch3), we isolated a new allele, ceh-14(ch2), in which only one LIM domain is disrupted. The latter mutant allele is partially defective for thermosensation. Analysis of both mutant alleles showed that they are defective in phasmid dye-filling. However, the cell body, dendritic outgrowth and ciliated endings of PHA and PHB appear normal, indicating that ceh-14 is not required for growth. The loss of a LIM domain in the ceh-14(ch2) allele causes a partial loss-of-function phenotype. Examination of the neurites of ALA and tail neurons using a ceh-14::GFP reporter shows abnormal axonal outgrowth and pathfinding.

Distinct roles of the Gcn5 histone acetyltransferase revealed during transient stress-induced reprogramming of the genome.

BACKGROUND: Gcn5 belongs to a family of histone acetyltransferases (HATs) that regulate protein function by acetylation. Gcn5 plays several different roles in gene transcription throughout the genome but their characterisation by classical mutation approaches is hampered by the high degree of apparent functional redundancy between HAT proteins. RESULTS: Here we utilise the reduced redundancy associated with the transiently high levels of genomic reprogramming during stress adaptation as a complementary approach to understand the functions of redundant protein families like HATs. We show genome-wide evidence for two functionally distinct roles of Gcn5. First, Gcn5 transiently re-localises to the ORFs of long genes during stress adaptation. Taken together with earlier mechanistic studies, our data suggests that Gcn5 plays a genome- wide role in specifically increasing the transcriptional elongation of long genes, thus increasing the production efficiency of complete long transcripts. Second, we suggest that Gcn5 transiently interacts with histones close to the transcription start site of the many genes that it activates during stress adaptation by acetylation of histone H3K18, leading to histone depletion, probably as a result of nucleosome loss as has been described previously. CONCLUSIONS: We show that stress adaptation can be used to elucidate the functions of otherwise redundant proteins, like Gcn5, in gene transcription. Further, we show that normalization of chromatin-associated protein levels in ChIP experiments in relation to the histone levels may provide a useful complement to standard approaches. In the present study analysis of data in this way provides an alternative explanation for previously indicated repressive role of Gcn5 in gene transcription.

Comprehensive characterization of the embryonic factor LEUTX.

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

Homeodomain subtypes and functional diversity.

The homeodomain is a protein domain of about 60 amino acids that is encoded by homeobox genes. The homeodomain is a DNA binding domain, and hence homeodomain proteins are essentially transcription factors (TFs). They have been shown to play major roles in many developmental processes of animals, as well as fungi and plants. A primary function of homeodomain proteins is to regulate the expression of other genes in development and differentiation. Thousands of homeobox genes have been identified, and they can be grouped into many different classes. Often other conserved protein domains are found linked to a homeodomain. Several particular types of homeobox genes are organized into chromosomal clusters. The best-known cluster, the HOX cluster, is found in all bilaterian animals. Tetrapods contain four HOX clusters that arose through duplication in early vertebrate evolution. The genes in these clusters are called Hox genes. Lower chordates, insects and nematodes tend to have only one HOX cluster. Of particular interest is that many of the HOX cluster genes function in the process of pattern formation along the anterior-posterior body axis. Many other types of homeodomain proteins play roles in the determination of cell fates and cell differentiation. Homeobox genes thus perform key roles for all aspects of the development of an organism.

Tracking and characterization of partial and full epithelial-mesenchymal transition cells in a mouse model of metastatic breast cancer.

The various stages of epithelial-mesenchymal transition (EMT) generate phenotypically heterogeneous populations of cells. Here, we detail a dual recombinase lineage tracing system using a transgenic mouse model of metastatic breast cancer to trace and characterize breast cancer cells at different EMT stages. We describe analytical steps to label cancer cells at an early partial or a late full EMT state, followed by tracking their behavior in tumor slice cultures. We then characterize their transcriptome by five-cell RNA sequencing. For complete details on the use and execution of this protocol, please refer to Luond et al. (2021).

DUX4 is a multifunctional factor priming human embryonic genome activation.

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

Dye-filling of the amphid sheath glia: implications for the functional relationship between sensory neurons and glia in Caenorhabditis elegans.

The nervous system is composed of cells including neurons and glia. It has been believed that the former cells play central roles in various neural functions while the latter ones have only supportive functions for neurons. However, recent findings suggest that glial cells actively participate in neural activities, and the cooperation between neurons and glia is important for nervous system functions. In Caenorhabditis elegans, amphid sensory organs in the head also consist of sensory neurons and glia-like support cells (amphid socket and amphid sheath cells). Ciliary endings of some sensory neurons exposed to the environment detect various chemicals, molecules and signals, and the cilia of some neurons can also take up fluorescent dyes such as DiI. Here, we show that the amphid sheath glia are also stained with DiI and that its uptake by the amphid sheath cells correlates with DiI-filling of sensory neurons, suggesting that the amphid sheath glia might interact with sensory neurons. Furthermore, the localization of the amphid sheath cell reporter F52E1.2SP::YFP is abnormal in che-2 mutants, which have defective cilia. These findings imply that sensory neurons might affect amphid sheath glia functions in the amphid sensory organ of C. elegans.

Mitochondrial DNA level, but not active replicase, is essential for Caenorhabditis elegans development.

A number of studies showed that the development and the lifespan of Caenorhabditis elegans is dependent on mitochondrial function. In this study, we addressed the role of mitochondrial DNA levels and mtDNA maintenance in development of C. elegans by analyzing deletion mutants for mitochondrial polymerase gamma (polg-1(ok1548)). Surprisingly, even though previous studies in other model organisms showed necessity of polymerase gamma for embryonic development, homozygous polg-1(ok1548) mutants had normal development and reached adulthood without any morphological defects. However, polg-1 deficient animals have a seriously compromised gonadal function as a result of severe mitochondrial depletion, leading to sterility and shortened lifespan. Our results indicate that the gonad is the primary site of mtDNA replication, whilst the mtDNA of adult somatic tissues mainly stems from the developing embryo. Furthermore, we show that the mtDNA copy number shows great plasticity as it can be almost tripled as a response to the environmental stimuli. Finally, we show that the mtDNA copy number is an essential limiting factor for the worm development and therefore, a number of mechanisms set to maintain mtDNA levels exist, ensuring a normal development of C. elegans even in the absence of the mitochondrial replicase.

Spatio-temporal reference model of Caenorhabditis elegans embryogenesis with cell contact maps.

The nematode Caenorhabditis elegans has been used as a model for developmental biology for decades. Still, the few publicly available spatio-temporal (4D) data sets have conflicting information regarding variability of cell positions and are not well-suited for a standard 4D embryonic model, due to compression. We have recorded six uncompressed embryos, and determined their lineage and 4D coordinates, including nuclear radii, until the end of gastrulation. We find a remarkable degree of stability in the cell positions, as well as little rotational movement, which allowed us to combine the data into a single reference model of C. elegans embryogenesis. Using Voronoi decomposition we generated the list of all predicted cell contacts during early embryogenesis and calculated these contacts up to the approximately 150 cell stage, and find that about 1500 contacts last 2.5 min or longer. The cell contact map allows for comparison of multiple 4D data sets, e.g., mutants or related species, at the cellular level. A comparison of our uncompressed 4D model with a compressed embryo shows that up to 40% of the cell contacts can be different. To visualize the 4D model interactively we developed a software utility. Our model provides an anatomical resource and can serve as foundation to display 4D expression data, a basis for developmental systems biology.

Genome-wide characterisation of the Gcn5 histone acetyltransferase in budding yeast during stress adaptation reveals evolutionarily conserved and diverged roles.

BACKGROUND: Gcn5 is a transcriptional coactivator with histone acetyltransferase activity that is conserved with regard to structure as well as its histone substrates throughout the eukaryotes. Gene regulatory networks within cells are thought to be evolutionarily diverged. The use of evolutionarily divergent yeast species, such as S. cerevisiae and S. pombe, which can be studied under similar environmental conditions, provides an opportunity to examine the interface between conserved regulatory components and their cellular applications in different organisms. RESULTS: We show that Gcn5 is important for a common set of stress responses in evolutionarily diverged yeast species and that the activity of the conserved histone acetyltransferase domain is required. We define a group of KCl stress response genes in S. cerevisiae that are specifically dependent on Gcn5. Gcn5 is localised to many Gcn5-dependent genes including Gcn5 repressed targets such as FLO8. Gcn5 regulates divergent sets of KCl responsive genes in S. cerevisiae and S. pombe. Genome-wide localization studies showed a tendency for redistribution of Gcn5 during KCl stress adaptation in S. cerevisiae from short genes to the transcribed regions of long genes. An analogous redistribution was not observed in S. pombe. CONCLUSIONS: Gcn5 is required for the regulation of divergent sets of KCl stress-response genes in S. cerevisiae and S. pombe even though it is required a common group of stress responses, including the response to KCl. Genes that are physically associated with Gcn5 require its activity for their repression or activation during stress adaptation, providing support for a role of Gcn5 as a corepressor as well as a coactivator. The tendency of Gcn5 to re-localise to the transcribed regions of long genes during KCl stress adaptation suggests that Gcn5 plays a specific role in the expression of long genes under adaptive conditions, perhaps by regulating transcriptional elongation as has been seen for Gcn5 in S. pombe. Interestingly an analogous redistribution of Gcn5 is not seen in S. pombe. The study thus provides important new insights in relation to why coregulators like Gcn5 are required for the correct expression of some genes but not others.

A comprehensive classification and evolutionary analysis of plant homeobox genes.

The full complement of homeobox transcription factor sequences, including genes and pseudogenes, was determined from the analysis of 10 complete genomes from flowering plants, moss, Selaginella, unicellular green algae, and red algae. Our exhaustive genome-wide searches resulted in the discovery in each class of a greater number of homeobox genes than previously reported. All homeobox genes can be unambiguously classified by sequence evolutionary analysis into 14 distinct classes also characterized by conserved intron-exon structure and by unique codomain architectures. We identified many new genes belonging to previously defined classes (HD-ZIP I to IV, BEL, KNOX, PLINC, WOX). Other newly identified genes allowed us to characterize PHD, DDT, NDX, and LD genes as members of four new evolutionary classes and to define two additional classes, which we named SAWADEE and PINTOX. Our comprehensive analysis allowed us to identify several newly characterized conserved motifs, including novel zinc finger motifs in SAWADEE and DDT. Members of the BEL and KNOX classes were found in Chlorobionta (green plants) and in Rhodophyta. We found representatives of the DDT, WOX, and PINTOX classes only in green plants, including unicellular green algae, moss, and vascular plants. All 14 homeobox gene classes were represented in flowering plants, Selaginella, and moss, suggesting that they had already differentiated in the last common ancestor of moss and vascular plants.

Evolution of hedgehog and hedgehog-related genes, their origin from Hog proteins in ancestral eukaryotes and discovery of a novel Hint motif.

BACKGROUND: The Hedgehog (Hh) signaling pathway plays important roles in human and animal development as well as in carcinogenesis. Hh molecules have been found in both protostomes and deuterostomes, but curiously the nematode Caenorhabditis elegans lacks a bona-fide Hh. Instead a series of Hh-related proteins are found, which share the Hint/Hog domain with Hh, but have distinct N-termini. RESULTS: We performed extensive genome searches such as the cnidarian Nematostella vectensis and several nematodes to gain further insights into Hh evolution. We found six genes in N. vectensis with a relationship to Hh: two Hh genes, one gene with a Hh N-terminal domain fused to a Willebrand factor type A domain (VWA), and three genes containing Hint/Hog domains with distinct novel N-termini. In the nematode Brugia malayi we find the same types of hh-related genes as in C. elegans. In the more distantly related Enoplea nematodes Xiphinema and Trichinella spiralis we find a bona-fide Hh. In addition, T. spiralis also has a quahog gene like C. elegans, and there are several additional hh-related genes, some of which have secreted N-terminal domains of only 15 to 25 residues. Examination of other Hh pathway components revealed that T. spiralis - like C. elegans - lacks some of these components. Extending our search to all eukaryotes, we recovered genes containing a Hog domain similar to Hh from many different groups of protists. In addition, we identified a novel Hint gene family present in many eukaryote groups that encodes a VWA domain fused to a distinct Hint domain we call Vint. Further members of a poorly characterized Hint family were also retrieved from bacteria. CONCLUSION: In Cnidaria and nematodes the evolution of hh genes occurred in parallel to the evolution of other genes that contain a Hog domain but have different N-termini. The fact that Hog genes comprising a secreted N-terminus and a Hog domain are found in many protists indicates that this gene family must have arisen in very early eukaryotic evolution, and gave rise eventually to hh and hh-related genes in animals. The results indicate a hitherto unsuspected ability of Hog domain encoding genes to evolve new N-termini. In one instance in Cnidaria, the Hh N-terminal signaling domain is associated with a VWA domain and lacks a Hog domain, suggesting a modular mode of evolution also for the N-terminal domain. The Hog domain proteins, the inteins and VWA-Vint proteins are three families of Hint domain proteins that evolved in parallel in eukaryotes.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis.

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Redox properties and evolution of human glutaredoxins.

Glutaredoxins (Grxs) are glutathione-dependent oxidoreductases that belong to the thioredoxin superfamily catalyzing thiol-disulfide exchange reactions via active site cysteine residues. Focusing on the human dithiol glutaredoxins having a C-X-Y-C active site sequence motif, the redox potentials of hGrx1 and hGrx2 were determined to be -232 and -221 mV, respectively, using a combination of redox buffers, protein-protein equilibrium and thermodynamic linkage. In addition, a nonactive site disulfide was identified between Cys28 and Cys113 in hGrx2 using redox buffers and chemical digestion. This disulfide confers nearly five kcal mol(-1) additional stability by linking the C-terminal helix to the bulk of the protein. The redox potential of this nonactive site disulfide was determined to be -317 mV and is thus expected to be present in all but the most reducing conditions in vivo. As all human glutaredoxins contain additional nonactive site cysteine residues, a full phylogenetic analysis was performed to help elucidate their structural and functional roles. Three distinct groups were found: Grx1, Grx2, and Grx5, the latter representing a highly conserved group of monothiol glutaredoxins having a C-G-F-S active site sequence, with clear homologs from bacteria to human. Grx1 and Grx2 diverged from a common ancestor before the origin of vertebrates, possibly even earlier in animal evolution. The highly stabilizing nonactive site disulfide observed in hGrx2 is found to be a conserved feature within the deuterostomes and appears to be the only additional conserved intramolecular disulfide within the glutaredoxins.

A novel conserved RNA-binding domain protein, RBD-1, is essential for ribosome biogenesis.

Synthesis of the ribosomal subunits from pre-rRNA requires a large number of trans-acting proteins and small nucleolar ribonucleoprotein particles to execute base modifications, RNA cleavages, and structural rearrangements. We have characterized a novel protein, RNA-binding domain-1 (RBD-1), that is involved in ribosome biogenesis. This protein contains six consensus RNA-binding domains and is conserved as to sequence, domain organization, and cellular location from yeast to human. RBD-1 is essential in Caenorhabditis elegans. In the dipteran Chironomus tentans, RBD-1 (Ct-RBD-1) binds pre-rRNA in vitro and anti-Ct-RBD-1 antibodies repress pre-rRNA processing in vivo. Ct-RBD-1 is mainly located in the nucleolus in an RNA polymerase I transcription-dependent manner, but it is also present in discrete foci in the interchromatin and in the cytoplasm. In cytoplasmic extracts, 20-30% of Ct-RBD-1 is associated with ribosomes and, preferentially, with the 40S ribosomal subunit. Our data suggest that RBD-1 plays a role in structurally coordinating pre-rRNA during ribosome biogenesis and that this function is conserved in all eukaryotes.

Comprehensive analysis of gene expression patterns of hedgehog-related genes.

BACKGROUND: The Caenorhabditis elegans genome encodes ten proteins that share sequence similarity with the Hedgehog signaling molecule through their C-terminal autoprocessing Hint/Hog domain. These proteins contain novel N-terminal domains, and C. elegans encodes dozens of additional proteins containing only these N-terminal domains. These gene families are called warthog, groundhog, ground-like and quahog, collectively called hedgehog (hh)-related genes. Previously, the expression pattern of seventeen genes was examined, which showed that they are primarily expressed in the ectoderm. RESULTS: With the completion of the C. elegans genome sequence in November 2002, we reexamined and identified 61 hh-related ORFs. Further, we identified 49 hh-related ORFs in C. briggsae. ORF analysis revealed that 30% of the genes still had errors in their predictions and we improved these predictions here. We performed a comprehensive expression analysis using GFP fusions of the putative intergenic regulatory sequence with one or two transgenic lines for most genes. The hh-related genes are expressed in one or a few of the following tissues: hypodermis, seam cells, excretory duct and pore cells, vulval epithelial cells, rectal epithelial cells, pharyngeal muscle or marginal cells, arcade cells, support cells of sensory organs, and neuronal cells. Using time-lapse recordings, we discovered that some hh-related genes are expressed in a cyclical fashion in phase with molting during larval development. We also generated several translational GFP fusions, but they did not show any subcellular localization. In addition, we also studied the expression patterns of two genes with similarity to Drosophila frizzled, T23D8.1 and F27E11.3A, and the ortholog of the Drosophila gene dally-like, gpn-1, which is a heparan sulfate proteoglycan. The two frizzled homologs are expressed in a few neurons in the head, and gpn-1 is expressed in the pharynx. Finally, we compare the efficacy of our GFP expression effort with EST, OST and SAGE data. CONCLUSION: No bona-fide Hh signaling pathway is present in C. elegans. Given that the hh-related gene products have a predicted signal peptide for secretion, it is possible that they constitute components of the extracellular matrix (ECM). They might be associated with the cuticle or be present in soluble form in the body cavity. They might interact with the Patched or the Patched-related proteins in a manner similar to the interaction of Hedgehog with its receptor Patched.