MBE Advance Access originally published online on February 16, 2006
Molecular Biology and Evolution 2006 23(5):964-972; doi:10.1093/molbev/msj106
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005 |
Lineage-Specific Expansions and Contractions of the Bitter Taste Receptor Gene Repertoire in Vertebrates
Department of Biosystems Science, The Graduate University for Advanced Studies (Sokendai), Hayama, Kanagawa, Japan
E-mail: go_yasuhiro{at}soken.ac.jp.
| Abstract |
|---|
|
|
|---|
The sense of bitter taste plays a critical role in how organisms avoid generally bitter toxic and harmful substances. Previous studies revealed that there were 25 intact bitter taste receptor (T2R) genes in humans and 34 in mice. However, because the recent chicken genome project reported only three T2R genes, it appears that extensive gene expansions occurred in the lineage leading to mammals or extensive gene contractions occurred in the lineage leading to birds. Here, I examined the T2R gene repertoire in placental mammals (dogs, Canis familiaris; and cows, Bos taurus), marsupials (opossums, Monodelphis domestica), amphibians (frogs, Xenopus tropicalis), and fishes (zebrafishes, Danio rerio; and pufferfishes, Takifugu rubripes) to investigate the birth-and-death process of T2R genes throughout vertebrate evolution. I show that (1) the first extensive gene expansions occurred before the divergence of mammals from reptiles/birds but after the divergence of amniotes (reptiles/birds/mammals) from amphibians, (2) subsequent gene expansions continuously took place in the ancestral mammalian lineage and the lineage leading to amphibians, as evidenced by the presence of 15, 18, 26, and 49 intact T2R genes in the dog, cow, opossum, and frog genome, respectively, and (3) contractions of the gene repertoire happened in the lineage leading to chickens. Thus, continuous gene expansions have shaped the T2R repertoire in mammals, but the contractions subsequent to the first round of expansions have made the chicken T2R repertoire narrow. These dramatic changes in the repertoire size might reflect the daily intake of foods from an external environment as a driving force of evolution.
Key Words: bitter taste receptor genes gene expansion gene contraction sensory genes mammal chicken
| Introduction |
|---|
|
|
|---|
Senses of taste, vision, hearing, touch, and olfaction play critical roles in recognizing conditions of the external world and providing organisms with indispensable information for guiding behavior. It is thought that the repertoire of sensory receptors is intimately related to external environments in which organisms have lived and it must have been continuously remolded under ever-changing environments. Sense of taste itself is empirically classified into five basic modalities: sweet, sour, bitterness, salt, and umami (the taste of monosodium glutamate). Out of these, the tastes of sour and salt are perceived by ion channels, mainly proton-gated ion channels in sour (Ugawa et al. 1998
| Materials and Methods |
|---|
|
|
|---|
Samples
The subjects of this study are the dog (C. familiaris), cow (B. Taurus), opossum (M. domestica), chicken (Gallus gallus), frog (X. tropicalis), zebrafish (D. rerio), and pufferfish (T. rubripes). Sequences from each species were retrieved from the Ensembl database (Species-Ensembl version 31): dog (Build 1; 7.6 x genome coverage), cow (Btau_1.0; 3 x coverage), opossum (BROAD0.5; 7.2 x coverage), chicken (Build 1d; 6.6 x coverage), frog (Build 1; 7 x coverage), zebrafish (Build 4b; 5.7 x coverage), and pufferfish (Built 2d; 5.7 x coverage).
Procedure of Gene Identification
To identify the repertoire of T2R genes in each animal, TBlastN searches were conducted on the dog, cow, opossum, and chicken genome databases by using 25 human and 34 mouse intact genes (Conte et al. 2002
, 2003
; Shi et al. 2003
; Go et al. 2005
) as queries with a cutoff E value of 1010. Second, TBlastN searches were carried out to the amphibian and fish by using both 59 mammalian and three chicken intact T2R genes (Hillier et al. 2004
) as queries with a cutoff E value of 105. Among hits obtained from TBlastN searches, those that have more than 270 amino acid sequences with a presumed start and stop codon (average intact T2R genes are
300 amino acids) were selected as candidates for intact T2R genes. The candidate sequences were then used for TBlastN searches on the human genome database (Build 35) and those which showed best hits to the human T2R genes were chosen as intact T2R genes. Subsequently, a third round of TBlastN searches was performed by using the newly identified intact T2R genes from the seven vertebrate species (dogs, cows, opossums, chickens, frogs, zebrafishes, and pufferfishes) as queries to obtain unidentified T2R genes. The procedures for identifying partial T2R genes and T2R pseudogenes are described in the legend of supplementary figure 1 (see Supplementary Material online). In this study, partial T2R genes, which lack either a presumed start codon, a stop codon, or both in the amino acid sequence, were assumed to be intact. All nucleotide T2R sequences used in this study are presented in supplementary figure 2 (see Supplementary Material online).
Data Analyses
The T2R sequences were aligned by ClustalX (Thompson et al. 1997
) and the resulting alignment was manually adjusted. Phylogenetic trees were reconstructed by the Neighbor-Joining (NJ) (Saitou and Nei 1987
) methods implemented in MEGA3 (Kumar, Tamura, and Nei 2004
) with protein-Poisson distances (Nei and Kumar 2000
). Each node of the phylogenetic trees was evaluated by 500 bootstrap replications (Felsenstein 1985
). To date the pseudogenization events, all pseudogenes were also subjected to phylogenetic analysis by inferring the codon frames in the pseudogenes from those of closely related functional genes and conceptually translating to the amino acid sequence. To clarify the origin of each T2R genes under the assumption of molecular clock, a linearized tree was reconstructed after eliminating the T2R genes that evolve significantly faster or slower than the average rate (Takezaki, Rzhetsky, and Nei 1995
).
| Results and Discussion |
|---|
|
|
|---|
T2R Genes in Mammals and Chickens
TBlastN searches were performed on the dog, cow, and opossum genome databases using 25 human and 34 mouse intact T2R genes as queries with a cutoff E value 1010 (see also Materials and Methods and supplementary fig. 1). The results indicate 15 intact genes and four pseudogenes in dogs, 16 and 13 in cows, and 26 and 8 in opossums (table 1). The nucleotide length of intact mammalian T2R genes is from 876 to 1050 bp with the average being 933 bp and the amino acid sequence identity ranging from 19.1% to 89.1%. Due to incompleteness of the cow genome sequences, the total number of T2R genes in cows must be regarded as tentative. However, based on a comparison of syntenic regions among humans, mice, and dogs, the small number of dog T2R genes suggested that there must have been few lineage-specific expansions in this species (data not shown). Although gene conversion is a well-known phenomenon and plays an important role in homogenizing sequences in the genes from multigene families, no identifiable evidence was found for gene conversion, at least in mammalian T2R genes. After publication of the chicken draft genome sequence, TBlastN searches were again performed on those chicken sequences using 59 T2R genes in humans and mice and three in chickens as queries. The results showed that there are three intact T2R genes in chickens, which is consistent with the report from the previous chicken genome project (Hillier et al. 2004
|
Phylogeny of T2R Genes in Mammals and Chickens
To clarify the evolutionary history and relationships among mammalian and chicken T2R genes, the NJ tree was reconstructed on the basis of amino acid sequences of all T2R genes (36, human; 41, mouse; 19, dog; 27, cow; 33, opossum; 3, chicken) (fig. 1). The codon frame of two T2R pseudogenes in cows (Bota-T2R71p and -T2R72p) and one pseudogene in opossums (Modo-T2R32p) could not be determined due to their high sequence divergence; these three sequences were therefore excluded from the phylogenetic analysis. The evolutionary relationships are divided into two categories with respect to orthology. Category A (black lines in fig. 1) is defined as the orthologous gene set that has no extra copies in any species examined in this study. On the other hand, category B (gray lines in fig. 1) is defined such that at least one mammalian species has multiple copies that have resulted from lineage-specific gene duplications. The orthologous relationships of each T2R gene are summarized in table 2. Most of the genes (11/13) that are classified into the category of holding one-to-one orthologous relationships in the previous human and mouse comparison (category C in Go et al. 2005
|
|
The phylogenetic positions of three chicken T2R genes (Gaga-T2R1, -T2R2, and -T2R3) are depicted with black letters in figure 1. The Gaga-T2R1 gene has a monophyletic relationship with the T2R39 and T2R40 gene groups (fig. 1) and it is thought that Gaga-T2R1 is an orthologous gene to the T2R39 and T2R40 gene groups. On the other hand, given the phylogenetic topology and genetic distance, it is possible that neither Gaga-T2R2 nor -T2R3 have any orthologous gene in mammals. Also, as a logical outcome, all mammalian T2R genes except for the T2R39 and T2R40 gene groups do not show any orthologous relationships to chicken T2R genes. These results suggest that the T2R repertoire is quite different between mammals and chickens.
T2R genes appear to result from gene expansions in mammals (Hillier et al. 2004
). However, this assumption should be carefully examined by evolutionary relationships between mammal and chicken T2R genes. To show the origin and orthologous relationships of each T2R gene, the linearized tree was reconstructed after eliminating the T2R genes that evolve significantly faster or slower than the average rate (fig. 2). Among human and chicken intact T2R genes, four T2R genes (Hosa-T2R3, Hosa-T2R7, Hosa-T2R9, and Gaga-T2R2) showed a deviation from the average rate at a 5% significance level. As mentioned earlier, because Gaga-T2R1, and Hosa-T2R39 and -T2R40 probably have an orthologous relationship, the divergence point or speciation point between humans and chickens can be estimated at the dotted line in the linearized tree (fig. 2). According to the linearized tree, the assumption of gene expansions in mammalian lineages is likely to be true for the case of the lineage-specific gene duplications (those genes are classified into category B in fig. 1). However, several gene duplications (white diamonds in fig. 2) occurred before the divergence between mammals and birds/reptiles (dotted line in fig. 2). In other words, many orthologous gene sets already existed
310 MYA at the time of divergence between mammals and birds/reptiles (Benton 1997
; Hedges 2002
) and the number of pre-existing orthologous gene sets before the mammal and bird/reptile divergence can be estimated as 10
11 from the linearized tree (black circles in fig. 2). It can therefore be concluded that chickens have lost many orthologs to mammal T2R genes. In summary, the repertoire of mammalian T2R genes has been formed through continuous gene expansions that have started before the divergence from birds/reptiles (
310 MYA). As for chicken T2R genes, gene contractions subsequent to the first round of gene expansions appear to have made the chicken T2R repertoire narrow.
|
T2R Genes in Amphibians
By using 62 intact T2R genes detected in humans, mice, and chickens as queries, I identified 49 intact T2R genes and 5 pseudogenes in frogs. All of these were confirmed as T2R genes by reciprocal best-hit blast analyses to the human genome. The frog T2R amino acids sequences have similarity from 13.9% to 99.5% and the lengths range from 870 to 1008 bp with the average being 919 bp. The phylogenetic tree of 25 human, 34 mouse, 3 chicken, and 49 frog intact T2R genes shows that mammals and amphibians have different repertoires of T2R genes, indicating independent origins of mammal and amphibian T2R genes (fig. 3).
|
T2R Genes in Fishes
TBlastN searches detected four putative T2R genes in zebrafishes and three in pufferfishes. The nucleotide length and amino acid sequence similarity among fish T2R genes range from 918 to 963 bp and from 12.8% to 79.1%, respectively. Due to close relationships and similar sequences between T2R and V1R genes, one of the two types of pheromone receptor gene families, it would be difficult to annotate the candidate sequences as T2R genes with certainty. This is also consistent with the low bootstrap value for supporting the monophyletic relationship of mammals and fish T2R genes in the phylogenetic tree constructed from vertebrate T2R and V1R genes (fig. 4). However, it should be noted that both zebrafishes and pufferfishes are known to have only one copy of the V1R gene (Pfister and Rodriguez 2005
|
History of T2R Genes in Vertebrates
Using the results obtained from our survey of vertebrate genome sequences, the history of T2R genes is reconstructed and summarized in figure 5. According to the phylogenetic analysis in figure 4, there might be only a few ancestral T2R genes at the incipient stage of vertebrate evolution (
450 MYA), and through the evolutionary history in the lineage leading to the fish, this number of T2R genes did not change very much. About
360 MYA, at the divergence time between amphibians and amniotes (Hedges 2002
50 (fig. 3). Also, the linearized tree showed that the T2R repertoire has extensively expanded to harbor some 10 genes in the lineage of the common ancestors of amniotes, from 360 MYA to 310 MYA (fig. 2). After this initial expansion, however, the T2R repertoire contracted in the chicken lineage and the number of genes decreased to only three in the present chicken genome. In the lineage leading to mammals, on the contrary, the number of genes continuously increased (the numbers of ancestral T2R genes at each node were estimated by the linearized trees as shown in the human and chicken comparison), and after the mammalian radiation, lineage-specific expansions have extensively taken place in each mammalian species, leading to 3040 genes in each species, except for dogs. Most members of the order Carnivora including dogs are either carnivorous or omnivorous (Ray et al. 2005
|
Why Have Chickens Lost T2R Genes?
Although it is almost infeasible to show direct evidence for the cause of T2R gene contractions in chickens, it is tempting to speculate that the contraction is related to genome size reductions in chickens. To know the history of genome size transitions and determine the time of genome size reduction in the lineage leading to chickens, I collected vertebrate genome size data using the available database (Animal Genome Size Database available at http://www.genomesize.com/) and calculated the average genome size in each taxon. The results are summarized in table 3. While terrestrial vertebrates except birds have a genome size of roughly 3 pg (total amount of DNA in the haploid genome, also known as the C-value), the genome size of most birds is less than 2 pg except for ostriches (2.16 pg) with the average being 1.43 ± 0.17 pg from 243 avian species. Specifically, because crocodilians and perhaps also turtles are thought to be the closest living relatives to birds (Janke and Arnason 1997
|
Although we cannot address the issue precisely due to the incompleteness of reptile genome sequences, if genome size reduction is one cause for reduction of the T2R repertoire in chickens, there should be several or more than 10 T2R genes in the reptile genome. The T2R repertoire in reptiles should therefore be examined in order to elucidate further the evolutionary histories of vertebrate T2R genes.
| Conclusion |
|---|
|
|
|---|
In this study, it was revealed that great variation exists among the T2R repertoires of vertebrates, from a few T2R genes in chickens and fishes to 20
50 genes in mammals and amphibians. To comprehend the history of gene expansions and contractions (birth-and-death mode of evolution), phylogenetic analyses were performed and showed only one chicken T2R ortholog in mammals and no apparent mammalian orthologs in amphibians and fishes. After a precise analysis of the tempo and mode of evolution for each T2R gene, I concluded that continuous gene expansions have shaped the T2R repertoire in mammals, but gene contractions subsequent to the first expansions have made the chicken T2R repertoire narrow. These dramatic repertoire changes might reflect different evolutionary processes and different selective pressures in each organism in concert with ever-changing environments. | Supplementary Materials |
|---|
|
|
|---|
The sequences reported in this paper have been deposited in the DNA Data Bank of Japan/EMBL/GenBank database (accession numbers AB249684AB249829). Supplementary figures 1 and 2 are available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/).
| Acknowledgements |
|---|
|
|
|---|
The author thanks Naoyuki Takahata, Yoko Satta, and John Novembre for numerous discussions and helpful comments on the manuscript. This study was supported by the Japanese Society for Promotion of Science (JSPS) grant no. 17-02667 (to Y. G.).
| Footnotes |
|---|
Marta Wayne, Associate Editor
| References |
|---|
|
|
|---|
Adler, E., M. A. Hoon, K. L. Mueller, J. Chandrashekar, N. J. Ryba, and C. S. Zuker. 2000. A novel family of mammalian taste receptors. Cell 100:693702.[CrossRef][Web of Science][Medline]
Benton, M. J. 1997. Vertebrate paleontology. Chapman and Hall, New York.
Chandrashekar, J., K. L. Mueller, M. A. Hoon, E. Adler, L. Feng, W. Guo, C. S. Zuker, and N. J. Ryba. 2000. T2Rs function as bitter taste receptors. Cell 100:703711.[CrossRef][Web of Science][Medline]
Conte, C., M. Ebeling, A. Marcuz, P. Nef, and P. J. Andres-Barquin. 2002. Identification and characterization of human taste receptor genes belonging to the TAS2R family. Cytogenet. Genome Res. 98:4553.[CrossRef][Web of Science][Medline]
. 2003. Evolutionary relationships of the Tas2r receptor gene families in mouse and human. Physiol. Genomics 14:7382.
Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783791.[CrossRef][Web of Science]
Fischer, A., Y. Gilad, O. Man, and S. Pääbo. 2005. Evolution of bitter taste receptors in humans and apes. Mol. Biol. Evol. 22:432436.
Gilad, Y., V. Wiebe, M. Przeworski, D. Lancet, and S. Pääbo. 2004. Loss of olfactory receptor genes coincides with the acquisition of full trichromatic vision in primates. PLoS Biol. 2:120125.
Go, Y., Y. Satta, O. Takenaka, and N. Takahata. 2005. Lineage-specific loss of function of bitter taste receptor genes in humans and nonhuman primates. Genetics 170:313326.
Grus, W. E., P. Shi, Y. P. Zhang, and J. Zhang. 2005. Dramatic variation of the vomeronasal pheromone receptor gene repertoire among five orders of placental and marsupial mammals. Proc. Natl. Acad. Sci. USA 102:57675772.
Grus, W. E., and J. Zhang. 2004. Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats. Gene 340:303312.[CrossRef][Web of Science][Medline]
Hedges, S. B. 2002. The origin and evolution of model organisms. Nat. Rev. Genet. 3:838849.[CrossRef][Web of Science][Medline]
Hillier, L. W., W. Miller, E. Birney et al. (176 co-authors). 2004. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature 432:695716.[CrossRef][Medline]
Janke, A., and U. Arnason. 1997. The complete mitochondrial genome of Alligator mississippiensis and the separation between recent archosauria (birds and crocodiles). Mol. Biol. Evol. 14:12661272.[Abstract]
Kretz, O., P. Barbry, R. Bock, and B. Lindemann. 1999. Differential expression of RNA and protein of the three pore-forming subunits of the amiloride-sensitive epithelial sodium channel in taste buds of the rat. J. Histochem. Cytochem. 47:5164.
Kumar, S., K. Tamura, and M. Nei. 2004. MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief. Bioinform. 5:150163.
Kumazawa, Y., and M. Nishida. 1999. Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink: statistical evidence for archosaurian affinity of turtles. Mol. Biol. Evol. 16:784792.[Abstract]
Li, X., L. Staszewski, H. Xu, K. Durick, M. Zoller, and E. Adler. 2002. Human receptors for sweet and umami taste. Proc. Natl. Acad. Sci. USA 99:46924696.
Lin, W., T. E. Finger, B. C. Rossier, and S. C. Kinnamon. 1999. Epithelial Na+ channel subunits in rat taste cells: localization and regulation by aldosterone. J. Comp. Neurol. 405:406420.[CrossRef][Web of Science][Medline]
Matsunami, H., J. P. Montmayeur, and L. B. Buck. 2000. A family of candidate taste receptors in human and mouse. Nature 404:601614.[CrossRef][Medline]
Nei, M., X. Gu, and T. Sitnikova. 1997. Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc. Natl. Acad. Sci. USA 94:77997806.
Nei, M., and S. Kumar. 2000. Molecular evolution and phylogenetics. Oxford University Press, New York.
Nelson, G., M. A. Hoon, J. Chandrashekar, Y. Zhang, N. J. Ryba, and C. S. Zuker. 2001. Mammalian sweet taste receptors. Cell 106:381390.[CrossRef][Web of Science][Medline]
Niimura, Y., and M. Nei. 2005. Evolutionary dynamics of olfactory receptor genes in fishes and tetrapods. Proc. Natl. Acad. Sci. USA 102:60396044.
Pfister, P., and I. Rodriguez. 2005. Olfactory expression of a single and highly variable V1r pheromone receptor-like gene in fish species. Proc. Natl. Acad. Sci. USA 102:54895494.
Ray, J., K. Redford, R. Steneck, and J. Berger. 2005. Large carnivores and the conservation of biodiversity. Island Press, Chicago.
Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406425.[Abstract]
Shi, P., J. Zhang, H. Yang, and Y. P. Zhang. 2003. Adaptive diversification of bitter taste receptor genes in mammalian evolution. Mol. Biol. Evol. 20:805814.
Stevens, D. R., R. Seifert, B. Bufe, F. Muller, E. Kremmer, R. Gauss, W. Meyerhof, U. B. Kaupp, and B. Lindemann. 2001. Hyperpolarization-activated channels HCN1 and HCN4 mediate responses to sour stimuli. Nature 413:631635.[CrossRef][Medline]
Takezaki, N., A. Rzhetsky, and M. Nei. 1995. Phylogenetic test of the molecular clock and linearized trees. Mol. Biol. Evol. 12:823833.[Abstract]
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:48764882.
Ugawa, S., Y. Minami, W. Guo, Y. Saishin, K. Takatsuji, T. Yamamoto, M. Tohyama, and S. Shimaria. 1998. Receptor that leaves a sour taste in the mouth. Nature 395:555556.[CrossRef][Medline]
van Tuinen, M., and S. B. Hedges. 2001. Calibration of avian molecular clocks. Mol. Biol. Evol. 18:206213.
Wang, X., S. D. Thomas, and J. Zhang. 2004. Relaxation of selective constraint and loss of function in the evolution of human bitter taste receptor genes. Hum. Mol. Genet. 13:26712678.
Wong, G. T., K. S. Gannon, and R. F. Margolskee. 1996. Transduction of bitter and sweet taste by gustducin. Nature 381:796800.[CrossRef][Medline]
Young, J. M., C. Friedman, E. M. Williams, J. A. Ross, L. Tonnes-Priddy, and B. J. Trask. 2002. Different evolutionary processes shaped the mouse and human olfactory receptor gene families. Hum. Mol. Genet. 11:535546.
Young, J. M., M. Kambere, B. J. Trask, and R. P. Lane. 2005. Divergent V1R repertoires in five species: amplification in rodents, decimation in primates, and a surprisingly small repertoire in dogs. Genome Res. 15:231240.
Zhang, J., and D. M. Webb. 2003. Evolutionary deterioration of the vomeronasal pheromone transduction pathway in catarrhine primates. Proc. Natl. Acad. Sci. USA 100:83378341.
Zhang, X., and S. Firestein. 2002. The olfactory receptor gene superfamily of the mouse. Nat. Neurosci. 5:124133.[Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
X. Li, D. Glaser, W. Li, W. E. Johnson, S. J. O'Brien, G. K. Beauchamp, and J. G. Brand Analyses of Sweet Receptor Gene (Tas1r2) and Preference for Sweet Stimuli in Species of Carnivora J. Hered., July 1, 2009; 100(suppl_1): S90 - S100. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Oike, T. Nagai, A. Furuyama, S. Okada, Y. Aihara, Y. Ishimaru, T. Marui, I. Matsumoto, T. Misaka, and K. Abe Characterization of Ligands for Fish Taste Receptors J. Neurosci., May 23, 2007; 27(21): 5584 - 5592. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||






