The cys-loop ligand-gated ion channel gene superfamily of the red flour beetle, Tribolium castaneum
© Jones and Sattelle; licensee BioMed Central Ltd. 2007
Received: 11 April 2007
Accepted: 19 September 2007
Published: 19 September 2007
Members of the cys-loop ligand-gated ion channel (cys-loop LGIC) superfamily mediate chemical neurotransmission and are studied extensively as potential targets of drugs used to treat neurological disorders such as Alzheimer's disease. Insect cys-loop LGICs are also of interest as they are targets of highly successful insecticides. The red flour beetle, Tribolium castaneum, is a major pest of stored agricultural products and is also an important model organism for studying development.
As part of the T. castaneum genome sequencing effort, we have characterized the beetle cys-loop LGIC superfamily which is the third insect superfamily to be described after those of Drosophila melanogaster and Apis mellifera, and also the largest consisting of 24 genes. As with Drosophila and Apis, Tribolium possesses ion channels gated by acetylcholine, γ-amino butyric acid (GABA), glutamate and histamine as well as orthologs of the Drosophila pH-sensitive chloride channel subunit (pHCl), CG8916 and CG12344. Similar to Drosophila and Apis, Tribolium cys-loop LGIC diversity is broadened by alternative splicing although the beetle orthologs of RDL and GluCl possess more variants of exon 3. Also, RNA A-to-I editing was observed in two Tribolium nicotinic acetylcholine receptor subunits, Tcasα6 and Tcasβ1. Editing in Tcasα6 is evolutionarily conserved with D. melanogaster, A. mellifera and Heliothis virescens, whereas Tcasβ1 is edited at a site so far only observed in the beetle.
Our findings reveal that in diverse insect species the cys-loop LGIC superfamily has remained compact with only minor changes in gene numbers. However, alternative splicing, RNA editing and the presence of divergent subunits broadens the cys-loop LGIC proteome and generates species-specific receptor isoforms. These findings on Tribolium castaneum enhance our understanding of cys-loop LGIC functional genomics and provide a useful basis for the development of improved insecticides that target an important agricultural pest.
In insects, members of the cys-loop ligand-gated ion channel (cys-loop LGIC) superfamily mediate both fast excitatory and inhibitory synaptic transmission in the nervous system. The superfamily includes cation permeable nicotinic acetylcholine receptors (nAChRs) [1, 2], γ-amino butyric acid (GABA)-gated anion channels , glutamate-gated chloride channels (GluCls)  and histamine-gated chloride channels (HisCls) [5, 6]. Studies of Drosophila melanogaster and Apis mellifera have shown that cys-loop LGICs mediate important aspects of behaviour such as escape response , learning and memory [8–12]. Members of the cys-loop superfamily of ionotropic receptors are also of considerable interest as they are targets of widely used insecticides . For example, nAChRs are targets of neonicotinoids [14, 15], a class of insect control chemicals which include imidacloprid with worldwide annual sales of approximately one billion US dollars . Also, GABA receptors, GluCls and HisCls are targets of fipronil and avermectins [17, 18].
The red flour beetle, Tribolium castaneum, is a highly sophisticated genetic and developmental model organism , and is a major global pest of stored agricultural products. In the USA alone, T. castaneum contributes to over one billion US dollars worth of damage to wheat and corn every year . Tribolium has proven highly adaptable, developing resistance to a wide range of insecticides raising the need for the development of improved and novel control agents. To date, partial DNA sequences of only a few T. castaneum cys-loop LGIC subunits have been reported. These include a HisCl , transcript variants of a nAChR subunit orthologous to Drosophila melanogaster Dα6 (submitted to NCBI by Jin and colleagues, Acession Numbers EF127806–EF127810) and a GABA-gated ion channel which is an ortholog of Drosophila RDL . Interestingly, a mutation changing alanine 302 to serine in the Tribolium GABA receptor is associated with cyclodiene resistance [23, 24]. Indeed, the same mutation is found in RDL from cyclodiene resistant strains of a wide range of insect species [25, 26]. Characterizing the full complement of Tribolium cys-loop LGIC subunits represents a critical step in identifying key components of the beetle nervous system as well pinpointing particular insecticide targets. Here we have used sequence data from the T. castaneum genome project  to provide the first description of a complete cys-loop LGIC gene superfamily from an invertebrate pest species.
The T. castaneum cys-loop LGIC superfamily consists of 24 subunit members
Percentage identity/similarity between T. castaneum and A. mellifera nAChR subunit protein sequences. Proposed orthologs are shown in bold.
Percentage identity/similarity between T. castaneum and D. melanogaster non-nAChR subunit protein sequences. Proposed orthologs are shown in bold.
Tribolium nicotinic acetylcholine receptor subunits
We identified 12 candidate nAChR subunit-encoding genes in the T. castaneum genome. This is the largest insect nAChR gene family so far described as those of D. melanogaster, Anopheles gambiae and A. mellifera consist of 10, 10 and 11 subunits respectively [1, 32, 33]. Eleven of the Tribolium nAChR subunits possess the two adjacent cysteine residues in loop C (Fig. 2) which are important for acetylcholine (ACh) binding , defining them as α subunits. The remaining subunit was designated β since it lacks the vicinal cysteines.
Tribolium possesses core groups of nAChR subunits that are highly conserved between different insect species . Thus, subunit equivalents of Dα1–7, Dβ1 and Dβ2 are evident in the beetle genome (Fig. 5). As with Anopheles, Apis and several other insects , the Tribolium ortholog of Dβ2 is of the α type (Tcasα8). Interestingly, the beetle possesses an additional α subunit (Tcasα11) in the Dβ2 group, indicating a gene duplication in the Tribolium lineage. In line with this, both Tcasα8 and Tcasα11 genes are tightly clustered together within 8 kb of each other in the beetle genome suggesting that both subunits arose from a recent duplication event from a common gene. Both A. gambiae and D. melanogaster possess three subunits (Dα5, Dα6 and Dα7) that are very similar and show substantial homology to the vertebrate α7 nAChR subunit [32, 37]. In T. castaneum, two orthologs of these subunits (Tcasα6 and Tcasα7) are also similar to human α7, sharing 46% and 45% identity respectively at the protein level. The third subunit, Tcasα5, when compared with vertebrate nAChR subunits, is most similar to α7 but shares only 25% identity and, along with its Apis ortholog (Amelα5), departs strongly from Dα5 (Fig. 5). Tribolium nAChR subunits outside of the Dα5–7 group show 19%–39% identity with vertebrate subunits. As is the case for Dα1, Dα2, Dα3, Dα4, Dβ2 and their Anopheles and Apis orthologs, the corresponding Tribolium subunits (Tcasα1–α4, Tcasα8 and Tcasα11) have an insertion in loop F (Fig. 2), which may contribute to imidacloprid interactions . The Dα1, Dα2 and Dβ2 genes, as well as their Anopheles orthologs, Agamα1, Agamα2 and Agamα8, are similarly arranged and tightly clustered within 200 kb and 220 kb respectively . Immunohistochemical and coimmunoprecipitation studies show that Dα1, Dα2 and Dβ2 are integral components of certain nAChRs subtypes, leading to the suggestion that clustering may facilitate coordinated expression and co-assembly of the nAChR subunits . In Apis, Amelα1 and Amelα2 are clustered but are separated from the honey bee Dβ2 ortholog, Amelα8 , while in the Tribolium genome all equivalent beetle subunits, Tcasα1, Tcasα2 and Tcasα8, are located on different linkage groups (Table 1). The separation of these genes may thus result in diversification of receptor expression and coassembly. Tribolium does however, show clustering of Tcasα7 and Tcasβ1 (both genes lie within 8 kb of each other) which is conserved in the genomes of Anopheles and Apis but not Drosophila .
Analysis of D. melanogaster, A. gambiae and A. mellifera nAChR gene families has shown that each insect possesses at least one divergent subunit that shares relatively low sequence identity with other nAChR subunits . The four insect nAChR gene families described so far each contain a different complement of divergent subunits. Thus, Drosophila and Anopheles have one divergent subunit each but are of the β and α types respectively [32, 40], while Apis and Tribolium each have two subunits which are α and β in the honey bee and are both α subunits in the beetle  (Fig. 5). One of the Tribolium divergent subunits, Tcasα9, possesses an atypical FxCC amino acid motif instead of YxCC found in loop C of all other insect α nAChR subunits characterised to date (Fig. 2). The nematode, Caenorhabditis elegans, is the only other organism known to possess nAChR subunits with the FxCC motif . Since site-directed mutagenesis has shown that a Tyr to Phe substitution in the heterologously expressed vertebrate α7 subunit results in a tenfold lower affinity for ACh , Tcasα9 may have unusual ligand binding properties. In addition, Tcasα9 lacks the GEK motif characteristic of nAChR subunits which precedes TM2 (Fig. 2) and plays an important role in ion permeation and selectivity . Most notably, the absence of the highly conserved glutamic acid residue may have given rise to a receptor with distinct ion channel properties since a substitution of the equivalent glutamic acid residue in the vertebrate α7 nAChR abolishes permeability to calcium ions but not monovalent cations .
Tribolium GABA-gated ion channels
The Tribolium genome contains orthologs of the three known D. melanogaster GABA-gated ion channels, RDL, GRD and LCCH3  (Fig. 5 and Table 2). RDL possesses a PAR sequence before TM2 which is characteristic of cys-loop ligand-gated anion channels  while GRD and LCCH3 lack this sequence. This may underlie the findings that RDL forms homomeric anion channels  whereas GRD and LCCH3 form heteromultimeric cation channels when expressed in Xenopus laevis oocytes . The PAR motif is also present in Tcas_RDL and is absent in Tcas_GRD and Tcas_LCCH3 (Fig. 2) indicating that Tribolium may possess both ligand-gated anion and cation channels. However, whereas electrophysiology has clearly shown that GABA induces inhibitory chloride channels in insects, it remains to be established whether GABA-gated cation channels function in vivo .
As is the case for Drosophila RDL, exons 3 and 6 are alternatively spliced in Tcas_RDL (Fig. 6C) although, whereas Drosophila and Apis have two alternatives for exon 3 [29, 56], Tribolium has three. To accommodate this extra alternative exon and anticipate the possibility of additional exons in other insect species, we have revised the nomenclature of RDL splice variants. Thus, RDL exons "a" and "b" are denoted exon 3a and exon 3b respectively, while exons "c" and "d" are now called exon 6a and exon 6b. We have designated the extra exon in Tcas_RDL as exon 3c. Exons 3a and 3b are completely conserved in Drosophila RDL , Amel_RDL  and Tcas_RDL, while exon 3c of Tribolium introduces two novel amino acid residues near to loop D (Fig. 6C). RT-PCR [see Additional file 3 for primers used] show that all six possible combinations of alternate exons for Tcas_RDL are transcribed. Since studies of Drosophila RDL demonstrate that alternative splicing alters affinity for GABA , the extra choice of exon in Tcas_RDL may give rise to a receptor with extended functional range when compared to orthologs in other insects. The most variation seen between a Tcas_RDL alternate exon and its equivalent in Drosophila is in exon 6a which differ by eight residues, most of which are located in the vicinity of loop F (Fig. 6C). This may give rise to receptor variants that have distinct ligand binding characteristics in the two insect species.
Tribolium glutamate and histamine-gated chloride channels
As with D. melanogaster and A. mellifera, T. castaneum has one known glutamate-gated chloride channel (Tcas_GluCl) and two histamine-gated chloride channels (Tcas_HisCl1 and Tcas_HisCl2). Consistent with their putative function as anion channels, Tcas_GluCl, Tcas_HisCl1 and Tcas_HisCl2 all have the PAR motif preceding TM2 (Fig. 2) . Out of the ligand-gated anion channels, GluCl and HisCl2 are the most highly conserved between the fruit fly and beetle, sharing 83% and 79% identity respectively (Table 2).
Exon 3 of Tcas_GluCl is alternatively spliced as is the case with Drosophila GluCl and Apis GluCl [33, 50]. However, whereas the fruit fly and honey bee GluCls each have two alternatives for exon 3, Tribolium has three (Fig. 6D). Sequence chromatograms of Tcas_GluCl RT-PCR products [using primers shown in Additional file 1] showed mixed peaks corresponding to exon 3 indicating multiple variants are transcribed. In order to maintain consistency with other cys-loop LGIC subunits, we have revised the nomenclature of GluCl alternative exons so that Modules 1 and 2 correspond to exons 3a and 3b respectively. Similar to RDL, alternative splicing of GluCl alters amino acid residues in the vicinity of loop D which may introduce variation in ligand-binding properties. With the extra alternative exon (exon 3c), Tcas_GluCl may have a wider range in receptor functional properties when compared to GluCls of other insect species.
Tribolium pHCl and uncharacterized cys-loop LGIC subunits
The pH sensitive chloride channel (pHCl) first identified in Drosophila  is well conserved in T. castaneum, showing 68% identity (Table 2). The Drosophila pHCl has several splice variants, of which Variant 3 introduces an insertion in the intracellular region between TM3 and TM4. We detected a similar variant in Tcas_pHCl transcripts (Fig. 8B) where the peptide insertion differs by only one amino acid residue from that of Drosophila pHCl  and is completely identical to the equivalent insertion in Amel_pHCl . However, unlike Drosophila and Apis, the Variant 3 insertion does not introduce a potential protein kinase C phosphorylation site in Tcas_pHCl. This may lead to the beetle pHCl having distinct characteristics since phosphorylation of the large intracellular region modulates receptor assembly and function [60, 61], Interestingly, we detected a second insertion at the same site which has not been observed in Drosophila or Apis. Denoted Tcas_pHCl Variant 3a (Fig. 8B), this insertion introduces a putative casein kinase II phosphorylation site. Tcas_pHCl also has an insertion in loop C (Variant 4, Fig. 8C) which is likely to impact on ligand binding. The equivalent of this insertion has not been observed in Drosophila but has been detected in Apis with the insertion sequence and length differing considerably between beetle and honey bee (Fig. 8C).
Five Drosophila cys-loop LGIC subunits have yet to be functionally characterized. These are CG6927, CG7589, CG8916, CG11340 and CG12344. CG8916 and CG12344 appear to be closely related with GRD and HisCls respectively while CG6927, CG7589 and CG11340 forms a distinct subfamily of cys-loop LGIC subunits which, following on from a recent study , we have denoted Insect Group 1 of cys-loop LGIC subunits (Fig. 5). Tribolium also possesses five uncharacterized subunits. Two of these, Tcas_8916 and Tcas_12344, are candidate orthologs of CG8916 and CG12344 being notably similar to their Drosophila counterparts, both sharing 54% identity (Table 2). The remaining three subunits belong to Insect Group I (Fig. 5) and have been denoted Tcas_CLGC1, Tcas_CLGC2 and Tcas_CLGC3, standing for c ys-loop l igand g ated ion c hannel. Despite sharing highest sequence identities with Drosophila Insect Group I subunits (Table 2), their orthologous relationships are uncertain since, for example, Tcas_CLGC2 shows 36% identity with both CG6927 and CG7589 (Table 2). A. mellifera possesses only one subunit, Amel_6927, in Insect Group 1 , indicating that gene duplication occurred after the emergence of the Hymenoptera to give rise to the three subunits present in both Tribolium (Coleoptera) and Drosophila (Diptera). This is in line with recent findings suggesting that the Hymenoptera are basal to the Coleoptera [63, 64]. Consistent with gene duplication occurring within Insect Group 1 just before the emergence of Coleoptera, Tcas_CLGC1, Tcas_CLGC2 and Tcas_CLGC3 are tightly clustered together in the Tribolium genome, being located within 10 Kb in linkage group 8. In the more evolutionarily advanced D. melanogaster , the three genes are separated with CG7589 and CG11340 being respectively located on the left and right arms of chromosome 3 and CG6927 being present on chromosome X .
Insect genome sequencing projects have allowed the identification and comparison of gene superfamilies from diverse species. As part of the Tribolium sequencing project , we have described the beetle cys-loop LGIC superfamily which encodes for receptors that play major roles in the nervous system and are also targets of highly successful insecticides. This is the first complete cys-loop LGIC superfamily to be described from a Coleoptera and an agricultural pest species, and is the third to be reported after those of the Dipteran D. melanogaster [29, 62, 66] and the Hymenopteran A. mellifera . In the three insect species, which represents over 300 million years of evolution , the cys-loop LGIC superfamily has remained compact with only minor changes in gene numbers. However, alternative splicing and RNA A-to-I editing have considerably increased receptor diversity, effectively introducing changes in functionally significant and highly conserved regions to generate subunit isoforms particular to certain insect species. Also, it is becoming apparent that each insect possesses a distinct complement of highly divergent nAChR subunits whose sequences do not reflect a high degree of evolutionary constraint and thus may play diverse roles in different species. In addition, a group of cys-loop LGIC subunits that appear to be particular to insects, which was noted as D. melanogaster Group 1 , and denoted here as Insect Group 1 to accommodate sequences from other species (Fig. 5), may represent more recent members of the superfamily since only one subunit is found in Apis and three are present after the emergence of Coleoptera. It will be of interest to determine the functional role played by Insect Group 1 subunits and the ligands to which they respond.
The species-specific diversification arising from alternative splicing, RNA editing and divergent subunits, as well as insect-specific subunits, represents promising receptor differences to target for the future rational design of insecticides that control pest species while sparing beneficial insects. The use of heterologous expression systems such as Xenopus laevis oocytes has allowed the functional characterisation of several Drosophila cys-loop LGICs such as RDL , GRD and LCCH3 , GluCl , HisCl1 and HisCl2 [5, 6] and pHCl . Similar studies of heterologously expressed ion channels from other insect species including T. castaneum, in combination with the use of three-dimensional models based on the crystal structure of molluscan acetylcholine binding proteins [1, 68–70], will likely prove useful in the search for novel compounds that show selectivity for receptors of certain insect species as well as in determining the mechanisms of insecticide interaction with cys-loop LGICs. For insect nAChRs, functional expression in heterologous systems has so far proven elusive , although low levels of receptor activity have been observed for the locust Schistocerca gregaria αL1 subunit expressed in Xenopus oocytes . Nevertheless, Drosophila nAChR α subunits can form robust functional channels when coexpressed with a vertebrate β2 subunit  and studies on such hybrid receptors have provided insights into the selectivity of neonicotinoids for insect nAChRs over those of vertebrates, regions of subunit proteins involved in imidacloprid interactions and the actions of different neonicotinoids . These studies have highlighted Dα1 and Dα2 as being sensitive to imidacloprid. Also, Mpα2, which is the aphid Myzus persicae ortholog of Dα1, shows high levels of imidacloprid binding when coexpressed with the rat β2 subunit . The functional expression of nAChRs with insect β subunits has yet to be achieved but it is worth noting that the ortholog of Dβ2 is an α subunit (e.g. Tcasα8 in Tribolium) in all insects so far studied outside the Drosophila genus . Since members of the Dβ2 group are closely related to those of the Dα1 and Dα2 groups (Fig. 5) and share an insertion in loop F which may contribute to imidacloprid sensitivity , it would be of interest to determine whether Tcasα8 is sensitive to neonicotinoids. If this is the case, it would also be worth studying Tcasα11, which appears to be a product of gene duplication of Tcasα8, and assessing how both subunits contribute to neonicotinoid sensitivity particularly in light of the finding that gene duplication has given rise to insecticide resistance at another synaptic target site, acetylcholinesterase .
Parental RNAi, where RNA interference arising from double-stranded RNA introduced into the mother also spreads to the offspring, is highly efficient in Tribolium [27, 76]. The combination of genome information and the use of RNAi in C. elegans has proved considerably instructive in determining roles played by genes [77, 78]. Thus, the beetle provides a powerful tool for studying gene function in an insect pest species. For example, RNAi could be used to elucidate the roles played by Tribolium cys-loop LGIC subunits in various aspects of development, behaviour and response to insecticides. Recently, it has been shown that a Drosophila Dα6 knockout mutant is highly resistant to spinosad . Studies have shown that T. castaneum is susceptible to spinosad, although to a lesser degree than other insect pests of stored wheat . It will be of interest to determine the effect of silencing Tcasα6 on the beetle's susceptibility to spinosad and perhaps validate in a pest species findings based on the Drosophila genetic model organism.
Using information from the Tribolium castaneum genome sequencing project, we report, for the first time from an invertebrate pest species, a complete cys-loop LGIC superfamily, which encodes for receptors that play important roles in the nervous system as well as for targets of widely-used insecticides. The present study enhances our understanding of the functional genomics of the insect cys-loop LGIC superfamily. Our findings reveal an emerging consensus that in over 300 million years of insect evolution, the cys-loop LGIC superfamily has remained compact with only minor changes in gene numbers. However, alternative splicing, RNA editing and the presence of divergent subunits broadens the cys-loop LGIC proteome and generates species-specific receptor isoforms. Thus, the paper provides several new insights into the molecular diversity of cys-loop LGICs between different organisms and provides an important foundation for associating particular cys-loop LGIC subtypes with development as well as for the generation of improved insecticides that target the red flour beetle.
Identification of cys-loop LGIC subunits in the T. castaneum genome
To identify putative cys-loop LGIC subunits, we screened the T. castaneum genome (assembly version 2.0)  with cDNA sequences of every member of the D. melanogaster and A. mellifera cys-loop LGIC superfamilies using the tBLASTn algorithm . Candidate beetle cys-LGIC subunits were identified based on their considerable sequence homology with previously characterized subunits (sequences with lowest similarity had E Value 8e-19), particularly in the N-terminal ligand-binding domain and the four transmembrane regions. The highly variable N-terminal signal peptides, which are a feature of cys-loop LGIC subunits, were identified in the GLEAN consensus set of predicted genes . RT-PCRs were performed [see Additional file 1 for primers used] to verify and correct the open-reading frame sequences of each subunit.
Reverse transcription and polymerase chain reaction
Total RNA was extracted from 15 Tribolium castaneum adult beetles (Georgia GA2 strain) homogenized in Trizol (Invitrogen) using the RNeasy Mini Kit (Qiagen) and first strand cDNA was synthesized from 1 μg total RNA using Superscript™ III First-Strand Synthesis Super Mix (Invitrogen). Nested RT-PCR reactions were performed to detect transcripts of beetle cys-loop LGIC subunits as well as to detect transcript variants arising from alternative splicing. Primer pairs [see Additional files 1 and 3 for primer sequences] recognising different exons were used to allow identification of cDNA-specific products. The PCR reactions were performed in a total volume of 50 μl composed of Taq polymerase and 1 × PCR buffer (Sigma), 0.2 mM dNTP mix (Roche), 0.4 μM each primer and 2 μl first strand cDNA template. The PCR reaction conditions were 35 cycles of: 95°C for 30 s, 55°C for 30 s, 72°C for 90 s. The first PCR was used at a final dilution of 1 in 5000 as template for the second, nested PCR reaction. For Tcas_RDL and Tcas_GluCl, a 1 in 500 dilution was used to amplify enough DNA for sequencing. DNA sequence chromatograms for each cys-loop LGIC subunit were analyzed using Chromas 2 (Technelysium Pty Ltd) to detect single nucleotide polymorphisms (SNPs) or RNA editing sites as shown by mixed signal peaks. No SNPs were observed and the putative RNA editing sites detected in Tcasα6 and Tcasβ1 were verified by amplifying and sequencing genomic DNA present in the extracted total RNA, which was first treated with DNase-free RNase (Roche), using primers recognising intron DNA [see Additional file 4 for primers used]. Sequence chromatograms showing a defined region of mixed peaks indicated differential splicing. The corresponding cys-loop LGIC PCR products were cloned into the pGEM-T Easy vector (Promega) and between 10 to 20 transformants were sequenced to identify individual subunit isoforms. All PCR products were analyzed by electrophoresis in a TAE gel and then purified using the QIAquick Gel Extraction Kit (Qiagen) while subunits cloned in pGEM-T Easy were purified using the QIAprep Spin Miniprep Kit (Qiagen). Purified DNA was sequenced by the dye termination method at Cogenics . For analyzing RNA editing levels in Tcasα6 splice variants, RT-PCR was performed using forward primers specific to either exon 3a or exon 3b and reverse primers recognising one of the three alternatives for exon 8 [see Additional file 3 for primers used]. A nested PCR approach was adopted since two rounds of PCR were required to amplify enough variant containing exon 8c for visualization on an agarose gel. The proofreading Pfu Turbo DNA polymerase (Stratagene) was used in 2 × 30-cycle reactions on four independently made first-strand cDNAs. The sequence chromatograms of the amplification products were analyzed to give editing levels where the proportion edited = height of guanosine peak/(height of guanosine peak + height of adenosine peak). One-way and two-way ANOVA analyses were performed with Turkey's Multiple Comparison and Bonferroni tests respectively using Graphpad Prism 4 . Products amplified by either one or two PCR reactions gave similar editing levels.
The multiple protein sequence alignment was constructed with ClustalX  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 . The protein alignment was viewed using GeneDoc . Identity values between subunit sequences were calculated using the GeneDoc program. The neighbour-joining method  and bootstrap resampling , available with the ClustalX program, were used to construct a phylogenetic tree, which was then displayed using the TreeView application . Signal peptide cleavage sites were predicted using the SignalP 3.0 server  and membrane-spanning regions were predicted using the TMpred program . The PROSITE database  was used to identify potential phosphorylation sites.
cys-loop ligand-gated ion channel
γ-amino butyric acid
glutamate-gated chloride channel
GABAA and glycine receptor-like subunit of Drosophila
histaminegated chloride channel
ligand-gated chloride channel homolog 3
ligand-gated ion channel
nicotinic acetylcholine receptor
resistant to Dieldrin
We are indebted to the T. castaneum Genome Project (Human Genome Sequencing Center), which provided the starting point for this study. We thank Sue Brown for kindly providing us with Tribolium castaneum beetles. This work was funded by the Medical Research Council (MRC), Dupont and the Wellcome Trust Integrative Physiology Initiative in Ion Channels (OXION).
- Sattelle DB, Jones AK, Sattelle BM, Matsuda K, Reenan R, Biggin PC: Edit, cut and paste in the nicotinic acetylcholine receptor gene family of Drosophila melanogaster. Bioessays. 2005, 27 (4): 366-376. 10.1002/bies.20207.PubMedView Article
- Tomizawa M, Casida JE: Structure and diversity of insect nicotinic acetylcholine receptors. Pest Manag Sci. 2001, 57 (10): 914-922. 10.1002/ps.349.PubMedView Article
- Buckingham SD, Biggin PC, Sattelle BM, Brown LA, Sattelle DB: Insect GABA receptors: splicing, editing, and targeting by antiparasitics and insecticides. Mol Pharmacol. 2005, 68 (4): 942-951. 10.1124/mol.105.015313.PubMedView Article
- Vassilatis DK, Elliston KO, Paress PS, Hamelin M, Arena JP, Schaeffer JM, Van der Ploeg LH, Cully DF: Evolutionary relationship of the ligand-gated ion channels and the avermectin-sensitive, glutamate-gated chloride channels. J Mol Evol. 1997, 44 (5): 501-508. 10.1007/PL00006174.PubMedView Article
- Gisselmann G, Pusch H, Hovemann BT, Hatt H: Two cDNAs coding for histamine-gated ion channels in D. melanogaster. Nat Neurosci. 2002, 5 (1): 11-12. 10.1038/nn787.PubMedView Article
- Zheng Y, Hirschberg B, Yuan J, Wang AP, Hunt DC, Ludmerer SW, Schmatz DM, Cully DF: Identification of two novel Drosophila melanogaster histamine-gated chloride channel subunits expressed in the eye. J Biol Chem. 2002, 277 (3): 2000-2005. 10.1074/jbc.M107635200.PubMedView Article
- Fayyazuddin A, Zaheer MA, Hiesinger PR, Bellen HJ: The nicotinic acetylcholine receptor Dalpha7 is required for an escape behavior in Drosophila. PLoS Biol. 2006, 4 (3): e63-10.1371/journal.pbio.0040063.PubMed CentralPubMedView Article
- Gauthier M, Dacher M, Thany SH, Niggebrugge C, Deglise P, Kljucevic P, Armengaud C, Grunewald B: Involvement of alpha-bungarotoxin-sensitive nicotinic receptors in long-term memory formation in the honeybee (Apis mellifera). Neurobiol Learn Mem. 2006, 86 (2): 164-174. 10.1016/j.nlm.2006.02.003.PubMedView Article
- Thany SH, Gauthier M: Nicotine injected into the antennal lobes induces a rapid modulation of sucrose threshold and improves short-term memory in the honeybee Apis mellifera. Brain Res. 2005, 1039 (1-2): 216-219. 10.1016/j.brainres.2005.01.056.PubMedView Article
- Dacher M, Lagarrigue A, Gauthier M: Antennal tactile learning in the honeybee: effect of nicotinic antagonists on memory dynamics. Neuroscience. 2005, 130 (1): 37-50. 10.1016/j.neuroscience.2004.09.006.PubMedView Article
- Hosler JS, Buxton KL, Smith BH: Impairment of olfactory discrimination by blockade of GABA and nitric oxide activity in the honey bee antennal lobes. Behav Neurosci. 2000, 114 (3): 514-525. 10.1037/0735-7044.114.3.514.PubMedView Article
- El Hassani AK, Dacher M, Gauthier M, Armengaud C: Effects of sublethal doses of fipronil on the behavior of the honeybee (Apis mellifera). Pharmacol Biochem Behav. 2005, 82 (1): 30-39. 10.1016/j.pbb.2005.07.008.PubMedView Article
- Raymond-Delpech V, Matsuda K, Sattelle BM, Rauh JJ, Sattelle DB: Ion channels: molecular targets of neuroactive insecticides. Invert Neurosci. 2005, 5 (3-4): 119-133. 10.1007/s10158-005-0004-9.PubMedView Article
- Matsuda K, Buckingham SD, Kleier D, Rauh JJ, Grauso M, Sattelle DB: Neonicotinoids: insecticides acting on insect nicotinic acetylcholine receptors. Trends Pharmacol Sci. 2001, 22 (11): 573-580. 10.1016/S0165-6147(00)01820-4.PubMedView Article
- Millar NS, Denholm I: Nicotinic acetylcholine receptors: targets for commercially important insecticides. Invert Neurosci. 2007
- Tomizawa M, Casida JE: Neonicotinoid insecticide toxicology: mechanisms of selective action. Annu Rev Pharmacol Toxicol. 2005, 45: 247-268. 10.1146/annurev.pharmtox.45.120403.095930.PubMedView Article
- Bloomquist JR: Chloride channels as tools for developing selective insecticides. Arch Insect Biochem Physiol. 2003, 54 (4): 145-156. 10.1002/arch.10112.PubMedView Article
- Iovchev M, Kodrov P, Wolstenholme AJ, Pak WL, Semenov EP: Altered drug resistance and recovery from paralysis in Drosophila melanogaster with a deficient histamine-gated chloride channel. J Neurogenet. 2002, 16 (4): 249-261. 10.1080/01677060216293.PubMedView Article
- Brown SJ, Denell RE, Beeman RW: Beetling around the genome. Genet Res. 2003, 82 (3): 155-161. 10.1017/S0016672303006451.PubMedView Article
- Throne JE, Hallman GJ, Johnson JA, Follett PA: Post-harvest entomology research in the United States Department of Agriculture-Agricultural Research Service. Pest Manag Sci. 2003, 59 (6-7): 619-628. 10.1002/ps.690.PubMedView Article
- Iovchev M, Boutanaev A, Ivanov I, Wolstenholme A, Nurminsky D, Semenov E: Phylogenetic shadowing of a histamine-gated chloride channel involved in insect vision. Insect Biochem Mol Biol. 2006, 36 (1): 10-17. 10.1016/j.ibmb.2005.09.002.PubMedView Article
- Thompson M, Steichen JC, ffrench-Constant RH: Conservation of cyclodiene insecticide resistance-associated mutations in insects. Insect Mol Biol. 1993, 2 (3): 149-154.PubMedView Article
- Andreev D, Rocheleau T, Phillips TW, Beeman RW, ffrench-Constant RH: A PCR diagnostic for cyclodiene insecticide resistance in the red flour beetle Tribolium castaneum. Pestic Sci. 1994, 41: 345-349. 10.1002/ps.2780410410.View Article
- Miyazaki M, Matsumura F, Beeman RW: DNA sequence and site of mutation of the GABA receptor of cyclodiene-resistant red flour beetle, Tribolium castaneum. Comp Biochem Physiol B Biochem Mol Biol. 1995, 111 (3): 399-406. 10.1016/0305-0491(95)00007-U.PubMedView Article
- Ffrench-Constant RH, Anthony N, Aronstein K, Rocheleau T, Stilwell G: Cyclodiene insecticide resistance: from molecular to population genetics. Annu Rev Entomol. 2000, 45: 449-466. 10.1146/annurev.ento.45.1.449.PubMedView Article
- Ffrench-Constant RH, Rocheleau TA, Steichen JC, Chalmers AE: A point mutation in a Drosophila GABA receptor confers insecticide resistance. Nature. 1993, 363 (6428): 449-451. 10.1038/363449a0.PubMedView Article
- Tribolium Genome Sequencing Consortium: The genome sequence of the first Coleoptera, Tribolium castaneum: A model of insect development with agricultural insights. In Press. 2007
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView Article
- Jones AK, Sattelle DB: The cys-loop ligand-gated ion channel superfamily of the honeybee, Apis mellifera. Invert Neurosci. 2006, 6 (3): 123-132. 10.1007/s10158-006-0026-y.PubMedView Article
- Sine SM, Engel AG: Recent advances in Cys-loop receptor structure and function. Nature. 2006, 440 (7083): 448-455. 10.1038/nature04708.PubMedView Article
- Corringer PJ, Le Novere N, Changeux JP: Nicotinic receptors at the amino acid level. Annu Rev Pharmacol Toxicol. 2000, 40: 431-458. 10.1146/annurev.pharmtox.40.1.431.PubMedView Article
- Jones AK, Grauso M, Sattelle DB: The nicotinic acetylcholine receptor gene family of the malaria mosquito, Anopheles gambiae. Genomics. 2005, 85 (2): 176-187. 10.1016/j.ygeno.2004.09.001.PubMedView Article
- Jones AK, Raymond-Delpech V, Thany SH, Gauthier M, Sattelle DB: The nicotinic acetylcholine receptor gene family of the honey bee, Apis mellifera. Genome Res. 2006, 16 (11): 1422-1430. 10.1101/gr.4549206.PubMed CentralPubMedView Article
- Kao PN, Karlin A: Acetylcholine receptor binding site contains a disulfide cross-link between adjacent half-cystinyl residues. J Biol Chem. 1986, 261 (18): 8085-8088.PubMed
- Jones AK, Brown LA, Sattelle DB: Insect nicotinic acetylcholine receptor gene families: from genetic model organism to vector, pest and beneficial species. Invert Neurosci. 2007, 7 (1): 67-73. 10.1007/s10158-006-0039-6.PubMedView Article
- Jones AK, Marshall J, Blake AD, Buckingham SD, Darlison MG, Sattelle DB: Sgbeta1, a novel locust (Schistocerca gregaria) non-alpha nicotinic acetylcholine receptor-like subunit with homology to the Drosophila melanogaster Dbeta1 subunit. Invert Neurosci. 2005, 5 (3-4): 147-155. 10.1007/s10158-005-0007-6.PubMedView Article
- Grauso M, Reenan RA, Culetto E, Sattelle DB: Novel Putative Nicotinic Acetylcholine Receptor Subunit Genes, Dalpha5, Dalpha6 and Dalpha7, in Drosophila melanogaster Identify a New and Highly Conserved Target of Adenosine Deaminase Acting on RNA-Mediated A-to-I Pre-mRNA Editing. Genetics. 2002, 160 (4): 1519-1533.PubMed CentralPubMed
- Shimomura M, Yokota M, Matsuda K, Sattelle DB, Komai K: Roles of loop C and the loop B-C interval of the nicotinic receptor alpha subunit in its selective interactions with imidacloprid in insects. Neurosci Lett. 2004, 363 (3): 195-198. 10.1016/j.neulet.2003.12.115.PubMedView Article
- Chamaon K, Smalla KH, Thomas U, Gundelfinger ED: Nicotinic acetylcholine receptors of Drosophila: three subunits encoded by genomically linked genes can co-assemble into the same receptor complex. J Neurochem. 2002, 80 (1): 149-157. 10.1046/j.0022-3042.2001.00685.x.PubMedView Article
- Lansdell SJ, Millar NS: Dbeta3, an atypical nicotinic acetylcholine receptor subunit from Drosophila : molecular cloning, heterologous expression and coassembly. J Neurochem. 2002, 80 (6): 1009-1018. 10.1046/j.0022-3042.2002.00789.x.PubMedView Article
- Jones AK, Sattelle DB: Functional genomics of the nicotinic acetylcholine receptor gene family of the nematode, Caenorhabditis elegans. Bioessays. 2004, 26 (1): 39-49. 10.1002/bies.10377.PubMedView Article
- Galzi JL, Bertrand D, Devillers-Thiery A, Revah F, Bertrand S, Changeux JP: Functional significance of aromatic amino acids from three peptide loops of the alpha 7 neuronal nicotinic receptor site investigated by site-directed mutagenesis. FEBS Lett. 1991, 294 (3): 198-202. 10.1016/0014-5793(91)80668-S.PubMedView Article
- Jensen ML, Schousboe A, Ahring PK: Charge selectivity of the Cys-loop family of ligand-gated ion channels. J Neurochem. 2005, 92 (2): 217-225. 10.1111/j.1471-4159.2004.02883.x.PubMedView Article
- Bertrand D, Galzi JL, Devillers-Thiery A, Bertrand S, Changeux JP: Mutations at two distinct sites within the channel domain M2 alter calcium permeability of neuronal alpha 7 nicotinic receptor. Proc Natl Acad Sci U S A. 1993, 90 (15): 6971-6975. 10.1073/pnas.90.15.6971.PubMed CentralPubMedView Article
- Kondrashov FA, Koonin EV: Origin of alternative splicing by tandem exon duplication. Hum Mol Genet. 2001, 10 (23): 2661-2669. 10.1093/hmg/10.23.2661.PubMedView Article
- Lansdell SJ, Millar NS: Cloning and heterologous expression of Dalpha4, a Drosophila neuronal nicotinic acetylcholine receptor subunit: identification of an alternative exon influencing the efficiency of subunit assembly. Neuropharmacology. 2000, 39 (13): 2604-2614. 10.1016/S0028-3908(00)00111-8.PubMedView Article
- Green WN, Wanamaker CP: The role of the cystine loop in acetylcholine receptor assembly. J Biol Chem. 1997, 272 (33): 20945-20953. 10.1074/jbc.272.33.20945.PubMedView Article
- Imoto K, Busch C, Sakmann B, Mishina M, Konno T, Nakai J, Bujo H, Mori Y, Fukuda K, Numa S: Rings of negatively charged amino acids determine the acetylcholine receptor channel conductance. Nature. 1988, 335 (6191): 645-648. 10.1038/335645a0.PubMedView Article
- Hoopengardner B, Bhalla T, Staber C, Reenan R: Nervous system targets of RNA editing identified by comparative genomics. Science. 2003, 301 (5634): 832-836. 10.1126/science.1086763.PubMedView Article
- Semenov EP, Pak WL: Diversification of Drosophila chloride channel gene by multiple posttranscriptional mRNA modifications. J Neurochem. 1999, 72 (1): 66-72. 10.1046/j.1471-4159.1999.0720066.x.PubMedView Article
- Hoopengardner B: Adenosine-to-inosine RNA editing: perspectives and predictions. Mini Rev Med Chem. 2006, 6 (11): 1213-1216. 10.2174/138955706778742812.PubMedView Article
- Gehle VM, Walcott EC, Nishizaki T, Sumikawa K: N-glycosylation at the conserved sites ensures the expression of properly folded functional ACh receptors. Brain Res Mol Brain Res. 1997, 45 (2): 219-229. 10.1016/S0169-328X(96)00256-2.PubMedView Article
- Nishizaki T: N-glycosylation sites on the nicotinic ACh receptor subunits regulate receptor channel desensitization and conductance. Brain Res Mol Brain Res. 2003, 114 (2): 172-176. 10.1016/S0169-328X(03)00171-2.PubMedView Article
- Agrawal R, Stormo GD: Editing efficiency of a Drosophila gene correlates with a distant splice site selection. Rna. 2005, 11 (5): 563-566. 10.1261/rna.7280605.PubMed CentralPubMedView Article
- Gisselmann G, Plonka J, Pusch H, Hatt H: Drosophila melanogaster GRD and LCCH3 subunits form heteromultimeric GABA-gated cation channels. Br J Pharmacol. 2004, 142 (3): 409-413. 10.1038/sj.bjp.0705818.PubMed CentralPubMedView Article
- Ffrench-Constant RH, Rocheleau TA: Drosophila gamma-aminobutyric acid receptor gene Rdl shows extensive alternative splicing. J Neurochem. 1993, 60 (6): 2323-2326. 10.1111/j.1471-4159.1993.tb03523.x.PubMedView Article
- Hosie AM, Buckingham SD, Presnail JK, Sattelle DB: Alternative splicing of a Drosophila GABA receptor subunit gene identifies determinants of agonist potency. Neuroscience. 2001, 102 (3): 709-714. 10.1016/S0306-4522(00)00483-8.PubMedView Article
- Harvey RJ, Schmitt B, Hermans-Borgmeyer I, Gundelfinger ED, Betz H, Darlison MG: Sequence of a Drosophila ligand-gated ion-channel polypeptide with an unusual amino-terminal extracellular domain. J Neurochem. 1994, 62 (6): 2480-2483.PubMedView Article
- Schnizler K, Saeger B, Pfeffer C, Gerbaulet A, Ebbinghaus-Kintscher U, Methfessel C, Franken EM, Raming K, Wetzel CH, Saras A, Pusch H, Hatt H, Gisselmann G: A novel chloride channel in Drosophila melanogaster is inhibited by protons. J Biol Chem. 2005, 280 (16): 16254-16262. 10.1074/jbc.M411759200.PubMedView Article
- Cascio M: Modulating inhibitory ligand-gated ion channels. Aaps J. 2006, 8 (2): E353-61. 10.1208/aapsj080240.PubMed CentralPubMedView Article
- Thany SH, Lenaers G, Raymond-Delpech V, Sattelle DB, Lapied B: Exploring the pharmacological properties of insect nicotinic acetylcholine receptors. Trends Pharmacol Sci. 2007, 28 (1): 14-22. 10.1016/j.tips.2006.11.006.PubMedView Article
- Dent JA: Evidence for a diverse cys-loop ligand-gated ion channel superfamily in early bilateria. J Mol Evol. 2006, 62 (5): 523-535. 10.1007/s00239-005-0018-2.PubMedView Article
- Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443 (7114): 931-949. 10.1038/nature05260.View Article
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006, 16 (11): 1334-1338. 10.1101/gr.5204306.PubMed CentralPubMedView Article
- Flybase. [http://flybase.bio.indiana.edu/]
- Littleton JT, Ganetzky B: Ion channels and synaptic organization: analysis of the Drosophila genome. Neuron. 2000, 26 (1): 35-43. 10.1016/S0896-6273(00)81135-6.PubMedView Article
- Cully DF, Paress PS, Liu KK, Schaeffer JM, Arena JP: Identification of a Drosophila melanogaster glutamate-gated chloride channel sensitive to the antiparasitic agent avermectin. J Biol Chem. 1996, 271 (33): 20187-20191. 10.1074/jbc.271.33.20187.PubMedView Article
- Celie PH, Klaassen RV, van Rossum-Fikkert SE, van Elk R, van Nierop P, Smit AB, Sixma TK: Crystal structure of acetylcholine-binding protein from Bulinus truncatus reveals the conserved structural scaffold and sites of variation in nicotinic acetylcholine receptors. J Biol Chem. 2005, 280 (28): 26457-26466. 10.1074/jbc.M414476200.PubMedView Article
- Ernst M, Brauchart D, Boresch S, Sieghart W: Comparative modeling of GABA(A) receptors: limits, insights, future developments. Neuroscience. 2003, 119 (4): 933-943. 10.1016/S0306-4522(03)00288-4.PubMedView Article
- Smit AB, Brejc K, Syed N, Sixma TK: Structure and function of AChBP, homologue of the ligand-binding domain of the nicotinic acetylcholine receptor. Ann N Y Acad Sci. 2003, 998: 81-92. 10.1196/annals.1254.010.PubMedView Article
- Marshall J, Buckingham SD, Shingai R, Lunt GG, Goosey MW, Darlison MG, Sattelle DB, Barnard EA: Sequence and functional expression of a single alpha subunit of an insect nicotinic acetylcholine receptor. Embo J. 1990, 9 (13): 4391-4398.PubMed CentralPubMed
- Bertrand D, Ballivet M, Gomez M, Bertrand S, Phannavong B, Gundelfinger ED: Physiological properties of neuronal nicotinic receptors reconstituted from the vertebrate beta 2 subunit and Drosophila alpha subunits. Eur J Neurosci. 1994, 6 (5): 869-875. 10.1111/j.1460-9568.1994.tb00997.x.PubMedView Article
- Matsuda K, Shimomura M, Ihara M, Akamatsu M, Sattelle DB: Neonicotinoids show selective and diverse actions on their nicotinic receptor targets: electrophysiology, molecular biology, and receptor modeling studies. Biosci Biotechnol Biochem. 2005, 69 (8): 1442-1452. 10.1271/bbb.69.1442.PubMedView Article
- Tomizawa M, Millar NS, Casida JE: Pharmacological profiles of recombinant and native insect nicotinic acetylcholine receptors. Insect Biochem Mol Biol. 2005, 35 (12): 1347-1355. 10.1016/j.ibmb.2005.08.006.PubMedView Article
- Labbe P, Berthomieu A, Berticat C, Alout H, Raymond M, Lenormand T, Weill M: Independent Duplications of the Acetylcholinesterase Gene Conferring Insecticide Resistance in the Mosquito Culex pipiens. Mol Biol Evol. 2007
- Bucher G, Scholten J, Klingler M: Parental RNAi in Tribolium (Coleoptera). Curr Biol. 2002, 12 (3): R85-6. 10.1016/S0960-9822(02)00666-8.PubMedView Article
- Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC: Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans. Nature. 1998, 391 (6669): 806-811. 10.1038/35888.PubMedView Article
- Buckingham SD, Esmaeili B, Wood M, Sattelle DB: RNA interference: from model organisms towards therapy for neural and neuromuscular disorders. Hum Mol Genet. 2004, 13 Spec No 2: R275-88. 10.1093/hmg/ddh224.PubMedView Article
- Perry T, McKenzie JA, Batterham P: A Dalpha6 knockout strain of Drosophila melanogaster confers a high level of resistance to spinosad. Insect Biochem Mol Biol. 2007, 37 (2): 184-188. 10.1016/j.ibmb.2006.11.009.PubMedView Article
- Huang F, Subramanyam B, Toews MD: Susceptibility of laboratory and field strains of four stored-product insect species to spinosad. J Econ Entomol. 2004, 97 (6): 2154-2159.PubMedView Article
- Tribolium castaneum genome project. [http://www.hgsc.bcm.tmc.edu/projects/tribolium/]
- Cogenics. [http://www.cogenics.com/]
- Graphpad software. [http://www.graphpad.com]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralPubMedView Article
- Benner SA, Cohen MA, Gonnet GH: Amino acid substitution during functionally constrained divergent evolution of protein sequences. Protein Eng. 1994, 7 (11): 1323-1332. 10.1093/protein/7.11.1323.PubMedView Article
- Genedoc homepage. [http://www.nrbsc.org/gfx/genedoc/index.html]
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMed
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View Article
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12 (4): 357-358.PubMed
- Dyrlov Bendtsen J, Nielsen H, Von Heijne G, Brunak S: Improved Prediction of Signal Peptides: SignalP 3.0. J Mol Biol. 2004, 340 (4): 783-795. 10.1016/j.jmb.2004.05.028.View Article
- TMpred - prediction of transmembrane regions and orientation. [http://www.ch.embnet.org/software/TMPRED_form.html]
- Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, Langendijk-Genevaux PS, Pagni M, Sigrist CJ: The PROSITE database. Nucleic Acids Res. 2006, 34 (Database issue): D227-30. 10.1093/nar/gkj063.PubMed CentralPubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.