Physiol. Genomics Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Physiol. Genomics 27: 219-230, 2006. First published August 22, 2006; doi:10.1152/physiolgenomics.00322.2005
1094-8341/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Figures
Right arrow All Versions of this Article:
27/3/219    most recent
00322.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Saito, S.
Right arrow Articles by Shingai, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Saito, S.
Right arrow Articles by Shingai, R.
Received 28 December 2005; accepted in final form 17 August 2006.
Physiological Genomics 27:219-230 (2006)
1094-8341/06 $8.00 © 2006 American Physiological Society

Call For Papers: Comparative Genomics

Evolution of thermoTRP ion channel homologs in vertebrates

Shigeru Saito and Ryuzo Shingai

Laboratory of Bioscience, Faculty of Engineering, Iwate University, 4-3-5 Ueda, Morioka 020-8551, Japan


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
In mammalian thermosensation, nine temperature-sensitive ion channels that are activated by distinct temperature thresholds have been identified as thermosensors. These ion channels belong to the transient receptor potential (TRP) superfamily and are referred to as "thermoTRPs" (TRPV1, TRPV2, TRPV3, TRPV4, TRPM2, TRPM4, TRPM5, TRPM8, and TRPA1). To elucidate the evolutionary processes of thermoTRPs, we conducted comprehensive searches for mammalian thermoTRP gene homologs in the draft genome sequences of chicken (Gallus gallus), western clawed frog (Xenopus tropicalis), zebrafish (Danio rerio), and pufferfish (Fugu rubripes). Newly identified homologs were compared with known thermoTRPs, and phylogenetic analyses were conducted. Our comparative analyses revealed that most of the mammalian thermo-TRP members already existed in the common ancestor of fishes and tetrapods. Tetrapods shared almost the same repertoire, except that the western clawed frog expanded TRPV4s (six copies) and TRPM8s (two copies), which were diversified considerably. Comparisons of nonsynonymous and synonymous substitution rates among TRPV4s suggested that one copy of the TRPV4 channel in the western clawed frog retained its original function, while the other copies diversified and obtained slightly different properties. In fish lineages, several members of thermo-TRPs have duplicated in the whole genome duplication occurred in the ancestral ray-finned fish; however, some of the copies have subsequently been lost. Furthermore, fishes do not possess the three members of thermoTRPs existed in mammals, e.g., thermoTRPs activated by noxious heat, warm, and cool temperatures. Our results suggest that thermosensation mechanisms have changed through vertebrate evolution with respect to thermosensor repertoires.

temperature-sensitive ion channel; vertebrate; comparative genomics; gene duplication; phylogenetic analyses


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
VERTEBRATES SENSE AMBIENT temperature accurately and adapt to environmental temperature changes through behavioral and physiological responses. Temperature sensation is transmitted via the peripheral nerves that are located in the dorsal root ganglia and trigeminal ganglion (30). The axons of these neurons are extended to terminal organs such as skin and muscle where they are exposed to thermal and other sensory stimuli. Initiation of thermosensation signal transduction is thought to be mediated by temperature-sensitive ion channels expressed in sensory ganglion neurons (30). Currently, nine temperature-sensitive ion channels have been identified in mammals, mainly in humans, mice, and rats.

All nine known temperature-sensitive channels belong to the transient receptor potential (TRP) superfamily (5, 30, 31, 48, 49). Due to their ion channel properties, these nine channels are called "thermoTRPs." The TRP superfamily is further divided into subfamilies based on the similarity of their amino acid sequences. Among the nine thermoTRPs, four belong to the TRPV (V for vanilloid) subfamily (TRPV1, TRPV2, TRPV3, and TRPV4), four to the TRPM (M for melastatin) subfamily (TRPM2, TRPM4, TRPM5, and TRPM8), and one to the TRPA (A for ankyrin) subfamily (TRPA1). In mammals, the TRP superfamily currently consists of ~30 members that sense diverse physical and chemical stimuli (2, 5, 30, 31, 48, 49), suggesting that repeated gene duplications and subsequent functional diversifications have played fundamental roles in the evolution of TRP superfamily.

The nine thermoTRPs are activated by distinct temperature thresholds. TRPV1 is activated by temperatures >42°C, TRPV2 by >52°C, TRPV3 by >33°C, TRPV4 between 27 and 42°C, TRPM2 between 35 and 42°C, TRPM4 and TRPM5 between 15 and 35°C, TRPM8 by <25°C, and TRPA1 by <17°C, when overexpressed in cultured cells or Xenopus oocytes (7, 8, 22, 23, 3032, 39, 40, 44, 4851). Virtually the entire range of temperatures that mammals are exposed is covered by the temperature threshold of these nine thermo-TRPs. All thermoTRPs, except for the TRPM2, M4, and M5, are expressed in the peripheral nervous system. Therefore, six thermoTRPs (TRPV1–4, TRPM8, and TRPA1) are suggested to be involved in ambient and body temperature perception in mammals. For TRPV1, V3, and V4, knockout mice have been constructed, and behavioral responses have been examined (6, 21, 25). These mice showed deficits in responses to temperature consistent with temperature ranges of destructed thermo-TRPs. On the other hand, conflicting results have been reported for TRPA1 knockouts generated in two different laboratories (4, 20). One group claimed that TRPA1 knockout mice show deficit in response to noxious cold (20), while the other group did not observe such effects on cold temperature perception (4). In addition, such unequivocal observations have been reported when TRPA1 was overexpressed in cultured cells (17, 40). Therefore, the role for TRPA1 in cold temperature perception has been disputed.

Some of the thermoTRPs are expressed in tissues other than neuron such as skin in the case of TRPV1, V3, and V4 (10, 30, 33, 51). It has been suggested that keratinocytes may act as thermal receptors since synapses have not been found between keratinocytes and sensory termini. Putative signal transduction from keratinocytes has been proposed to be mediated by an ATP-gated channel (P2X3) that is present in sensory termini (30, 33). TRPM5 is expressed in the taste buds of the tongue and has been suggested to function downstream of the taste receptors and to be involved in the modulation of thermal sensitivity in sweet taste perception (41).

Interestingly, thermoTRPs has been reported to be activated by stimuli other than temperature. For example, TRPV1 is also activated by acidic pH and capsaicin, TRPV2 by growth factor molecules, TRPV4 by hypotonic stimulus, TRPM8 by menthol and icilin, and TRPA1 by icilin (2, 8, 22, 23, 3032, 40, 48). Therefore, it appears that thermoTRPs serve not only as thermosensors but can also transmit other kinds of sensations.

Mammalian thermoTRP homologs have recently been identified in invertebrates. For example, TRPA1 has been reported in sea squirt (Ciona intestinalis), fruit fly (Drosophila melanogaster), and nematode (Caenorhabditis elegans) (29, 49). In the case of D. melanogaster, TRPA1 has been shown to be activated by warmth (approximately >27°C, Ref. 45) and to be involved in thermotaxis (35). These findings indicate that TRPA1 arose and was involved in a thermosensation mechanism in an early stage of animal evolution. As mammalian TRPA1 is activated by cold (<17°C; Refs. 40, 45) rather than warmth, its activation temperature threshold must have changed in the course of animal evolution. On the other hand, mammalian orthologous genes for TRPV1–4 and TRPM2, M4, M5, M8 have not been found in C. intestinalis, D. melanogaster, and C. elegans genome sequences (29). These observations suggest that these thermoTRPs arose after the divergence of the urochordate and vertebrate lineages. However, the identification of very few thermoTRP homologs in nonmammalian species to date has impeded elucidation of the evolutionary processes of thermoTRPs in vertebrate lineages. Recently, however, the draft genome sequences have been released for several key model organisms from different vertebrate classes and have provided an opportunity for the genome-wide screening for thermoTRP homologs.

In the present study, we conducted comprehensive searches for the genes of mammalian thermoTRP homologs in the draft genome sequences of chicken (Gallus gallus), western clawed frog (Xenopus tropicalis, diploid species closely related to the tetraploid species Xenopus laevis), zebrafish (Danio rerio), and pufferfish (Fugu rubripes). Together with well-characterized thermoTRPs of humans and rodents, the thermoTRP homologs were compared among species belonging to the four major different vertebrate classes: mammals, birds, amphibians, and fishes. Phylogenetic analyses of thermoTRP homologs were then carried out with the aim of elucidating the evolutionary processes of thermoTRPs in vertebrates. Since the thermosensor repertoire directly affects thermosensation mechanism, comparative analyses of thermoTRP homologs supply basic information for understanding the evolutionary changes of the thermosensation mechanisms in vertebrate lineages.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Retrieving sequences of the thermoTRP homologs from the draft genome sequences.
The full-length cDNA sequences of nine known human thermoTRPs (TRPV1, TRPV2, TRPV3, TRPV4, TRPM2, TRPM4, TRPM5, TRPM8, and TRPA1) were retrieved from the GenBank (http://www.ncbi.nlm.nih.gov/Genbank/index.html; see Table 1 for their accession nos.). The amino acid sequences were obtained by translating cDNA sequences. Using these amino acid sequences as queries, we performed TBLASTN searches (3) with the E value set below 10–10 against the draft genome sequences of G. gallus (assembly built 1.1, released in July 2004) (13), X. tropicalis (assembly v. 4.1, released in August 2005), D. rerio (assembly Zv5, released in May 2005), and F. rubripes (assembly v. 4.0, released in October 2004) (14). The TBLASTN searches were performed based on the online genome sequence databases published by the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/mapview/) for G. gallus, by Ensembl (http://www.ensembl.org/) for D. rerio, and by the Department of Energy's Joint Genome Institute (JGI; http://genome.jgi-psf.org/) for X. tropicalis and F. rubripes. Sequences >200-amino-acid long were retained from the hits obtained by TBLASTN searches. Because all of the hits were located within the gene regions in the assembly annotations, we were able to retrieve the predicted mRNA sequences from the respective databases. In the case of TRPA1, hits that showed similarity only in the ankyrin (ANK) repeat domains residing in an ~600-amino-acid portion of the NH2-terminal were excluded since these genes do not possess the central transmembrane domains that are well conserved in the TRP ion channels (Supplementary Fig. 1). (The online version of this article contains supplemental data.) We thus obtained 36 predicted sequences from the assembly annotations of the respective draft genome sequences. These predicted mRNA sequences were aligned with previously reported thermoTRPs whose full-length cDNA sequences were available in the databases.


View this table:
[in this window]
[in a new window]

 
Table 1. Transient receptor potential genes used in the present study

 

Figure 1
View larger version (17K):
[in this window]
[in a new window]

 
Fig. 1. Phylogenetic tree of the vertebrate transient receptor potential-vanilloid (TRPV) subfamily. Amino acid sequences from ankyrin (ANK) repeat 1 to TRP domains (392 residues) were used for estimating the numbers of substitutions by applying the Jones-Taylor-Thornton (JTT) model. From the estimated numbers of amino acid substitutions among TRPVs, the phylogenetic tree was reconstructed by the minimum-evolution (ME) method. The statistical confidence (bootstrap value) is indicated next to each interior branch. Gene duplication event are designated by open diamonds. WC and AC frog denote western and African clawed frog, respectively.

 
To more accurately predict mRNA sequences, the following procedures were performed. When predicted sequences lacked the conserved portions of the TRP genes, we tried to find conserved sequences of the lacking portions in the respective draft genome sequences by TBLASTN search. When we obtained the hits that were located between the flanking exons in the genomic sequence, and exon/intron boundaries were demarcated by standard donor/acceptor splice sites without frame shift, the portions were included in the sequences to improve the annotations. Similarly, when predicted sequences possessed portions that did not exist in the counterparts of other species, and these portions were not supported by expressed sequence tags (ESTs) and could be omitted by standard donor/acceptor splice sites without resulting in a frame shift, we excluded the portions from the annotations. On the basis of these procedures, 18 out of 36 annotations were modified. The predicted mRNA sequences that were modified in the present study are shown in Supplementary Fig. 2. Using newly obtained and previously reported sequences as queries, we performed a second round of TBLASTN searches against the draft genome sequences of the respective species. There were no new hits in the second round of TBLASTN searches.


Figure 2
View larger version (7K):
[in this window]
[in a new window]

 
Fig. 2. Phylogenetic tree of vertebrate TRPV4s. Amino acid sequences from ANK repeat 1 to TRP domains (479 residues) were used for estimating the numbers of substitutions by applying the JTT model. From estimated numbers of amino acid substitutions, the phylogenetic tree was reconstructed by the ME method.

 
In the case of TRPA1, nonvertebrate species were used for phylogenetic analysis. The mRNA sequences of D. melanogaster, Anopheles gambiae, and C. elegans were retrieved from GenBank. The predicted mRNA sequence of C. intestinalis TRPA1 was retrieved from the genome sequence database at JGI.

In the present study, we used gene names according to the unified nomenclature for the TRP superfamily (24). The names of newly identified genes were assigned according to the names of the known genes contained in the same clusters in the phylogenetic tree. The thermoTRP homologs used in the present study are listed in Table 1.

Molecular phylogenetic and evolutionary analyses.
The amino acid sequences were deduced by translating mRNA sequences. Multiple sequence alignments were performed using the CLUSTAL W algorithm (12), with minor manual adjustments. Evolutionary distances between the amino acid or nucleotide sequences were calculated with the Jones-Taylor-Thornton (JTT) (16) or Tamura-Nei (42) models, respectively, after all alignment gap sites were eliminated. From the evolutionary distances thus determined, phylogenetic trees were reconstructed by the neighbor-joining (NJ) or minimum-evolution (ME) methods (36, 37). The statistical confidence of each branch in the phylogenetic tree was estimated by the bootstrap method (11) with 1,000 replications. The rates of nonsynonymous (dN) and synonymous (dS) substitutions, and their ratios (dN/dS) were estimated for all pairs of vertebrate TRPV4s by the modified Nei-Gojobori method (52) with the transition/transversion ratio equal to 1.5 and Jukes-Cantor correction (18). The standard error was estimated by the bootstrap method with 1,000 replications. For testing the purifying selection, Z-tests were performed under the null hypothesis of dS = dN. All of the data analyses were performed using MEGA 3 (19).

Synteny detection.
To identify the orthologous genes surrounding the thermoTRP genes, we performed best-hit gene analyses. For each comparison, the gene of one species (species 1) located in the vicinity of the thermoTRP gene was used as query, and a TBLASTN search was first performed against the genome sequence of another species (species 2). When the best-hit gene in the genome sequence of species 2 was located around the orthologous thermoTRP gene, we used the amino acid sequence of this gene as a query and performed a TBLASTN search against the genome sequence of species 1. If the best-hit gene in the second TBLASTN search was the same gene used as the query for the first TBLASTN search, we defined this gene pair as orthologous genes. To further examine whether the conserved synteny exists around the thermoTRP orthologs between different species, we subsequently performed the alignments of the genome sequences between the regions containing the thermoTRP orthologs using the PipMaker (38). When the results obtained by the PipMaker were consistent with those obtained by best hit gene analysis, we considered it a confirmation of our results.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
We first aligned the predicted amino acid sequences with those of known thermoTRPs. Some of the predicted sequences did not contain portions of exons that were apparently present in appropriate locations in the genome sequences. To obtain more accurate predictions, annotations were modified according to the procedures described in MATERIALS AND METHODS.

TRPV subfamily in vertebrates.
We identified four copies of mammalian thermoTRP homologs in chicken, nine in western clawed frog, two in zebrafish, and three in pufferfish (Table 1). The amino acid sequences of these TRPVs were aligned with those of TRPVs available in the databases (Supplementary Fig. 3), and the phylogenetic trees were reconstructed using conserved portion containing predicted ANK repeat, transmembrane (TM), pore loop (PL), and TRP domains (Supplementary Figs. 1 and 3). A phylogenetic tree reconstructed from amino acid sequences using the ME method is shown in Fig. 1. The vertebrate TRPVs were divided into four major clusters that were supported by 100% bootstrap values. The clusters are designated as TRPV1/2, TRPV3, TRPV4, and TRPV5/6 according to the known TRPs included in the respective clusters. The newly predicted genes identified in the present study were named according to the known genes contained in the respective clusters. Although TRPV5s and TRPV6s are not thermo-TRPs, we included these genes in the phylogenetic analysis; inclusion of TRPV5s and TRPV6s in the analysis enables them to be used as out-groups since they diverged first from the three remaining clusters.


Figure 3
View larger version (22K):
[in this window]
[in a new window]

 
Fig. 3. Gene orders surrounding TRPV1, TRPV2, and TRPV3 of 4 species belonging to different vertebrate classes. The lengths of the genomic regions are shown on the right. The arrow indicates the gene with the direction. The orthologous genes are connected by the lines. The open, filled, and striped arrows indicate TRPV1, TRPV2, and TRPV3 genes, respectively. Chr, chromosome.

 
In the TRPV1/2 cluster, tetrapod TRPV1s and TRPV2s were separated into different monophyletic clusters, and fish TRPV1/2s were clustered with the tetrapod TRPV1 cluster. However, when the phylogenetic tree was reconstructed from amino acid sequences using the NJ method, the fish TRPV1/2 cluster diverged first from the tetrapod TRPV1 and TRPV2 clusters in the TRPV1/2 cluster (Supplementary Fig. 4A). Consequently, the phylogenetic position of fish TRPV1/2s was not clearly resolved. Two copies of the pufferfish TRPV1/2s did not form a monophyletic cluster, and pufferfish TRPV1/2a clustered with zebrafish TRPV1/2 in both the NJ and ME trees. The TRPV3 cluster contained one of each copy of the tetrapod TRPV3s, and none of the fish TRPV homologs were contained in this cluster. The TRPV4 cluster diverged before the divergence between TRPV1/2 and TRPV3 clusters (the bootstrap value supporting this branch was 75% in the ME tree). The TRPV4 cluster contained one of each copy of the TRPV4s of mammals, chicken, and fishes, and six copies of TRPV4s of western clawed frog. Notably, one copy of TRPV4 of western clawed frog (designated as TRPV4a) was clustered with other vertebrate TRPV4s, and the other five copies of western clawed frog (TRPV4b–f) formed a monophyletic cluster that diverged first in the TRPV4 cluster in both the ME and NJ trees. However, when, a phylogenetic tree was reconstructed from the nucleotide sequences using the ME method, six copies of TRPV4s of the western clawed frog formed a monophyletic cluster with a 100% bootstrap value, suggesting that these six TRPV4s were acquired in the lineage leading to the western clawed frog (Supplementary Fig. 4B). Therefore, in the phylogenetic tree which was inferred from amino acid sequences (Fig. 1 and Supplementary Fig. 4A), the ancestral position of the monophyletic cluster containing TRPV4b–f of western clawed frog in the TRPV4 cluster were likely to be the result of considerable divergence between TRPV4b–f and the remaining vertebrate TRPV4s. To infer an accurate phylogenetic relationship of TRPV4 proteins, we reconstructed a phylogenetic tree from the amino acid sequences of vertebrate TRPV4s using fish TRPV4s as out-groups (Fig. 2). In this phylogenetic tree, the six copies of TRPV4s of the western clawed frog formed a monophyletic cluster, and TRPV4a diverged first in this cluster. The evolutionary rate accelerated considerably after the divergence between TRPV4a and the ancestral gene of TRPV4b–f, suggesting the relaxation of the functional constraint in the lineages leading to TRPV4b–f.


Figure 4
View larger version (22K):
[in this window]
[in a new window]

 
Fig. 4. Phylogenetic tree of vertebrate transient receptor potential-melastatin (TRPM) 2, TRPM4, TRPM5, and TRPM8 homologs. Amino acid sequences of the central conserved portion (from position 340 to 1283, 673 residues) were used for estimating the numbers of substitutions by the JTT model. The phylogenetic tree was reconstructed by the ME method based on the estimated numbers of amino acid substitutions.

 
To examine the extent of the functional constraint of TRPV4s of the western clawed frog, we estimated the dN and dS substitutions and their ratios (dN/dS) for each pair of TRPV4s (Table 2). The dN/dS ratios among human, mouse, rat, and chicken TRPV4s reported to be the functional channels (22), ranged from 0.023 to 0.141. The dN/dS ratios between TRPV4a of the western clawed frog and other vertebrate TRPV4s ranged from 0.063 to 0.081, suggesting that the functional constraint of TRPV4a is similar to those of other vertebrate TRPV4s. On the other hand, dN/dS ratios between TRPV4s of human, mouse, rat, and chicken and TRPV4b–f of the western clawed frog ranged from 0.118 to 0.167 (data not shown), and dN/dS ratios among six copies of TRPV4s of the western clawed frog ranged from 0.318 to 1.015 (Table 2), suggesting relaxation of the functional constraint among TRPV4b–f of the western clawed frog. In 11 of 15 pairs among TRPV4s of western clawed frog, dSs were significantly larger than dNs (Table 2). The six TRPV4 proteins of the western clawed frog were predicted to be at least 591 amino acids long (for TRPV4f) without intervening stop codons. Furthermore, core structural domains were predicted such as three ANK repeat, six TM, PL, and TRP domains (Supplementary Figs. 1 and 3), and thus these five TRPV4 proteins retained conserved structural features of functional TRPV channels. In addition, portions of the exons of TRPV4a, V4b, and V4d were supported by ESTs.


View this table:
[in this window]
[in a new window]

 
Table 2. Rates of nonsynonymous and synonymous substitutions among vertebrate TRPV4s

 
Most of the predicted TRPV proteins used in the present study possessed a conserved core structure such as predicted ANK repeat, TM, PL, and TRP domains (Supplementary Figs. 1 and 3). However, pufferfish TRPV1/2b protein lacked a portion between TM2 and TM3 domains, thus possessing only five TM domains, suggesting that its secondary structure may differ from that of other TRP channels.

Phylogenetic analysis of the vertebrates TRPV subfamily revealed that the chicken possesses one copy each of TRPV1, V2, V3, and V4, that the western clawed frog possesses one copy each of TRPV1, V2, and V3, and six copies of TRPV4, that the zebrafish possesses one copy each of TRPV1/2 and V4, and that the pufferfish possesses two copies of TRPV1/2 and one copy of TRPV4 (Fig. 1). The numbers of thermoTRPV homologs in each species are summarized in Table 3. Note that TRPV5 and TRPV6 are not thermoTRP and thus not included in the table.


View this table:
[in this window]
[in a new window]

 
Table 3. Number of thermoTRP homologs in vertebrates

 
To gain insight into the evolutionary processes that generated the TRPV homologs, the genomic locations and the gene orders surrounding the TRPVs were compared. In the genome sequences of the human, chicken, and western clawed frog, TRPV1, V2, and V3 were located within the same chromosome or scaffold. TRPV4, on the other hand, were located in a different chromosome or scaffold in these species. In the zebrafish genome sequence, TRPV1/2 and TRPV4s were located in the same chromosome, however, ~36.4 Mb away from each other. Therefore, TRPV4s and the other TRPVs were located at different genomic locations in all species examined. The gene orders around TRPV4 among human, chicken, and zebrafish were well conserved (data not shown). In the genome sequence of the western clawed frog, TRPV4a and TRPV4b–f were located in scaffold 17 and 330, respectively. The gene order around TRPV4a and other vertebrate TRPV4s were well conserved (data not shown). TRPV4b–f were located within an ~167-kb region and organized tandemly without intervening genes in the following order TRPV4b, V4c, V4d, V4e, and V4f. Several genes were located at each side of the TRPV4a and TRPV4b–f in the respective scaffolds, thus TRPV4a was not organized in tandem with TRPV4b–f in the genome sequence of western clawed frog.

TRPV1, V2, and V3 were located close to each other. Figure 3 shows the gene orders surrounding TRPV1, V2, and V3 among the human, chicken, western clawed frog, and zebrafish. In the genome sequences of the human, chicken, and western clawed frog, TRPV1 and TRPV3 were located adjacently. Furthermore, in the chicken and western clawed frog, TRPV2 was located near TRPV1 and TRPV3 (within 0.58 Mb). On the other hand, in human chromosome 17, TRPV2 was located ~12 Mb away from TRPV1 and TRPV3, further away from the latter two TRPVs. The gene orders around TRPVs were well conserved among the human, chicken, and western clawed frog. The gene order around TRPV1/2 of zebrafish was also conserved, although we could not find TRPV2 and TRPV3 orthologs in the fish genomes. The pufferfish TRPV1/2a and TRPV1/2b were located in scaffold 511 and 188, respectively, thus located in different genomic regions (Supplementary Fig. 5A). The gene orders around pufferfish TRPV1/2a and zebrafish TRPV1/2 were conserved (data not shown), suggesting these two genes are orthologs as suggested by phylogenetic analysis (Fig. 1).


Figure 5
View larger version (11K):
[in this window]
[in a new window]

 
Fig. 5. Phylogenetic tree of animal transient receptor potential-ankyrin (TRPA) 1 homologs. Amino acid sequences from ANK repeat 1 to transmembrane (TM) 6 domains (826 residues) were used. The numbers of amino acid substitutions were estimated using the JTT model, and the phylogenetic tree was reconstructed by the ME method.

 
TRPM2, TRPM4, TRPM5, and TRPM8 homologs of vertebrates.
The TBLASTN searches identified three copies of thermoTRPMs in chicken and pufferfish, four in western clawed frog, and five in zebrafish (Table 1). The amino acid sequences of these genes were aligned with those of previously reported thermoTRPMs (Supplementary Fig. 6). Then, a phylogenetic tree was reconstructed from the amino acid sequences using the ME method (Fig. 4). The vertebrate TRPMs were divided into four major clusters that were supported by ≥97% bootstrap values. The clusters were designated as TRPM2, TRPM4, TRPM5, and TRPM8 according to the names of the known genes contained in the respective clusters. TRPM2 and TRPM8 clusters diverged first from remaining two clusters. All four clusters contained tetrapod and fish TRPMs, except that the TRPM8 cluster did not contain fish TRPM8s. The TRPM8 cluster contained two genes of western clawed frog that formed a monophyletic cluster. The extent of the divergence between these two copies corresponded to that of the divergence between TRPM8s of chicken and mammals. TRPM8a and TRPM8b of the western clawed frog were arranged tandemly in scaffold 15, suggesting that these genes were produced by tandem gene duplication event.


Figure 6
View larger version (29K):
[in this window]
[in a new window]

 
Fig. 6. Evolutionary processes of the vertebrate thermoTRP homologs. The inferred evolutionary events (gene duplication and gene loss) are indicated on the respective branches. The gene duplication event of TRPV1/2 can be assigned either before or after the divergence between the fish and tetrapod lineages, thus these alternative hypotheses are shown with dashed arrows. Note that if duplication of TRPV1/2 occurred in the common ancestors of fishes and tetrapods, TRPV2 must have been lost in the common ancestor of pufferfish and zebrafish.

 
TRPM4 and TRPM5 clusters did not contain TRPMs of chicken and western clawed frog, respectively (Fig. 4). On the other hand, the TRPM4 cluster contained three copies of zebrafish TRPM4s, and one copy (TRPM4a) formed monophyletic cluster with pufferfish TRPM4. When a phylogenetic tree was reconstructed using the NJ method, zebrafish TRPM4b and TRPM4c formed a monophyletic cluster with pufferfish TRPM4 (bootstrap value supporting this branch was 59%). Therefore, we examined the gene order around TRPM4s and found conserved synteny between zebrafish TRPM4a and pufferfish TRPM4 (data not shown), suggesting that the pufferfish TRPM4 and zebrafish TRPM4a are the orthologous genes as inferred by the ME tree in Fig. 4. Zebrafish TRPM4a was located in chromosome 3, whereas zebrafish TRPM4b and TRPM4c were in chromosome 12 (Supplementary Fig. 5B). Zebrafish TRPM4b and TRPM4c, which were the most closely related among the three copies of zebrafish TRPM4 sequences, were located relatively close to one another with just four intervening genes between them.

The numbers of vertebrate thermoTRPM homologs are summarized in Table 3. All of the vertebrate TRPM homologs examined in the present study possessed six predicted TM domains and one TRP domain (Supplementary Figs. 1 and 6).

Several TRP channels (TRPM8, TRPV1, TRPV3) have recently been reported to exhibit a shift in the voltage dependence of their channel activations (shift their current-voltage relationships) in response to temperature stimuli; and this trait has been proposed to be coupled with the temperature sensitivity of thermoTRPs (27, 47). Although the voltage sensor has not yet been identified in thermoTRPs, the TM4 domain may be a good candidate since, in the case of the Shaker K+ channel, positively charged amino acids in the TM4 domain were shown to act as a voltage sensor (1). In this respect, it is worth noting that the positively charged arginine residue located in the predicted TM4 domain of the TRPMs and the TRPVs was entirely conserved (Supplementary Figs. 3 and 6).

TRPA1 homologs of vertebrates.
The searches for TRPA1 homologs identified one copy each in the chicken, western clawed frog, and pufferfish, and two in zebrafish (Table 1). The amino acid sequences of these TRPA1s were aligned with those of TRPA1 of mammals, sea squirt (C. intestinalis), insects (D. melanogaster and A. gambiae), and nematode (C. elegans) (Supplementary Fig. 7). Phylogenetic trees were then reconstructed from the amino acid sequences by the NJ and ME methods (Fig. 5). The topologies of the phylogenetic trees inferred by both methods were consistent. The phylogenetic relationship of vertebrate TRPA1s was consistent with the established speciation process of vertebrates. Two copies of zebrafish TRPA1s did not form a monophyletic cluster, and zebrafish TRPA1b clustered with pufferfish TRPA1. Zebrafish TRPA1a and TRPA1b were located in chromosome 2 and 24, respectively. It is worth noting that sea squirt TRPA1 clustered with insect TRPA1s with high bootstrap values (93% and 86% in the NJ and ME trees, respectively) consistent with Okamura et al.'s findings (29), although the sea squirt is considered to be more closely related to vertebrates than to insects. All of the vertebrate TRPA1s possessed the predicted 15 ANK repeat domains and six TM domains (Supplementary Figs. 1 and 7). The numbers of vertebrate TRPA 1 homologs are summarized in Table 3.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Evolution of thermoTRP homologs in vertebrates.
Our comprehensive searches for mammalian thermoTRP homologs in the draft genome sequences of the chicken, western clawed frog, zebrafish, and pufferfish revealed that the repertories of the homologs differed among species belonging to different classes (Table 3). We can be certain that our inability to find fish TRPV2, TRPV3, and TRPM8 was not due to inefficiency of the detection method used or incompleteness of the draft genome sequences for the following reasons. First, the criteria used for the TBLASTN searches were sensitive enough to detect other members belonging to the same subfamily. For example, when TRPV1 was used as a query, our search detected TRPV1, V2, V3, and V4 if they existed in the genome sequence. Second, we searched against the draft genome sequences of two fish species; thus it is unlikely that TRPV2, TRPV3, and TRPM8 are located in the genomic regions which coincidentally have not been sequenced yet in either of the species. Finally, the gene orders around TRPV1, V2, and V3 in the genome sequences of tetrapods and that of TRPV1/2 of zebrafish were well conserved (Fig. 3). To be missed under these circumstances, both TRPV2 and TRPV3 would have to have independently translocated to the regions that have not yet been sequenced in both fish genomes, a coincidence that is highly improbable.

Using the thermoTRP homologs detected in the present study, we performed phylogenetic analyses to elucidate evolutionary processes in vertebrate lineages. In the phylogenetic tree of the TRPV subfamily, the TRPV4 cluster was located at a position ancestral to the TRPV3 cluster (Fig. 1). Comparison of the genomic locations of the TRPVs revealed that TRPV4 was located at a different genomic location in tetrapods and fish, although TRPV1, V2, and V3 were closely located (Fig. 3). This arrangement is consistent with the hypothesis that TRPV4s diverged before the divergence of TRPV1, V2, and V3.

Although the TRPV4 cluster contained tetrapod and fish genes, the TRPV3 cluster did not contain fish genes. In addition, both the ME and NJ trees strongly indicate that fish TRPV1/2s were produced after the divergence between TRPV3 and TRPV1/2 clusters (Fig. 1 and Supplementary Fig. 4A). Therefore, we conclude that TRPV3 emerged in a common ancestor of fishes and tetrapods and has been lost in fish lineages. The inferred evolutionary processes of thermoTRP homologs of vertebrates are shown in Fig. 6.

The time of the gene duplication event producing TRPV1 and TRPV2 could not be unambiguously determined. In the NJ tree inferred from amino acid sequences, the gene duplication event appeared to occur after the divergence between fish and tetrapod lineages (Supplementary Fig. 4A). In the ME tree, however, the fish TRPV1/2s clustered with tetrapod TRPV1s (Fig. 1), suggesting that TRPV1 and TRPV2 already existed in the common ancestor of fishes and tetrapods and that TRPV2 was subsequently lost in the fish lineage. When the genomic location of zebrafish TRPV1/2 was examined, its genomic position corresponded to that of tetrapod TRPV1 (Fig. 3). In the genome sequences of the chicken and western clawed frog, TRPV1 and TRPV3 were located adjacently, whereas TRPV2 was located several genes away from TRPV1 and TRPV3. Taking into consideration the findings that the gene orders around the tetrapod TRPV2 were also conserved near the zebrafish TRPV1/2, we find it likely that three TRPVs may have already existed in the common ancestor of fishes and tetrapods and that they were organized in tandem, TRPV2-TRPV1-TRPV3, in the genome sequence. If so, then several genes were later inserted between TRPV1 and TRPV2 before the divergence of fish and tetrapod lineages, and TRPV2 and TRPV3 were subsequently lost in the common ancestor of zebrafish and pufferfish. Accordingly, the genomic organization of TRPVs is consistent with the phylogenetic relationship reconstructed by the ME method.

In the phylogenetic tree of the TRPM subfamily, all clusters except for the TRPM8 cluster contained both tetrapod and fish genes (Fig. 4). This situation suggests that the gene duplication events producing TRPM2, M4, M5, and M8 occurred before the divergence of fish and tetrapod lineages and that TRPM8s were lost in the common ancestor of zebrafish and pufferfish (Fig. 6). Meanwhile our findings indicate that TRPM4 and TRPM5 have been lost in the chicken and western clawed frog, respectively. Since TRPM4 and TRPM5 are reported to exhibit similar ion channel properties (41), they might be able to compensate for each other's functions if either were lost. The zebrafish possesses three copies of TRPM4s (termed a, b, and c). Amino acid sequence similarity and close proximity of the genomic location of TRPM4b and TRPM4c suggest they were produced by a tandem gene duplication event. On the other hands, TRPM4a did not cluster together with TRPM4b and TRPM4c but, rather, clustered with pufferfish TRPM4 (Fig. 4), suggesting that the gene duplication event occurred in the common ancestor of zebrafish and pufferfish and that one copy was lost in the pufferfish lineage (Fig. 6).

With regard to fish lineages, the timing of gene duplication events appeared to be similar for TRPV1/2, TRPM4, and TRPA1. Specifically, these gene duplication events occurred just before the divergence of zebrafish and pufferfish lineages (Figs. 1, 4, and 5). In this respect, it has been proposed that the whole genome duplication event occurred in the common ancestor of zebrafish and pufferfish and greatly affected the fish genome structures (9, 43). To examine whether three gene duplications of thermoTRPs were associated with genome duplication event, we compared the gene orders around each pair of the duplicated copies of TRPV1/2, TRPM4, and TRPA1. We found that several genes located adjacent to the duplicated copies of TRPs showed high amino acid sequence similarity to one another (from 47 to 88%, Supplementary Fig. 5) in the case of TRPV1/2s and TRPM4s, although a similar relationship was not found in the case of TRPA1a and TRPA1b. These findings suggest that duplicated copies of TRPV1/2s and TRPM4s are located in paralogous genomic regions. Therefore, it is likely that TRPV1/2s and TRPM4s were produced by the whole genome duplication event occurring in the common ancestor of zebrafish and pufferfish.

Although gene duplication events have occurred in the fish lineage, fishes have also subsequently lost TRPV1/2, TRPM4, and TRPA1. In the phylogenetic trees of the TRPV, TRPM, and TRPA1 subfamilies, the two duplicated copies from the same species did not form monophyletic clusters (Figs. 1, 4, and 5), suggesting that gene loss occurred after the gene duplication events (Fig. 6). Furthermore, pufferfish TRPV1/2b has deletion in the predicted TM domains; thus this TRPV may be a pseudogene, or at least it may possess the different ion channel properties (Supplementary Fig. 3). In addition, TRPV2, TRPV3, and TRPM8 were lost in the common ancestor of zebrafish and pufferfish (Fig. 6). Therefore, in the fish lineages, the thermoTRP homologs seem to have a tendency to be lost. If the activation temperature thresholds of fish thermoTRPs are similar to those of mammalian thermoTRPs, the activation temperature ranges of three thermoTRPs that fishes possess (TRPV1/2, TRPV4, and TRPA1) could roughly cover the entire range of ambient water temperature (30, 31, 48). In this respect, fishes may possess reduced sets of thermoTRP channels compared with those of tetrapods. Alternatively, other as yet uncharacterized channels may serve as the thermosensors in fishes; if so fishes would possess a different thermoTRP repertoire from mammals. Hence, examination of the activation temperature thresholds of fish thermoTRP homologs is an interesting subject for future investigation.

In the present study, we found that most of the thermoTRPs that mammals possess already existed in the common ancestor of tetrapods and ray-finned fishes and that thermoTRP homolog repertoires have changed in the lineages leading to different vertebrate classes (Table 3 and Fig. 6). The absence of orthologous genes of mammalian thermoTRPs beyond TRPA1 in urochordates (sea squirts) (29) suggests that most of the vertebrate thermoTRPs emerged after the divergence of vertebrate and urochordate lineages but before that of ray- and lobe-finned fish lineages. In this respect, identification and comparison of the thermoTRP homologs in organisms that diverged in the early stage of vertebrate evolution, such as cartilaginous fishes (sharks), and jawless fishes (lampreys), and Cephalochorda (amphioxuses), will help us to better understand of the origin of vertebrate thermoTRPs.

Implications for functional evolution of vertebrate thermo-TRPs.
The only mammalian thermoTRP ortholog known to be present in sea squirt (C. intestinalis) is TRPA1 (29), and TRPA1 orthologs exist in fruit fly (D. melanogaster) and nematode (C. elegans). Thus it can be concluded that TRPA1 emerged quite early in animal evolution. The finding that C. intestinalis TRPA1 is more similar to insect TRPA1 than to vertebrate TRPA1s is inconsistent with the speciation process as reported by Okamura et al. (29; see also Fig. 5). The activation temperature thresholds of TRPA1s are different between D. melanogaster and human when they were overexpressed in Xenopus oocytes (45). The functional portions responsible for this interspecies difference are not yet known. Examination of the temperature threshold of C. intestinalis TRPA1 may supply insight into change in activation temperature threshold of TRPA1 in the course of animal evolution.

Our phylogenetic analysis of vertebrate thermoTRP homologs revealed that all of the thermoTRPs that mammals possess already existed in the common ancestor of tetrapods (Fig. 6), indicating that the genetic basis for thermosensation of warm-blooded animals, in terms of their thermosensor repertoires, was acquired in the common ancestor of tetrapods. Acquisition of homeothermism occurred after the divergence of amphibians and bird/mammal lineages and may have significantly affected the thermosensation mechanism. Since amphibians diverged just before the acquisition of homeothermism characteristic of mammals and birds, comparisons of thermoTRPs between amphibians and warm-blooded animal (birds and mammals) are valuable for better understanding the mechanism of thermosensation. In this respect, it is worth noting that the NH2- and COOH-terminal portions of the predicted amino acid sequences of TRPV3 of the western clawed frog were truncated compared with those of mammalian TRPV3s (Supplementary Fig. 3). Mammalian TRPV3 channels are activated by temperatures >33°C and are involved in sensing temperatures around the body temperature of warm-blooded animals (25, 33, 39, 51).

Since the predicted amino acid sequence of TRPV3 of the western clawed frog was not supported by the ESTs, prediction of the NH2 and COOH-terminal portions may not be correct. Therefore, we performed TBLASTN searches with the E value set <10–5 against the draft genome sequences of the western clawed frog using the amino acid sequences of the NH2 and COOH-terminal portions of the mammalian and chicken TRPV3s as queries; however, no hits were obtained. These results suggest that the NH2 and COOH-terminal portions of TRPV3 of the western clawed frog differ from those of TRPV3s of mammals and birds. Truncations of the COOH-terminal portion of TRPV1 reduced its temperature threshold when expressed in cultured cells (46), thus the COOH-terminal portion is suggested to be responsible for determining the temperature thresholds of TRPV channels. Determination of the cDNA sequences of amphibian TRPV3s and examination of their temperature thresholds may provide intriguing insights into the evolution of thermosensory function relevant to the evolutionary acquisition of homeothermism of mammals and birds.

Interestingly, western clawed frog not only possesses homologs for most of the mammalian thermoTRPs but has also gained and diversified TRPM8 and TRPV4 homologs (Figs. 2, 4 and 6). TRPM8a and TRPM8b, as well as TRPV4b–f, are located in tandem in the western clawed frog genome sequence, respectively; thus tandem duplication is the main mechanism that produced diversified copies of TRPV4 and TRPM8 homologs.

Phylogenetic analysis and comparisons of the dN/dS ratios among vertebrate TRPV4s revealed that the functional constraint of TRPV4a is similar to that of other vertebrate TRPV4s, which have been reported to be activated by temperature and hypotonic stimuli in in vitro expression experiments (22, 50) and thus are considered to be functional channels (Table 2). The other five copies (TRPV4b–f) have diverged considerably, and their dN/dS ratios suggest that substantial relaxation of the functional constraint has occurred in these lineages. However, in most of the comparisons, dSs are significantly larger than dNs, suggesting that most, if not all, of the TRPV4s have been subjected to purifying selection. In addition, all of the predicted amino acid sequences of TRPV4s in the western clawed frog possessed conserved structural domains, and EST data suggest that the mRNAs of TRPV4a, TRPV4b, and TRPV4d genes are transcribed. These observations suggest that diversified copies of TRPV4s of the western clawed frog maintain their function. After the divergence of TRPV4a and the ancestral gene of TRPV4b–f, considerable amino acid substitutions accumulated in the ancestral gene of TRPV4b–f (Fig. 2), which may contribute to the acquisition of different properties, while TRPV4a retained its original function. Subsequently, the ancestral gene of TRPV4b–f gained copies through repeated tandem gene duplications, and further accumulations of the amino acid substitutions in each of the copies led to the diversified copies of TRPV4s in the western clawed frog.

Taking into consideration the fact that other vertebrate species possess only one copy of TRPV4, TRPV4s of the western clawed frog may have conferred some advantages that might be related to the traits specific to amphibians. TRPV4s of mammals and chickens are activated not only by temperature but also by hypotonic stimuli and are thought to serve as an osmosensor (22, 28). As amphibians have a unique subterrestrial life cycle that require them to be well adapted to both aquatic and terrestrial environments, the western clawed frog may utilize diversified copies of TRPV4s as osmosensors to accurately sense the osmotic pressure of its skin. Further examination of the functional properties of TRPV4s in amphibians is essential for understanding the functional evolution of TRPV4s.

In summary, in the present study we found that the repertoires of thermoTRP homologs have changed through vertebrate evolution. Since thermoTRPs are involved in thermosensation as well as other kinds of sensory detection (2, 30, 31, 48), variability of repertoires may be associated with adaptation of the organisms to their respective habitat environments. Further examination of the functions of thermoTRP homologs in organisms other than mammals will enhance our understanding of how the functional changes of thermoTRPs have contributed to the vertebrate adaptation mechanisms.


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This work was supported by a Grant-in-Aid for the 21st Century Center of Excellence Program from the Ministry of Education, Culture, Sports, Science and Technology of Japan.


    ACKNOWLEDGMENTS
 
We thank Masafumi Nozawa, Claire T. Saito, and two anonymous reviewers for fruitful discussions and critical reading of the manuscript. The data for the G. gallus genome were provided by the Genome Sequencing Center, Washington University School of Medicine (St. Louis). The genome sequences data of X. tropicalis were provided by the Department of Energy's JGI. The sequence data for the Danio rerio genome were produced by the Zebrafish Sequencing Group at the Sanger Institute and can be obtained from http://www.ensembl.org/Danio_rerio. The data of the F. rubripes genome have been provided freely by the Fugu Genome Consortium for use in this publication/correspondence only.


    FOOTNOTES
 
Address for reprint requests and other correspondence: S. Saito, Laboratory of Bioscience, Faculty of Engineering, Iwate Univ., 4-3-5 Ueda, Morioka 020-8551, (e-mail: shigeru{at}iwate-u.ac.jp)

Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Aggarwal SK, MacKinnon R. Contribution of the S4 segment to gating charge in the Shaker K+ channel. Neuron 16: 1169–1177, 1996.[CrossRef][ISI][Medline]
  2. Alexander SP, Mathie A, Peters JA. Guide to receptors and channels, 1st edition (2005 revision). Br J Pharmacol 144: S1–S128, 2005.[CrossRef][ISI][Medline]
  3. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402, 1997.[Abstract/Free Full Text]
  4. Bautista DM, Jordt SE, Nikai T, Tsuruda PR, Read AJ, Poblete J, Yamoah EN, Basbaum AI, Julius D. TRPA1 mediates the inflammatory actions of environmental irritants and proalgesic agents. Cell 124: 1269–1282, 2006.[CrossRef][ISI][Medline]
  5. Birnbaumer L, Yidirim E, Abramowitz J. A comparison of the genes coding for canonical TRP channels and their M, V and P relatives. Cell Calcium 33: 419–432, 2003.[CrossRef][ISI][Medline]
  6. Caterina MJ, Leffler A, Malmberg AB, Martin WJ, Trafton J, Petersen-Zeitz KR, Koltzenburg M, Basbaum AI, Julius D. Impaired nociception and pain sensation in mice lacking the capsaicin receptor. Science 288: 306–313, 2000.[Abstract/Free Full Text]
  7. Caterina MJ, Rosen TA, Tominaga M, Brake AJ, Julius D. A capsaicin-receptor homologue with a high threshold for noxious heat. Nature 398: 436–441, 1999.[CrossRef][Medline]
  8. Caterina MJ, Schumacher MA, Tominaga M, Rosen TA, Levine JD, Julius D. The capsaicin receptor: a heat-activated ion channel in the pain pathway. Nature 389: 816–824, 1997.[CrossRef][Medline]
  9. Christoffels A, Koh EG, Chia JM, Brenner S, Aparicio S, Venkatesh B. Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol 21: 1146–1151, 2004.[Abstract/Free Full Text]
  10. Denda M, Fuziwara S, Inoue K, Denda S, Akamatsu H, Tomitaka A, Matsunaga K. Immunoreactivity of VR1 on epidermal keratinocyte of human skin. Biochem Biophys Res Commun 285: 1250–1252, 2001.[CrossRef][ISI][Medline]
  11. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783–791, 1985.[CrossRef][ISI]
  12. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680, 1994.[Abstract/Free Full Text]
  13. International Chicken Genome Sequencing Consortium. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature 432: 659–777, 2004.[CrossRef][Medline]
  14. International Fugu Genome Consortium. Whole-Genome Shotgun assembly and analysis of the genome of Fugu rubripes. Science 297: 1301–1310, 2002.[Abstract/Free Full Text]
  15. Jaquemar D, Schenker T, Trueb B. An ankyrin-like protein with transmembrane domains is specifically lost after oncogenic transformation of human fibroblasts. J Biol Chem 274: 7325–7333, 1999.[Abstract/Free Full Text]
  16. Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 8: 275–282, 1992.[Abstract/Free Full Text]
  17. Jordt SE, Bautista DM, Chuang H, McKemy DD, Zygmunt PM, Högestätt ED, Meng ID, Julius D. Mustard oils and cannabinoids excite sensory nerve fibres through the TRP channel ANKTM1. Nature 427: 260–265, 2004.[CrossRef][Medline]
  18. Jukes TH, Cantor CR. Evolution of protein molecules. In: Mammalian Protein Metabolism III. New York, Academic, 1969, p. 21–132.
  19. Kumar S, Tamura K, Nei M. MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 5: 150–163, 2004.[Abstract/Free Full Text]
  20. Kwan KY, Allchorne AJ, Vollrath MA, Christensen AP, Zhang DS, Woolf CJ, Corey DP. TRPA1 contributes to cold, mechanical, and chemical nociception but is not essential for hair-cell transduction. Neuron 50: 277–289, 2006.[CrossRef][ISI][Medline]
  21. Lee H, Iida T, Mizuno A, Suzuki M, Caterina MJ. Altered thermal selection behavior in mice lacking transient receptor potential vanilloid 4. J Neurosci 25: 1304–1310, 2005.[Abstract/Free Full Text]
  22. Liedtke W, Choe Y, Martí-Renom MA, Bell AM, Denis CS, Sali A, Hudspeth AJ, Friedman JM. Vanilloid receptor-related osmotically activated channel (VR-OAC), a candidate vertebrate osmoreceptor. Cell 103: 525–523, 2000.[CrossRef][ISI][Medline]
  23. McKemy DD, Neuhausser WM, Julius D. Identification of a cold receptor reveals a general role for TRP channels in thermosensation. Nature 416: 52–58, 2002.[CrossRef][Medline]
  24. Montell C, Birnbaumer L, Flockerzi V, Bindels RJ, Bruford EA, Caterina MJ, Clapham DE, Harteneck C, Heller S, Julius D, Kojima I, Mori Y, Penner R, Prawitt D, Scharenberg AM, Schultz G, Shimizu N, Zhu MX. A unified nomenclature for the superfamily of TRP cation channels. Mol Cell 9: 229–231, 2002.[CrossRef][ISI][Medline]
  25. Moqrich A, Hwang SW, Earley TJ, Petrus MJ, Murray AN, Spencer KSR, Andahazy M, Story GM, Patapoutian A. Impaired thermosensation in mice lacking TRPV3, a heat and camphor sensor in the skin. Science 307: 1468–1472, 2005.[Abstract/Free Full Text]
  26. Nilius B, Mahieu F, Prenen J, Janssens A, Owsianik G, Vennekens R, Voets T. The Ca2+-activated cation channel TRPM4 is regulated by phosphatidylinositol 4,5-biphosphate. EMBO J 25: 467–478, 2006.[CrossRef][ISI][Medline]
  27. Nilius B, Talavera K, Owsianik G, Prenen J, Droogmans G, Voets T. Gating of TRP channels: a voltage connection? J Physiol 567: 35–44, 2005.[Abstract/Free Full Text]
  28. Nilius B, Vriens J, Prenen J, Droogmans G, Voets T. TRPV4 calcium entry channel: a paradigm for gating diversity. Am J Physiol Cell Physiol 286: C195–C205, 2004.[Abstract/Free Full Text]
  29. Okamura Y, Nishino A, Murata Y, Nakajo K, Iwasaki H, Ohtsuka Y, Tanaka-Kunishima M, Takahashi N, Hara Y, Yoshida T, Nishida M, Okado H, Watari H, Meinertzhagen IA, Satoh N, Takahashi K, Satou Y, Okada Y, Mori Y. Comprehensive analysis of the ascidian genome reveals novel insights into the molecular evolution of ion channel genes. Physiol Genomics 22: 269–282, 2005.[Abstract/Free Full Text]
  30. Patapoutian A, Peier M, Story GM, Viswanath V. ThermoTRP channels and beyond: mechanisms of temperature sensation. Nat Rev Neurosci 4: 529–539, 2003.[CrossRef][ISI][Medline]
  31. Pedersen SF, Owsianik G, Nilius B. TRP channels: an overview. Cell Calcium 38: 233–252, 2005.[CrossRef][ISI][Medline]
  32. Peier AM, Moqrich A, Hergarden AC, Reeve AJ, Andersson DA, Story GM, Earley TJ, Dragoni I, McIntyre P, Bevan S, Patapoutian A. A TRP channel that senses cold stimuli and menthol. Cell 108: 705–715, 2002.[CrossRef][ISI][Medline]
  33. Peier AM, Reeve AJ, Andersson DA, Moqrich A, Earley TJ, Hergarden AC, Story GM, Colley S, Hogenesch JB, McIntyre P, Bevan S, Patapoutian A. A heat-sensitive TRP channel expressed in keratinocytes. Science 296: 2046–2049, 2002.[Abstract/Free Full Text]
  34. Qiu A, Hogstrand C. Functional characterisation and genomic analysis of an epithelial calcium channel (ECaC) from pufferfish, Fugu rubripes. Gene 342: 113–123, 2004.[CrossRef][ISI][Medline]
  35. Rosenzweig M, Brennan KM, Tayler TD, Phelps PO, Patapoutian A, Garrity PA. The Drosophila ortholog of vertebrate TRPA1 regulates thermotaxis. Genes Dev 19: 419–424, 2005.[Abstract/Free Full Text]
  36. Rzhetsky A, Nei M. A simple method for estimating and testing minimum evolution trees. Mol Biol Evol 9: 945–967, 1992.[ISI]
  37. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425, 1987.[Abstract]
  38. Schwartz S, Zhang Z, Frazer KA, Smit A, Riemer C, Bouck J, Gibbs R, Hardison R, Miller W. PipMaker-A web server for aligning two genomic DNA sequences. Genome Res 10: 577–586, 2000.[Abstract/Free Full Text]
  39. Smith GD, Gunthorpe MJ, Kelsell RE, Hayes PD, Reilly P, Facer P, Wright JE, Jermank JC, Walhin JP, Ooi L, Egerton J, Charles KJ, Smart D, Randall AD, Anand P, Davis JB. TRPV3 is a temperature-sensitive vanilloid receptor-like protein. Nature 418: 186–190, 2002.[CrossRef][Medline]
  40. Story GM, Peier AM, Reeve AJ, Eid SR, Mosbacher J, Hricik TR, Earley TJ, Hergarden AC, Andersson DA, Hwang SW, McIntyre P, Jegla T, Bevan S, Patapoutian A. ANKTM1, a TRP-like channel expressed in nociceptive neurons, is activated by cold temperatures. Cell 112: 819–829, 2003.[CrossRef][ISI][Medline]
  41. Talavera K, Yasumatsu K, Voets T, Droogmans G, Shigemura N, Ninomiya Y, Margolskee RF, Nilius B. Heat activation of TRPM5 underlies thermal sensitivity of sweet taste. Nature 483: 1022–1025, 2005.
  42. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: 512–526, 1993.[Abstract]
  43. Taylor JS, Braasch I, Frickey T, Meyer A, Peer YV. Genome duplication, a trait shared by 22,000 species of ray-finned fish. Genome Res 13: 382–390, 2003.[Abstract/Free Full Text]
  44. Togashi K, Hara Y, Tominaga T, Higashi T, Konishi Y, Mori Y, Tominaga M. TRPM2 activation by cyclic ADP-ribose at body temperature is involved in insulin secretion. EMBO J 25: 1804–1815, 2006.[CrossRef][ISI][Medline]
  45. Viswanath V, Story GM, Peier AM, Petrus MJ, Lee VM, Hwang SW, Patapoutian A, Jegla T. Opposite thermosensor in fruitfly and mouse. Nature 423: 822–823, 2003.[CrossRef][Medline]
  46. Vlachová V, Teisinger J, Susánková K, Lyfenko A, Ettrich R, Vyklicky L. Functional role of C-terminal cytoplasmic tail of rat vanilloid receptor 1. J Neurosci 23: 1340–1350, 2003.[Abstract/Free Full Text]
  47. Voets T, Droogmans G, Wissenbach U, Janssens A, Flockerzi V, Nilius B. The principle of temperature-dependent gating in cold- and heat-sensitive TRP channels. Nature 430: 748–754, 2004.[CrossRef][Medline]
  48. Voets T, Talavera K, Owsianik G, Nilius B. Sensing with TRP channels. Nature Chem Biol 1: 85–92, 2005.
  49. Vriens J, Owsianik G, Voets T, Droogmans G, Nilius B. Invertebrate TRP proteins as functional models for mammalian channels. Pflügers Arch 449: 213–226, 2004.[ISI][Medline]
  50. Watanabe H, Vriens J, Suh SH, Benham CD, Droogmans G, Nilius B. Heat-evoked activation of TRPV4 channels in a HEK293 cell expression system and in native mouse aorta endothelial cells. J Biol Chem 277: 47044–47051, 2002.[Abstract/Free Full Text]
  51. Xu H, Ramsey IS, Kotecha SA, Moran MM, Chong JA, Lawson D, Ge P, Lilly J, Silos-Santiago I, Xie Y, DiStefano PS, Curtis R, Clapham DE. TRPV3 is a calcium-permeable temperature-sensitive cation channel. Nature 418: 181–186, 2002.[CrossRef][Medline]
  52. Zhang J, Rosenberg HF, Nei M. Positive Darwinian selection after gene duplication in primate ribonuclease genes. Proc Natl Acad Sci USA 95: 3708–3713, 1998.[Abstract/Free Full Text]
  53. Zhang Z, Okawa H, Wang Y, Liman ER. Phosphatidylinositol 4,5-bisphosphate rescues TRPM4 channels from desensitization. J Biol Chem 280: 39185–39192, 2005.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
J. Neurosci.Home page
D. A. Prober, S. Zimmerman, B. R. Myers, B. M. McDermott Jr, S.-H. Kim, S. Caron, J. Rihel, L. Solnica-Krezel, D. Julius, A. J. Hudspeth, et al.
Zebrafish TRPA1 Channels Are Required for Chemosensation But Not for Thermosensation or Mechanosensory Hair Cell Function
J. Neurosci., October 1, 2008; 28(40): 10102 - 10110.
[Abstract] [Full Text] [PDF]


Home page
FASEB J.