The cys-loop ligand-gated ion channel gene superfamily of the nematode, Caenorhabditis elegans

The nematode, Caenorhabditis elegans, possesses the most extensive known superfamily of cys-loop ligand-gated ion channels (cys-loop LGICs) consisting of 102 subunit-encoding genes. Less than half of these genes have been functionally characterised which include cation-permeable channels gated by acetylcholine (ACh) and γ-aminobutyric acid (GABA) as well as anion-selective channels gated by ACh, GABA, glutamate and serotonin. Following the guidelines set for genetic nomenclature for C. elegans, we have designated unnamed subunits as lgc genes (ligand-gated ion channels of the cys-loop superfamily). Phylogenetic analysis shows that several of these lgc subunits form distinct groups which may represent novel cys-loop LGIC subtypes.


Introduction
Members of the cysteine-loop ligand-gated ion channel (cys-loop LGIC) superfamily are neurotransmitter receptors that mediate synaptic transmission in vertebrates and invertebrates. Cys-loop LGICs are made up of five homologous subunits arranged around a central ion channel (Sine and Engel 2006). The characteristic cys-loop motif is situated in the N-terminal extracellular ligand-binding region of each subunit and consists of two disulphide-bond forming cysteines separated by 13 amino acid residues, several of which are highly conserved to form a signature sequence (Gene Ontology ID GO:0005230). The cys-loop plays roles in receptor assembly (Green and Wanamaker 1997) and gating of the ion channel (Grutter et al. 2005).
The subunit composition determines the functional and pharmacological properties of the cys-loop LGIC, thus receptor diversity is generated by multiple subunit-encoding genes. For instance, the human cys-loop LGIC superfamily consists of 45 genes encoding cation-permeable nicotinic acetylcholine receptors (nAChRs) (Kalamida et al. 2007), cation-permeable serotonin or 5-hydroxytryptamine type 3 (5-HT 3 ) receptors (Reeves and Lummis 2002), anion-permeable c-aminobutyric acid (GABA) type A and C (GABA A and GABA C ) receptors (Bormann 2000;Darlison et al. 2005) and anion-permeable glycine receptors (Kirsch 2006). As demonstrated in the fruit fly (Drosophila melanogaster), honey bee (Apis mellifera) and red flour beetle (Tribolium castaneum), insect cys-loop LGIC superfamilies are more compact, consisting of just over 20 subunits Littleton and Ganetzky 2000). Approximately 15 of these genes are known to encode cation-permeable nAChRs , anion and possibly cation channels gated by GABA (Buckingham et al. 2005), glutamate-gated anion channels (Vassilatis et al. 1997) and histamine-gated anion channels (Gisselmann et al. 2002;Zheng et al. 2002). The largest known cys-loop LGIC superfamily belongs to the simple nematode, Caenorhabditis elegans, where an early analysis of its completed genome revealed 90 ligand-gated ion channel genes (Bargmann 1998). In common with vertebrates, members of this large cys-loop LGIC superfamily include cationpermeable nAChRs (Jones and Sattelle 2004) and anionpermeable GABA receptors (Schuske et al. 2004). As is the case for insects, C. elegans has glutamate-gated anion channels (Wolstenholme and Rogers 2005) but so far no nematode histamine-gated ion channels have been found. However, C. elegans possesses cys-loop LGIC receptors not identified in vertebrates and insects which include anion channels gated by acetylcholine (Putrenko et al. 2005) and serotonin (Ranganathan et al. 2000). The considerable diversity of receptor subtypes in C. elegans may be broadened further as many cys-loop LGIC subunits have yet to be characterised. Here, we present an update on the genes present in the C. elegans cys-loop LGIC superfamily.

Methods
All C. elegans protein sequences bearing motifs particular to cys-loop LGICs were assembled from WormBase (http:// www.wormbase.org release WS185 Dec 23 2007). These motifs and their source databases are neurotransmitter-gated ion-channel ligand-binding INTERPRO:IPR006202, neurotransmitter-gated ion-channel INTERPRO:IPR006201, neurotransmitter-gated ion-channel ligand binding domain PFAM:PF02931, neurotransmitter-gated ion-channel transmembrane region PFAM:PF02932 and ligand-gated ion channel KOG3644. Cys-loop and TM2 sequences were identified from an alignment of protein sequences of whole subunits constructed with ClustalX (Thompson et al. 1997) using the slow-accurate mode with a gap opening penalty of 10 and a gap extension penalty of 0.1 as well as applying the Gonnet 250 protein weight matrix (Benner et al. 1994). The protein alignment was viewed using GeneDoc (http://www. nrbsc.org/gfx/genedoc/index.html). The neighbour-joining method (Saitou and Nei 1987), available with the ClustalX program, was used to construct a phylogenetic tree, which was then displayed using the TreeView application (Page 1996).

Results and discussion
Our analysis of the C. elegans genome revealed the largest cys-loop LGIC superfamily so far characterised consisting of 102 genes. As we have previously reported the C. elegans nAChR and orphan subunits (Brown et al. 2006;Jones and Sattelle 2004;Mongan et al. 1998), we concentrate here on the remaining cys-loop LGIC genes. For unnamed subunit genes, we continued to use the nomenclature first used for orphan nAChR subunits ). Thus, 24 genes were designated lgc for ligand-gated ion channel of the cys-loop LGIC superfamily (Table 1). As shown in Fig. 1, the C. elegans cys-loop LGIC subunits were grouped according to sequence homology. As with the C. elegans nAChR subunits (Mongan et al. 1998), each group was named, if possible, after a well characterised subunit. For a detailed comparison of C. elegans cys-loop LGIC subunits with those of other organisms, see (Dent 2006).

Ungrouped subunits
Three subunits, LGC-32, LGC-33 and LGC-34, were not placed into any groups as they are highly divergent (Fig. 1), only showing up to 10, 10 and 15% identity respectively with other C. elegans, human and D. melanogaster cys-loop LGICs. These subunits have not been functionally characterised so it is unknown what ligands act on them. The sequence of LGC-34 and a close homologue in the parasitic nematode, Dirofilaria immitis, has been previously described where it has been noted that even though the subunit possesses features common to cys-loop LGICs, such as conserved residues in the N-terminal extracellular region and four transmembrane domains, the cys-loop is absent (Yates and Wolstenholme 2004). Interestingly, a bacterial protongated ion channel (Glvi) resembling cys-loop LGICs but also lacking the cys-loop has been recently reported (Bocquet et al. 2007;Tasneem et al. 2005).
LGC-34 may thus represent more ancestral members of the cys-loop LGIC superfamily. Unlike Glvi, LGC-34 possesses two cysteine residues which are conserved in glutamate-, histamine-and vertebrate glycine-gated anion channels (Dent 2006;Yates and Wolstenholme 2004). Known as the second cys-loop, these two cysteines flank loop C which is a region that contributes to ligand binding (Corringer et al. 2000).

The EXP-1 Group
This group consists of only two subunits, one of which is EXP-1, a GABA-gated cation channel that mediates enteric muscle contraction (Beg and Jorgensen 2003). While EXP-1 overall resembles more closely anion permeable cys-loop LGICs such as UNC-49 as opposed to cation-permeable nAChRs, it lacks the PAR motif preceding the second transmembrane domain (TM2) which is important for anion selectivity (Jensen et al. 2005). Instead it possesses the residues ETE where the presence of glutamic acid residues has been implicated in determining cation selectivity (Jensen et al. 2005;Wotring and Weiss 2008). The other subunit in this group, LGC-35, also has a glutamate residue preceding TM2 (Table 1) and thus may be a cation channel.

The AVR-14 Group
This group consists of six subunits (Fig. 1) that make up glutamate-gated chloride channels that are targeted by the avermectin and milbemycin anthelmintics (Wolstenholme and Rogers 2005). Examples of GluCl function include AVR-15 and GLC-2 in the inhibition of pharyngeal pumping (Dent et al. 1997;Laughton et al. 1997;Pemberton et al. 2001) and AVR-14, AVR-15, GLC-1 and GLC-3 in regulating locomotion (Cook et al. 2006). Recently, it has been shown GLC-3 plays a role in olfactory behaviour (Chalasani et al. 2007).

The UNC-49 Group
Products of the unc-49 gene form a GABA-gated anion channel that mediates body muscle inhibition during locomotion (Bamber et al. 1999;Richmond and Jorgensen 1999). There are four other subunits included in this group (Fig. 1).
With regard to other species, members of the UNC-49 group are clearly most closely related to mammalian and insect GABA-gated anion channels (Dent 2006). GAB-1 forms functional GABA receptors when coexpressed in Xenopus laevis oocytes with either HG1A or HG1E, which are putative GABA receptor subunits from the parasitic nematode Haemonchus contortus (Feng et al. 2002). The remaining subunits have yet to be characterised but it is worth noting that LGC-36 possesses the amino acids AER instead of PAR before TM2 (Table 1). The presence of a glutamic acid residue at this position and the absence of a proline may result in a cation-selective channel (Jensen et al. 2005). LGIC subunits that have been given names in this report. Amino acid residues preceding TM2 that are important for determining the charge selectivity of the ion channel (Jensen et al. 2005) are given. The cys-loop sequence is shown with the cysteine residues highlighted in black shading while gray shading indicates conserved residues. The presence of a putative second cys-loop (Dent 2006 LGC-50 T20B12.9 MOD-1 SAR LGC-51 F12B6.3 GGR-3 PAR LGC-52 Y73F8A.2 GGR-3 PAR LGC-53 T21F2.1 GGR-3 PAR LGC-54 T15B7.16 GGR-3 PAR LGC-55 Y113G7A.5 GGR-3 PAR Invert Neurosci (2008) 8:41-47 43 The GGR-1 Group This group includes six subunits (Fig. 1), the functional characterisation of which has so far not been published. The neurotransmitters/ligands to which GGR-1 group subunits respond to cannot be easily inferred through phylogenetic analysis. For example, when comparing with human cys-loop LGICs they bear closest resemblance to glycine a subunits with approximately 25% identity while LGC-42 is slightly more similar to histamine-gated anion channels with 26% identity when comparing with insect cys-loop LGICs. Members of the GGR-1 group have been previously denoted as ''ce_Group I'' subunits which most closely resemble insect histamine-gated anion channels (Dent 2006). It would be of interest to determine if members of the GGR-1 group respond to glycine, histamine or other ligands not so far known to act on C. elegans cysloop LGICs. The subunits of the GGR-1 group are likely to be anion-selective as they all possess the PAR motif before TM2 (Table 1) which is important for anion selectivity (Jensen et al. 2005). The only exception is GGR-1 which instead has the amino acid residues PGR. With the proline still present and an absence of an acidic residue (particularly glutamic acid) in this region, it is likely that GGR-1 will also be anion selective.

The LGC-45 Group
This group consists of three subunits (Fig. 1), the functional characterisation of which has so far not been published. As is the case for members of the GGR-1 group, LGC-43, LGC-44 and LGC-45 most closely resemble glycine a subunits and histamine-gated anion channels when compared with human and insect cys-loop LGICs respectively. However, unlike the GGR-1 group, the LGC-45 group subunits are likely to be cation selective as they lack the PAR motif before TM2 (Table 1) which is important for anion selectivity (Jensen et al. 2005). In particular, LGC-44 and LGC-45 possess glutamic acid instead of the proline residue which is likely to result in a cation-selective channel (Wotring and Weiss 2008). While LGC-44 shows notable sequence identity with other cys-loop LGIC subunits (up to 37% with LGC-45), it lacks the cys-loop (Table 1), thus it may represent an ancestral cys-loop LGIC subunit as has been noted for LGC-34.

Nicotinic acetylcholine receptors
Orphan Group

GGR-1 Group
LGC-45 Group ACC-1 Group Fig. 1 Tree showing the cysloop LGIC gene superfamily of C. elegans. This figure focuses on cys-loop LGIC subunits other than nAChRs although representative subunits of the nAChR groups are shown. For the complete C. elegans nAChR gene family, see ). Based on protein sequence homology, the subunits are divided into eight groups with three highly divergent subunits, LGC-32,

GGR-3 Group
LGC-33 and LGC-34, not belonging to any group et al. 2005). Acetylcholine-gated anion channels have also been characterised from the snail, Lymnae stagnalis, which appear to have evolved from cation channels through amino acid substitutions in the ion channel pore (van Nierop et al. 2005). In contrast, acetylcholine-gated anion channels of the ACC-1 group may have arisen from substitutions in the ligand-binding domain of anion channels (Putrenko et al. 2005;van Nierop et al. 2005).

The GGR-3 Group
This group includes six subunits (Fig. 1), the functional characterisation of which has so far not been published. All six subunits are likely to form anion-selective channels as they have the PAR amino acid motif before TM2 (Table 1) (Jensen et al. 2005). When comparing with human cys-loop LGICs they bear closest resemblance to GABA a, c and e subunits with approximately 22% identity while they most closely resemble histamine-gated chloride channels (21% identity) when comparing with insect cys-loop LGICs. As with the GGR-1 group, it would be of interest to determine if the GGR-3 group consists of ion channels responding to ligands not so far known to act on C. elegans cys-loop LGICs.

Cys-loop LGIC superfamilies of other nematode species
Genome sequencing projects are allowing the comparison of cys-loop LGIC superfamilies from different nematode species. For instance, Caenorhabditis briggsae, which is seemingly identical to C. elegans (Gupta et al. 2007;Stein et al. 2003), also has an extensive cys-loop LGIC gene superfamily although there are a few exceptions where there appears to be gene expansions in C. elegans, mainly within the nAChR orphan group. For example, lgc-23, lgc-24 and lgc-28 in C. elegans ) may be paralogues resulting from gene duplication in the elegans species lineage since only one homologue (Cbrlgc-28) is observed in C. briggsae. In line with this, lgc-23 and lgc-24 are situated close together within 6 kb in the C. elegans genome which may reflect a recent gene duplication event. Likewise, lgc-11 and lgc-12, as well as lgc-25 and lgc-26 have single homologues in C. briggsae, Cbr-lgc-11 and Cbr-lgc-26 respectively. Six C. elegans cys-loop LGIC subunits, lgc-13 to lgc-18, are tightly clustered within 20 kb of the T01H10 cosmid whereas in C. briggsae there are seven similar genes (CBG01219, CBG01221, CBG02706, CBG16212, CBG16226, CBG16227 and CBG16229), all of which are not clustered together. Outside of the nAChR orphan group, avr-15 and glc-1 may be paralogues arising in the elegans lineage as there is only a single homologue (CBG06688) in C. briggsae.
The cys-loop LGIC gene families of the parasitic nematodes, Brugia malyai and Trichinella spiralis, are much smaller than those of C. elegans and C. briggsae consisting of 30 and 19 subunits respectively (Williamson et al. 2007). The cys-loop LGIC subunits of both parasites have clear orthologous relationships with those of C. elegans with one exception, a nAChR subunit (ACR-26) which so far appears particular to B. malayi. As discussed by Williamson and colleagues, the striking difference in the cysloop LGIC gene family sizes may reflect a free-living versus a parasitic lifestyle. Thus, the large complement of subunits in C. elegans and C. briggsae may be required for a small nervous system to be able to respond to many environmental cues which may be more restricted in the specific environmental niche occupied by parasitic nematodes. This shows that the extensive cys-loop LGIC gene superfamily observed for C. elegans is not a feature of all nematodes. Wotring VE, Weiss DS (2008)