| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Genetics, Vol. 173, 131-149, May 2006, Copyright © 2006
doi:10.1534/genetics.106.055749
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




* Howard Hughes Medical Institute, Department of Biochemistry and Molecular Biophysics, Columbia University Medical Center, New York, New York 10032,
Howard Hughes Medical Institute, Department of Biochemistry and Molecular Biophysics, Center for Computational Biology and Bioinformatics, Columbia University Medical Center, New York, New York 10032 and
Institute of Neuroscience, University of Oregon, Eugene, Oregon 97403
1 Corresponding author: Columbia University, 701 W. 68th St., HHSC 724, New York, NY.
E-mail: or38{at}columbia.edu
| ABSTRACT |
|---|
|
|
|---|
The nematode Caenorhabditis elegans provides a simple model organism to study the lateralization of nervous system function (HOBERT et al. 2002). Such lateralization can be observed in the chemosensory system of the nematode. The best-studied chemosensory neurons are a group of 12 classes of neurons called the amphid sensory neurons (Figure 1A). Each class consists of one pair of two bilaterally symmetric and morphologically indistinguishable neurons, most of which are chemosensory neurons (Figure 1B). At least two classes of these chemosensory neurons, the AWC odor sensory neuron class and the ASE gustatory neuron class are functionally lateralized, allowing the animal to sense and discriminate different sensory cues with the left and right neuron (PIERCE-SHIMOMURA et al. 2001; WES and BARGMANN 2001) (Figure 1B). Functional lateralization of AWCL/R and ASEL/R correlates with the left/right asymmetric expression of putative chemoreceptors (YU et al. 1997; TROEMEL et al. 1999; CHANG et al. 2004) but the extent of lateralization of chemoreceptor gene expression in these neurons is still unclear. In the AWCL/R sensory neurons, only one left/right asymmetrically expressed chemoreceptor has been reported, a seven-transmembrane receptor (TROEMEL et al. 1999) (Figure1B). In the ASEL/R sensory neurons, a total of four asymmetrically expressed, putative chemoreceptors, which all belong to the family of receptor guanylyl cyclases, were known prior to this study (YU et al. 1997; JOHNSTON et al. 2005) (Figure 1B). gcy-5 and gcy-22 are expressed in the right ASE neuron (ASER), whereas gcy-6 and gcy-7 are expressed exclusively in the left ASE neuron (ASEL).
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
The following transgenes were used: otIs3: Is[gcy-7prom::gfp; lin-15(+)] (CHANG et al. 2003), expressed in ASEL and the excretory canal cell; otIs151: Is[ceh-36prom::dsRed2; rol-6(d)] (JOHNSTON and HOBERT 2003), expressed in ASEL/R and AWCL/R; otIs133: Is[ttx-3promB::rfp; pNC4.2(unc-4(+)] (WENICK and HOBERT 2004), expressed in AIYL/R; oyIs17: Is[gcy-8prom::gfp; lin-15(+)], expressed in AFDL/R; and oyIs51: Is[srh-142prom::rfp; lin-15(+)], expressed in ADFL/R, both gifts from Piali Sengupta.
Sequence analysis:
To identify GCY sequences, the sets of predicted proteins for C. elegans and C. briggsae were obtained from the Sanger Institute (http://www.sanger.ac.uk/Projects/C_elegans/WORMBASE/current/wormpep_download.shtml and ftp://ftp.sanger.ac.uk/pub/wormbase/cbriggsae/cb25.agp8/). Representative GCYs were used as position-specific iterated basic local alignment search tool (PSI-BLAST) queries to search the two proteomes. The HMMER 2.3.2 package (EDDY 1998) was used to construct a hidden Markov model from an alignment of GCY catalytic domains and to search for additional GCYs.
The intracellular regions of the transmembrane GCYs and the complete soluble GCY sequences were aligned with T-coffee version 2.03 (NOTREDAME et al. 2000). Maximum parsimony phylogenetic trees were found via heuristic search with PAUP* version 4.0 beta 10 (SWOFFORD 2003). The trees were generated in 100 repeated searches with random addition of taxa to obtain the starting tree. Robustness of the tree partitions was evaluated by constructing a bootstrap consensus tree with 1000 replicates. The trees were visualized with TreeView version 1.6.6 (PAGE 1996).
Nomenclature of gcy genes and gcy gene predictions:
Most but not all gcy names were previously assigned (http://www.wormbase.org). We named two previously unnamed gcy genes gcy-28 (T01A4.1) and gcy-29 (C04H5.3). Both code for receptor-type proteins. A few gcy genes have been double named in the past. The most current names are (with old names in parentheses): gcy-18 (gcy-26), gcy-20 (gcy-16), gcy-17 (gcy-24), odr-1 (gcy-10), gcy-28 (gcy-38 in WS149), gcy-29 (gcy-39 in WS149).
Since we detected several cases where individual parts of one C. elegans gene are homologous to separate, adjacent predicted C. briggsae genes (gene prediction in Wormbase WS149), we suspected that C. briggsae genes may have been incorrectly predicted. We therefore ran the FGENESH program at http://www.softberry.com (SALAMOV and SOLOVYEV 2000) on chromosomal regions that contained the following predicted C. briggsae genes (from WS149): CBG07423(CBP15915) + CBG07424(CBP15916) + CBG07425(CBP15917), CBG20867(CBP04902) + CBG20868(CBP04903), and CBG19454(CBP11205) + CBG19453(CBP11204). In each of these three cases, we found that FGENESH predicted only one gene, whose product was homologous over its entire length with putative C. elegans orthologs. However, in two cases, the revised gene prediction overlooked exons predicted in the original prediction, which contained homology to GCY proteins. We therefore assembled alternate gene prediction by hand, on the basis of the similarity to known GCY proteins. We named these revised predictions CBP15915*, CBP04902*, and CBP11205* (see supplemental Figures 1 and 2 at http://www.genetics.org/supplemental/). On the basis of primary sequence homology, we also suspect that one of the two C. briggsae orthologs of C. elegans, gcy-35, has been incorrectly predicted as the two separate genes CBG20390 (CBP04780) and CBG20392 (CBP04781). Similarly, CBG10472 (CBP08561) and CBG10474 (CBP08562) likely constitute one gene. However, neither of these two suspicions could be corroborated by FGENESH.
Generation of gfp reporter gene fusions, transgenic animals, and identification of reporter gene expressing cells:
Most reporter genes were created by PCR fusion (HOBERT 2002) and some were generated by subcloning PCR amplicons into pPD95.75 (see Table 1). Primer sequences and resulting transgenic arrays are shown in Table 1. DNA was injected at
1050 ng/µl using either unc-122::gfp or rol-6 as injection marker. Cell identifications were done on the basis of overall cell position and morphology and was significantly aided by the uses of the following four colabeling procedures.
|
| RESULTS |
|---|
|
|
|---|
|
GCY-27 is an unusual receptor-type GCY protein (Figure 2). While containing all the intracellular signature motifs of receptor-type GCYs, the predicted GCY-27 protein is several hundred amino acids shorter than all other predicted receptor-type GCY proteins and does not contain a predicted signal sequence (SS), transmembrane domain (TM), or other extracellular motifs (Figure 2). In the absence of complete cDNA/EST sequences, we cannot rule out a gene prediction error, but we consider this to be unlikely for two reasons: (a) The upstream gene adjacent to gcy-27 is relatively close (Figure 4), leaving little room for such a prediction oversight, particularly since in other gcy genes the extracellular domains are large and composed of many exons; and (b) a C. briggsae ortholog (described in more detail in a later section) also lacks the extracellular domain. While receptor-type GCY proteins that lack a TM and SS domain have been described before (MORTON 2004), they do not appear to contain the protein kinase domain that is present in GCY-27 and all other receptor-type GCY proteins. It is interesting to note that GCY-27 is most closely related to the intracellular domain of ODR-1 (Figure 3A), a transmembrane GCY protein whose extracellular domain (which contains no canonical RFLBP domain) has previously been shown to be dispensable for its function in chemosensory signal transduction (L'ETOILE and BARGMANN 2000).
|
|
|
Notably, in a substantial number of cases, the degree of sequence relation correlates with proximity in the genome sequence (Figure 3B). For example, all six members of one subgroup ("gcy-5 subgroup") are located within an
7-Mb interval on chromosome II, five of these (gcy-1gcy-5) map to an
800-kb interval, and three of them are directly adjacent genes, separated by only a few hundred base pairs (Figure 3B; Figure 4). This suggests that these genes arose by relatively recent gene duplications. Likewise, all seven members of what we term the "gcy-7 subgroup" of gcy genes (Figure 3A) reside on a single chromosome and five of them in a <8-Mb interval (Figure 3B; Figure 4). As we will demonstrate below, the degree of sequence relation and chromosomal location also correlates with similarities in gene expression patterns.
Expression patterns of receptor-type gcy genes:
While the expression of all seven soluble gcy genes has already been described (YU et al. 1997; CHEUNG et al. 2004; GRAY et al. 2004), the expression patterns of only 8 of the 27 receptor-type gcy genes was previously reported, 6 in the context of a preliminary and incomplete genome analysis (gcy-5, gcy-6, gcy-7, gcy-8, gcy-10/odr-1, gcy-12, and gcy-22) and 2 in the course of a functional analysis (odr-1 and daf-11) (YU et al. 1997; BIRNBY et al. 2000; L'ETOILE and BARGMANN 2000; JOHNSTON et al. 2005) (summarized in Table 2). We generated gfp reporter fusions to the putative cis-regulatory regions of the remaining 19 gcy genes (schematically shown in Figure 4). In all except one case, sequences to the next upstream genes were included in the reporter gene fusions (Figure 4). While such upstream sequences are most likely to harbor gene regulatory elements, it needs to be kept in mind that additional regulatory elements may be located elsewhere and that therefore the expression patterns of the gfp reporter genes can provide only a first approximation of endogenous gene expression profiles. Since gcy genes are likely to primarily function as chemoreceptors in the mature nervous system, we restricted our gene expression analysis of transgenic animals that harbor the respective reporter genes to larval and adult stages. To allow us to reliably detect expression in the ASEL/R gustatory neurons, the neurons in which our laboratory is most interested in, we used a transgenic reporter array in the background in which the ASEL/R neurons are labeled with DsRed2 (otIs151; see MATERIALS AND METHODS), thereby enabling us to score for an overlap between green and red fluorescent signals. In most cases, multiple lines were analyzed per construct (see Table 1) and we did not observe any notable differences between individual lines. Reporter gene expression was observed for each of the gcy genes examined. Expression patterns are shown in Figure 5 and summarized in Table 2, which also lists previously described expression patterns. The expression patterns of all receptor-type gcy genes can be summarized as follows.
|
|
Cell types:
The expression of receptor-type gcy genes is strongly biased toward the nervous system. With the exception of two non-neuronally expressed gcy genes (gcy-9 and gcy-11), all gcy genes are expressed in a restricted subset of neurons. Of 27 genes, 21 are exclusively expressed in the nervous system and 4 are expressed in restricted sets of both neuronal and non-neuronal cells (Table 2; Figure 5).
Expression within the nervous system:
Within the nervous system, most but not all, C. elegans gcy genes are expressed in sensory neurons. All except two pairs (ASH and ADF) of the 12 amphid sensory neuron classes shown in Figure 1A express at least one gcy gene. Unlike what the previous analyses of eight gcy genes seemed to indicate (YU et al. 1997; BIRNBY et al. 2000; L'ETOILE and BARGMANN 2000; JOHNSTON et al. 2005), neuronal gcy gene expression is, however, not restricted to sensory neurons. Nonsensory neurons that express gcy genes include the AIY, AIM, AVK, RIA, and PVT interneuron classes and the RIM neuron, which is both a motor neuron and an interneuron. In addition, as mentioned above, many neurons in the nervous system, including, for example, ventral cord motor neurons, express the widely expressed gcy-28 gene. As summarized in Table 3, 25 of 27 of the C. elegans receptor-type gcy genes are expressed in the nervous system, 25 of 27 are expressed in various types of sensory neurons plus other neurons, 15 of 27 are expressed exclusively in sensory neurons, and 9 of 27 are restricted to single neuron classes. The expression of mammalian gcy genes, of which there are only 7, show roughly comparable patterns; some are expressed in non-neuronal cells and those that are expressed in neurons are strongly biased toward expression in sensory structures (WEDEL and GARBERS 2001).
|
|
gcy genes that are asymmetrically expressed in either ASER or ASEL are not asymmetrically expressed in other sensory neurons. For example, the gcy-1 and gcy-3 genes, two ASER-expressed genes, are expressed bilaterally in other pairs of neurons.
In addition to the asymmetrically expressed gcy genes, we found two bilaterally, albeit weakly, expressed gcy genes in ASEL/R, gcy-19 and gcy-29, increasing the percentage of ASE-expressed gcy genes to 41% (11/27) of all gcy genes. The possible lack of regulatory elements in gfp reporter genes may lead to the oversight of perhaps even more ASE-expressed gcy genes.
We have not observed any other obvious left/right asymmetric gcy gene patterns in bilaterally symmetric neurons, including the AWCL/R neurons (which we could easily identify with the DsRed2-expressing otIs151 transgene). The AWCL/R neurons are the only other known neuron pair displaying functional laterality (WES and BARGMANN 2001) and although they express five gcy genes (Table 4), none of them is obviously lateralized.
Similarity of chromosomal position, primary sequence, and gene expression patterns:
Strikingly, the coexpression of gcy genes in ASEL or ASER correlates extensively with the primary sequence similarity and chromosomal location of the gcy genes. The ASER-expressed gcy-1, gcy-3, gcy-4, and gcy-5 genes fall into one sequence subgroup (Figure 3A) and localize within
800 kb and two of them are directly adjacent to one another (Figure 3B; Figure 4). Notably, their coexpression is not due simply to a joint cis-regulatory element; as our reporter constructs clearly demonstrate (Figure 4), each gene contains separable cis-regulatory elements.
Similarly, all four ASEL-expressed gcy genes (gcy-6, gcy-7, gcy-14, gcy-20) fall into one sequence subgroup (Figure 3A) and localize at the center of chromosome V (Figure 3B). These observations argue that gene duplication events that lead to the generation of these paralogs also duplicated their regulatory regions. The only exception to this pattern is the gcy-22 gene, which is predominantly expressed in ASER, but, in terms of primary sequence identity, clusters more closely with the ASEL-expressed genes than with the ASER-expressed genes (Figure 3).
A similar correlation of sequence relation and gene expression can also be observed in non-ASE-expressed gcy genes. The gcy-8, gcy-18, and gcy-23 genes are related by sequence (Figure 3A), localize within a 15.6-Mb interval (Figure 3B), and are all coexpressed in the AFD sensory neurons (Figure 5). Additionally, the AFD-expressed gcy-29 gene is closely related to gcy-8, gcy-18, and gcy-23 by primary sequence, but it localizes to a different chromosome. Finally, the sequence-related gcy-15 and gcy-21 genes (Figure 3A) are in close chromosomal proximity (
530 kb; Figure 3B) and are coexpressed in the ASG chemosensory neurons (Figure 5).
Taken together, all these similarities indicate that related gcy genes arose by local gene duplication events that duplicated not only the protein-coding region but also their cis-regulatory control regions.
Regulation of left/right asymmetric expression of gcy genes:
How is the laterality of the newly characterized ASEL/R-expressed gcy genes controlled? We have previously identified a complex network of transcription factors and microRNAs (miRNAs) that control ASE laterality after animals have passed through an initial hybrid precursor state (CHANG et al. 2003, 2004; JOHNSTON and HOBERT 2003, 2005; JOHNSTON et al. 2005). These regulatory factors fall into two broad categories: (1) factors that control the activity of a bistable feedback loop that determines whether an ASE neuron adopts the ASEL or ASER state and (2) factors that act outside the regulatory loop to determine specific subsets of terminal differentiation features. To test whether the newly identified, asymmetrically expressed gcy genes are subject to regulation by this network or controlled by a different set of regulatory factors, we analyzed gfp reporter expression profiles in two null mutant backgrounds that are representative for each category, lsy-6 and lim-6. Animals that lack the lsy-6 miRNA display a complete switch from the ASEL fate to the ASER fate (JOHNSTON and HOBERT 2003) and animals that lack the lim-6 LIM homeobox gene fail to activate a defined subset of ASEL features and fail to repress a subset of ASER features (HOBERT et al. 1999; JOHNSTON et al. 2005). We find that all newly identified gcy genes are components of the ASEL and ASER state that is controlled by the lsy-6-dependent regulatory feedback loop. The ASER-specific gcy genes are ectopically activated in ASEL upon loss of the ASEL inducer lsy-6, and the ASEL-specific gcy genes are lost in the ASEL neuron of lsy-6 null mutant animals (Figure 6). We note that in the case of the ASEL-biased gcy-14 gene, which is normally strongly expressed in ASEL and weakly expressed in ASER, the conversion of the ASEL to the ASER fate does not entail a loss of gcy-14 expression, but rather the bilateralization of weak expression (Figure 6B).
How does the lim-6 LIM homeobox gene contribute to the expression of the newly identified gcy genes? Similarly to the previously described ASEL-specific gcy-6 and gcy-7 genes, we found that the ASEL bias of gcy-14 and gcy-20 are unaffected by loss of lim-6 (Figure 6B), demonstrating that ASEL inducers such as lsy-6 independently regulate the expression of lim-6 and ASEL-specific gcy genes (Figure 6C). Similarly to the previously described ASER-specific gcy-5 and gcy-22 genes, we found that lim-6 represses the expression of the ASER-specific gcy-1, gcy-3, and gcy-4 genes in ASEL (Figure 6, A and B). In summary, all known asymmetrically expressed gcy genes are controlled by the same categories of gene regulatory factors.
Functional analysis of gcy-5:
Three putative loss-of-function alleles have been generated by the C. elegans knockout consortia in one of the nine asymmetrically expressed gcy genes, gcy-5 (Figure 7A). At least one of them is a putative molecular null allele (see legend to Figure 7A). We analyzed all three mutant gcy-5 alleles in a chemotaxis assay that measures the functionality of the ASE neurons (BARGMANN and HORVITZ 1991). Our previous work has demonstrated that the ASEL neuron preferentially senses sodium, but not chloride, whereas the ASER neuron senses chloride and weakly contributes to sodium sensation (PIERCE-SHIMOMURA et al. 2001; CHANG et al. 2004). Neither gcy-5 mutant allele shows any significant defects in attraction to sodium or chloride (Figure 7B). We conclude that gcy-5 is not required for generic aspects of ASER development or function.
|
100 million years ago (STEIN et al. 2003), provided us with the opportunity to examine the evolution of the gcy gene family. To identify the complete set of guanylyl cyclase genes in C. briggsae, we again employed PSI-BLAST to search the latest release of the complete genome databases of C. briggsae, using a set of known GCY proteins as queries. We identified a total of 33 gcy genes in C. briggsae, one less than in C. elegans (Figure 3A). Compared to C. elegans, C. briggsae contains one additional soluble gcy gene (total = 8) and two fewer receptor-type gcy genes (total = 25). A total of 23 nematode receptor and soluble gcy genes show clear one-to-one ortholog matches between C. elegans and C. briggsae (Figure 3A). All the other gcy genes show species-specific gene duplication events. For example, the C. elegans gcy-6 gene has duplicated to produce two paralogs, while the directly adjacent C. elegans gcy-1, gcy-2, and gcy-3 genes are paralogs with only a single ortholog in C. briggsae (Figure 3A). We have mentioned above the existence of a receptor-type GCY protein, GCY-27, that lacks a predicted extracellular domain. C. elegans GCY-27 is closely related to the intracellular domain of the receptor-type ODR-1; both proteins are more closely related to each other than to the two closest C. briggsae homologs (Figure 3A). One of these homologs, CBP15915* (see MATERIALS AND METHODS for refined gene predictions), is clearly orthologous to ODR-1 due to synteny in the chromosomal region; that is, all predicted proteins in the direct neighborhood of ODR-1 are homologous to the neighboring proteins of the putative C. briggsae ODR-1 ortholog. The protein designated CBP11205* (see MATERIALS AND METHODS for refined gene prediction) is the likely ortholog of GCY-27. Both proteins share the unusual feature of not containing an extracellular domain. As in the case of GCY-27, the absence of an extracellular domain in the C. briggsae ortholog is not likely caused by a failure to predict more exons because neighboring genes are in close proximity. We conclude that the extracellular domain of an ancestral protein was lost before the C. briggsae/C. elegans split, giving rise to GCY-27 and its ortholog CBP11205*.
We sought to investigate the issue of orthology and paralogy in more detail by (a) analyzing the degree of synteny among gcy genes and (b) by determining the expression patterns of some orthologous C. briggsae genes. For the syteny analysis, we focused on one subgroup of left/right asymmetrically expressed gcy subfamilies, the "gcy-5 subfamily" (Figure 3A).
We found that the chromosomal arrangement of members of the gcy-5 subfamily differ significantly between C. elegans and C. briggsae. The gcy-19 orthologs are located in distinct environments in C. elegans and C. briggsae (Figure 8A), suggesting a translocation event after the lineage split. The case of the gcy-1, gcy-2, gcy-3, gcy-4, and gcy-5 genes is more complex and provides a fascinating glimpse into evolutionary divergence. While gcy-4 and gcy-5 are direct neighbors in both C. elegans and C. briggsae, they are translocated together with their neighbors to a distinct region in C. briggsae (Figure 8B). Curiously, this distinct region contains the single C. briggsae ortholog of three C. elegans paralogs from the gcy-5 subfamily, gcy-1, gcy-2, and gcy-3 (Figure 8B). Since all paralogous members of the gcy-5 subgroup likely arose by gene duplication, the C. briggsae gene organization is probably more reflective of the ancestral situation than C. elegans. One conceivable scenario is that a common ancestor of C. elegans and C. briggsae contained, like C. briggsae, three adjacent gcy genes (which themselves arose by gene duplication). In the C. elegans lineage, this cluster split up, with gcy-4 and gcy-5 translocating to a distinct chromosomal region and the gcy-1/2/3 ortholog also translocating to a distinct location and then subsequently duplicating to generate gcy-1, gcy-2, and gcy-3 in the C. elegans lineage.
|
We found that C. briggsae gcy-4prom::gfp is expressed bilaterally in both ASEL and ASER (Figure 9A). This is in striking contrast to C. elegans gcy-4prom::gfp, whose expression is biased to ASER in C. elegans (Figure 5C; Figure 6B). Such a difference could be caused by two different mechanisms. The gcy-4 loci may contain distinct cis-regulatory information in C. elegans and C. briggsae, or, alternatively, the cis-regulatory information may be the same, but different trans-acting factors in C. elegans and C. briggsae interpret this information differentially. To distinguish between these possibilities, we injected the C. briggsae gcy-4prom::gfp construct into C. elegans. If the C. elegans and C. briggsae gcy-4 reporter constructs contained the same cis-regulatory information, C. briggsae gcy-4prom::gfp should be expressed in the same pattern in C. briggsae and C. elegans (bilateral expression in ASEL and ASER). In contrast, if there are differences in trans-acting factors, then the C. briggsae gcy-4prom::gfp reporter should be expressed in C. elegans in a manner similar to that of C.elegans gcy-4prom::gfp (biased to ASER). We found that C. briggsae gcy-4prom::gfp expression becomes ASER biased when injected into C. elegans, which mimics the expression of C. elegans gcy-4prom::gfp (Figure 9A). The difference between gcy-4 expression in C. elegans and gcy-4 expression in C. briggsae therefore does not appear to reflect a difference in their cis-regulatory architecture but rather indicates differences in the trans-acting factors that control gcy-4 expression.
|
| DISCUSSION |
|---|
|
|
|---|
The function of nematode receptor-type gcy genes:
C. elegans and C. briggsae contain an unusual number of receptor-type gcy genes. Insects such as Drosophila melanogaster or Anopheles gambiae contain six receptor-type guanylyl cyclases (MORTON 2004), mammals contain seven (four orphan and three peptide-binding receptors) (LUCAS et al. 2000; WEDEL and GARBERS 2001), but C. elegans contains 27 and C. briggsae 25 (this study). The physiological function of insect gcy genes is entirely unknown, although the expression of the only two analyzed receptors in sensory neurons (among other neurons) has been noted (MORTON 2004). Vertebrate gcy genes are expressed in several different tissue types, including chemosensory neurons (WEDEL and GARBERS 2001).
We propose that the significant expansion of receptor-type gcy genes in the nematode lineage is a reflection of their employment as chemoreceptors used to assess and navigate through their natural habitat. This hypothesis, which was also put forward by YU et al. (1997), is mainly based on the observation that almost 90% (24/27) of gcy genes are expressed in sensory neurons (lack of sensory neuron expression of three gcy genes may be caused merely by a lack of the complete set of cis-regulatory elements in the reporter genes used). Moreover, 41% (11/27) of gcy genes are expressed in the main gustatory neuron class of C. elegans, ASE. This neuron class has previously been shown to be functionally lateralized in that it can sense different chemosensory cues (PIERCE-SHIMOMURA et al. 2001). Consistent with a role of gcy genes as chemoreceptors, we find that 9/11 ASE-expressed gcy genes are expressed in a left/right asymmetric manner, thereby providing molecular correlates to functional lateralization.
Amino acids are among the several classes of chemicals that can be sensed by the ASE neurons (BARGMANN and HORVITZ 1991). The presence of a domain in the extracellular parts of many GCY proteins that is homologous to bacterial amino-acid-binding proteins, the RFLBR domain, makes receptor-type GCY proteins good candidates for amino acid receptors. More sensitive assays (WICKS et al. 2000; FAUMONT et al. 2005; MILLER et al. 2005; FAUMONT and LOCKERY 2006) than those used in the initial large scale survey of chemosensory cues (WARD 1973) will be required to establish the full spectrum of amino acids and other possible sensory cues that signal through ASE. Such a systematic cataloging of sensory cues needs to be followed by a systematic analysis of strains harboring deletions in gcy genes or misexpressing gcy genes to establish their roles as amino acid receptors.
A role for GCY proteins as salt receptors is also conceivable but highly speculative at present. There is as yet no consensus about the molecular identity of salt receptors in the vertebrate gustatory system. Notably, the crystal structure of the ANP receptor, a mammalian GCY protein that, like most C. elegans GCY proteins, contains an extracellular RFLBR domain, revealed a high-affinity chloride-binding site (VAN DEN AKKER et al. 2000).
Our attempt to establish a mutant phenotype for a gcy gene, the ASER-expressed gcy-5 gene, has failed so far, but we note that we have tested only one of the two known cues sensed in an ASER-specific manner, namely chloride ions. The failure to detect a mutant phenotype, however, does allow the conclusion that gcy genes such as gcy-5 are unlikely to control a fundamental, nonredundant aspect of the development or overall function of the neuron. It appears more likely that gcy-5 and, by inference, other ASE-expressed gcy genes fulfill a sensory-modality-specific function in ASE, such as being a receptor for a specific class of gustatory cues.
Whereas more than half of C. elegans gcy genes (15/27) are expressed exclusively in sensory neurons, one-third are also expressed in nonsensory neurons, including interneurons and non-neuronal cells, suggesting that GCY proteins also respond to endogenously produced ligands. Since all known ligands for mammalian GCY proteins are peptidergic-signaling molecules (LUCAS et al. 2000) and since C. elegans contains scores of neuropeptide-encoding genes (LI et al. 1999), we propose that nonsensory C. elegans GCY proteins may be receptors for peptidergic ligands.
The function of receptor-type GCY proteins may not be restricted to a receptor function. The C. elegans ODR-1 protein does not require its extracellular domain to fulfill its function in transducing odorsensory signals (L'ETOILE and BARGMANN 2000). Moreover, one GCY protein that we describe here, the ODR-1-related GCY-27 protein, entirely lacks an extracellular domain. These proteins may heterodimerize with ligand-binding GCY receptors to constitute a receptor complex and/or they may serve as second-messenger-producing signaling proteins that are embedded in signal transduction cascades triggered by other receptor systems.
Coexpression of gcy genes:
Another notable feature of gcy gene expression patterns is their degree of coexpression. Six gcy genes are coexpressed in the ASI sensory neuron class, five are coexpressed in ASER, four in ASEL, five in the AWC olfactory neuron class, four in the AFD thermosensory neuron class, and two gcy genes are each coexpressed in the ASG and PHA phasmid sensory neuron classes. In addition, several gcy genes expressed in multiple cell types show similar and nonintuitive combinations of expression patterns. gcy-7 and gcy-20 are coexpressed in ASEL and the excretory cell (two cells of no obvious relation) and daf-11 and odr-1 show a precisely overlapping expression pattern in five amphid sensory neurons (BIRNBY et al. 2000; L'ETOILE and BARGMANN 2000). Coexpression of receptor-type GCY proteins raises at least two different possibilities:
Laterality in the nematode nervous system:
With the identification of a total of nine asymmetrically expressed gcy genes, the ASE neurons present so far the most striking example of a lateralized chemosensory neuron class. As mentioned above, previous work has provided a conceptual framework for the functional relevance of laterality in the ASE neurons (PIERCE-SHIMOMURA et al. 2001). We expect that other neuron classes may similarly employ the principle of lateralizing chemosensory function. However, our expression analysis has not revealed further examples of laterality in gcy gene expression profiles in other bilaterally symmetric amphid neurons, thereby leaving only the AWC odorsensory and ASE gustatory neurons as neuron pairs with lateralized functions (Figure 1B). Considering the potential usefulness of lateralizing chemoreceptor function, it is our expectation that the analysis of expression of the hundreds of chemoreceptors of the seven-transmembrane receptor family (BARGMANN 1998) will reveal more examples of laterality in the nervous system.
Our analysis of the laterality of gcy gene expression in C. briggsae revealed several striking examples of evolutionary plasticity in the laterality of the gustatory system of nematodes. The variation that we observe appears to be caused by a variation in both cis-acting elements and trans-acting factors, a conclusion that we can draw from our comparison of C. briggsae gcy promoter activity in C. briggsae and C. elegans. Variations in cis-regulatory control have been recognized as a major feature of evolutionary processes (CARROLL et al. 2001) but our cross-species analysis also provides strong support for the more conventional view of the evolution of trans-acting factors. The clearly distinct nature of cis-regulatory architecture of at least some gcy genes in C. elegans and C. briggsae is a strong reminder that the use of bioinformatic tools that use phylogenetic conservation to identify cis-regulatory elements in genomic sequences (e.g., BIGELOW et al. 2004; for review see BULYK 2003) may easily lead to false-negative predictions. It is difficult at this point to speculate about differences in the trans-acting factors that control laterality in C. elegans and C. briggsae. We have demonstrated that so far all known left/right asymmetrically expressed terminal differentiation markers in the ASEL/R neurons, including all gcy genes, are under control of a previously described bistable feedback loop that is composed of several transcription factors and miRNAs (JOHNSTON et al. 2005). The activity of transcription factors that are controlled by the output of the loop, for example, the lim-6 LIM homeobox gene (Figure 6B), may have functionally diverged in C. briggsae. This can perhaps be best illustrated with the gcy-4 gene, which in C. elegans is repressed by lim-6 in ASEL. In C. briggsae, gcy-4 is not repressed in ASEL, yet the C. briggsae gcy-4 promoter can be repressed in C. elegans. These observations are consistent with the presence of a lim-6 responsive repressor element in both C. elegans and C. briggsae gcy-4 promoters but with an inability of C. briggsae lim-6 to control this element in C. briggsae.
Mutant screening approaches that identify the complete set of trans-acting factors controlling left/right asymmetric gcy gene expression in C. elegans and experimental promoter dissection approaches that identify cis-regulatory elements through which these trans-acting factors act are currently ongoing in our laboratory and are likely to reveal the molecular control of laterality in the gustatory system of C. elegans and its divergence in C. briggsae.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|