- Research article
- Open Access
Variation in ligand responses of the bitter taste receptors TAS2R1 and TAS2R4 among New World monkeys
© The Author(s). 2016
- Received: 4 July 2016
- Accepted: 30 September 2016
- Published: 12 October 2016
New World monkeys (NWMs) are unique in that they exhibit remarkable interspecific variation in color vision and feeding behavior, making them an excellent model for studying sensory ecology. However, it is largely unknown whether non-visual senses co-vary with feeding ecology, especially gustation, which is expected to be indispensable in food selection. Bitter taste, which is mediated by bitter taste receptors (TAS2Rs) in the tongue, helps organisms avoid ingesting potentially toxic substances in food. In this study, we compared the ligand sensitivities of the TAS2Rs of five species of NWMs by heterologous expression in HEK293T cells and calcium imaging.
We found that TAS2R1 and TAS2R4 orthologs differ in sensitivity among the NWM species for colchicine and camphor, respectively. We then reconstructed the ancestral receptors of NWM TAS2R1 and TAS2R4, measured the evolutionary shift in ligand sensitivity, and identified the amino acid replacement at residue 62 as responsible for the high sensitivity of marmoset TAS2R4 to colchicine.
Our results provide a basis for understanding the differences in feeding ecology among NWMs with respect to bitter taste.
- Bitter taste receptor
- G protein-coupled receptor
- New World monkey
- Interspecific functional variation
- Molecular evolution
Among the five basic taste qualities (sweet, bitter, sour, salty, and umami), bitter taste is thought to be particularly important for the survival of animals because it contributes to the avoidance of potentially toxic substances in food. Bitter taste is mediated by bitter taste receptors (TAS2Rs) [1, 2], which are mainly expressed in taste buds in the tongue. The repertoire of TAS2Rs in the genome varies considerably among animal species . TAS2Rs are G protein-coupled receptors that recognize a wide variety of bitter substances as ligands.
The responses of human TAS2Rs to various bitter substances have been thoroughly investigated . However, little is known about variation in the response properties of TAS2Rs among species. It is reasonable to speculate that the properties of bitter taste reflect adaptation to species-specific feeding environments. Variation in the response sensitivity of TAS2R16, a well-studied TAS2R, to salicin and its derivatives has been observed among humans, chimpanzees, and macaques [5, 6]. However, little is known about interspecific functional variation in TAS2R orthologs other than TAS2R16 and 38 in primates.
Among TAS2Rs, TAS2R1 and TAS2R4 are promising targets for interspecific comparative studies because they are relatively highly conserved among a wide range of species  and some functional studies have been performed [7–9]. In this study, we investigated interspecific variation in TAS2R1 and TAS2R4 in primates.
New World monkeys (NWMs) exhibit remarkable interspecific variation in color vision and diet (including fruits, nuts, leaves, flowers, plant exudates, and insects) [7–9], making them an excellent model for studying the relationships among senses . For example, owl monkeys are nocturnal, lack cone-based color vision , and have primarily frugivorous diets supplemented with insects, flowers, nectar, and leaves . Marmosets are polymorphic with respect to color vision (dichromacy in males and some females and trichromacy in other females)  and some species (e.g., Callithrix jacchus) have unique exudativore (gum-feeder) diets supplemented with insects .
In the present study, we performed calcium assays to compare the responses of NWM TAS2R1 and TAS2R4 to various bitter substances among species. We found that owl monkey TAS2R1 has a high sensitivity to camphor and marmoset TAS2R4 has a high sensitivity to colchicine. Furthermore, we inferred the amino acid sequences of a set of ancestral receptors of NWM TAS2R1 and TAS2R4. A functional assay of the ancestral receptors revealed the evolutionary timing of the shift in the ligand sensitivity of NWM TAS2R1 and TAS2R4. We also identified candidate amino acid residues responsible for the evolutionary shifts in ligand sensitivity. Furthermore, using site-directed mutagenesis experiments, we demonstrated that one of the candidate residues, F62, is responsible for the high sensitivity of marmoset TAS2R4 to colchicine. These results suggest a possible mechanism by which NWMs adapt to species-specific feeding environments by modifying the sensitivity of bitter taste receptors.
Construction of TAS2R expression vectors
Genomic DNA derived from blood samples of a captive common marmoset (Callithrix jacchus) and a captive Azara’s owl monkey (Aotus azarae) were provided by the Cooperative Research Program of the Primate Research Institute of Kyoto University under their Guidelines for the Care and Use of Laboratory Primates. Feces-derived genomic DNA of a wild white-faced capuchin (Cebus capucinus), a wild black-handed spider monkey (Ateles geoffroyi), and a wild mantled howler (Alouatta palliata) were obtained in previous studies by SK [11, 12]. TAS2R1 and TAS2R4 genes were cloned from these genome DNA samples using anthropoid genome-based sequences as PCR and sequencing primers .
TAS2R1 and TAS2R4 of marmosets, capuchin monkeys, owl monkeys, spider monkeys, and howler monkeys were tagged with the first 45 amino acids of rat somatostatin receptor 3 (sst tag)  and the last 8 amino acids of bovine rhodopsin (1D4 tag)  at the N- and C-terminal ends, respectively, and cloned into the mammalian expression vector pEAK10 (Edge BioSystems, Inc., Gaithersburg, MD, USA).
Transfection of HEK293T cells
HEK293T cells were transfected with the TAS2R expression constructs together with the chimeric G protein subunit Gα16gust44  in the mammalian expression vector pcDNA3.1 (Life Technologies, Inc., Carlsbad, CA, USA) using Lipofectamine 2000 (Life Technologies, Inc.). The empty pEAK10 vector was used as a negative control.
Calcium assays were performed as previously reported [4, 5], with some modifications. Briefly, one day after transfection, culture medium was exchanged with the assay buffer (10 mM HEPES, 130 mM NaCl, 10 mM glucose, 5 mM KCl, 2 mM CaCl2, 1.2 mM MgCl2, pH 7.4) and loaded with a calcium indicator, Calcium 4 (Molecular Devices, Inc., Eugene, OR, USA). Ligand solutions were prepared by dissolving bitter substances (camphor or colchicine) in the assay buffer and were added to cells. Fluorescence was measured at 525 nm following excitation at 485 nm with the FlexStation 3 Microplate Reader (Molecular Devices, Inc.). The calcium response amplitudes were expressed as ΔF/F, which is the ratio of the ligand-dependent increase in fluorescence to the fluorescence before the ligand addition. The response of cells that were transfected with the empty pEAK10 vector and Gα16gust44 was defined as the mock response (TAS2R-independent response) and subtracted from all responses.
Ancestral sequence reconstruction
Ancestral amino acid sequences of NWM TAS2R1 and TAS2R4 were inferred using the maximum likelihood method implemented in MEGA5 [16, 17] with the Dayhoff and JTT [18, 19] amino acid substitution models . The reconstructed ancestral TAS2Rs were tagged with sst and 1D4 tags at the N- and C-terminal regions, respectively, cloned into the pEAK10 vector, and expressed in HEK293T cells. Calcium assays were performed following the same protocol as that used for extant NWM TAS2Rs.
Site-directed mutagenesis was performed using the QuikChange Multisite-Directed Mutagenesis Kit (Agilent Technologies, Inc., Santa Clara, CA, USA).
Single nucleotide polymorphisms were detected in TAS2R1of the capuchin monkey, and the spider monkey samples, and in TAS2R4 of the mantled howler monkey sample (Additional file 1: Figure S1). Because there were no significant difference between haplotypes subjected to calcium imaging assays to explore the intraspecific variation in ligand sensitivity, we used representative haplotypes for each species to explore the interspecific variation.
All of the ancestral TAS2R4 receptors exhibited relatively low sensitivity to colchicine, comparable to that of TAS2R4 of the capuchin monkey, spider monkey, and howler monkey. However, they had significantly lower sensitivities than that of marmoset TAS2R4. Specifically, “ancestor 2” had a significantly lower sensitivity than that of marmoset TAS2R4, indicating that the high sensitivity of TAS2R4 to colchicine was acquired during evolution from “ancestor 2” to the marmoset. The sensitivity of “ancestor 1” was significantly higher than that of the owl monkey TAS2R4, indicating that the sensitivity was reduced during the evolution from “ancestor 1” to the owl monkey (Fig. 2b).
Candidate amino acid residues responsible for the evolutionary shifts in ligand sensitivity of NWM TAS2R1 and TAS2R4 can be detected. Only one amino acid residue differed between “ancestor 2” and “ancestor 1” of TAS2R4, position 82, which was occupied by a cysteine in “ancestor 2” and replaced with phenylalanine in “ancestor 1” (Fig. 2a). This indicates that the C82F replacement is responsible for the increase in ligand sensitivity from “ancestor 2” to “ancestor 1.” When we compared the amino acid sequence of the marmoset TAS2R1 with that of “ancestor 1,” we detected 11 substitutions (I8L, I15M, F22S, L95F, V149A, H164R, N158K, I174A, Y248C, I258T, S271C), one or some of which may be responsible for the reduction in the ligand sensitivity of TAS2R1 from “ancestor 1” to the marmoset. Similarly, the candidate amino acid residues responsible for the shifts in the ligand sensitivity of TAS2R4 can also be detected; L7F, S24N, L62F, T143N, V156F, F189L, Q208R, H249N, or M257L may explain the increase from “ancestor 2” to the marmoset and A9G, V70I, A135T, T153A, V156F, L177V, R205K, or F253T may explain the reduction from “ancestor 1” to the owl monkey.
We measured the responses of NWM TAS2R1 and TAS2R4 receptors to bitter substances to explore variation in ligand sensitivity among extant species and reconstructed ancestors. We detected variation in the sensitivity of TAS2R1 and TAS2R4 to camphor or colchicine, respectively. TAS2R1 of the owl monkey and TAS2R4 of the marmoset exhibited significantly higher response amplitudes than the corresponding TAS2R1 and TAS2R4 receptors in other species (Fig. 1). Importantly, our results demonstrate how an ancestral receptor that is intermediately responsive to a variety of compounds can be shaped by lineage-specific evolution to respond more strongly to different compounds in different extant taxa.
The ecological implications of the variation in ligand sensitivity among species suggest that functional variation in compounds of dietary items or antifeedants are shaping the taste receptor genes of primates. Our findings emphasize the importance of investigating the role of bitter substances in the diets of NWMs. Recent studies have suggested the presence of TAS2Rs not only in the tongue, but also in other organs, like the nasal cavity and airway . These receptors might therefore have a role in the detection of chemical compounds like camphor and colchicine based on their bitter taste, but also via other senses, such as olfaction.
Based on an analysis of ancestral receptors, we identified candidate amino acid substitutions responsible for the evolutionary shift in ligand sensitivity among species (Fig. 2). Specifically, we detected the following substitutions: C82F for the increase between “Ancestor 2” and “Ancestor 1” of TAS2R1; I8T, A13S, H144R, S145T, M157L, S161P, I166V, F177S, S178F, A181T, V192L, V194L, T205A, S249C, and I278F for the increase between “Ancestor 1” and owl monkey TAS2R1; I8L, I15M, F22S, L95F, V149A, H164R, N158K, I174A, Y248C, I258T, and S271C for the decrease between “Ancestor 1” and marmoset TAS2R1; L7F, S24N, L62F, T143N, V156F, F189L, Q208R, H249N, and M257L for the increase between “Ancestor 2” and marmoset TAS2R4; A9G, V70I, A135T, T153A, V156F, L177V, R205K, and F253T for the decrease between “Ancestor 1” and owl monkey TAS2R4. Homology modeling based on rhodopsin structures (or a recently resolved smoothened structure ) as a template would help reveal the mechanism underlying the interspecific variation. Additionally, site-directed mutagenesis experiments for these candidate residues are needed to identify the precise amino acid residues responsible for the observed interspecific variation in ligand sensitivity. Towards this goal, we conducted site-directed mutagenesis experiments and found that the amino acid difference at position 62 is one of the responsible sites for the difference in response amplitude between TAS2R4 of marmoset and other species. Because amino acid residue at position 62 would be situated in the extracellular side of second transmembrane domain , it might affect the ligand binding reaction of TAS2R4. Further investigations of additional candidate mutations may explain receptor function differences between NWMs.
In summary, we compared the ligand sensitivity of the bitter taste receptors of NWMs. We found that sensitivity of TAS2Rs varies among species. These results raise the possibility that variation among TAS2Rs in ligand sensitivity is associated with ecological adaptations. Using an ancestral analysis, we detected evolutionary shifts in ligand sensitivity during the molecular evolution of TAS2R1 and TAS2R4, and we identified amino acid residues responsible for the shifts. Most of the ancestral receptors showed intermediate responsiveness between the most and least responsive existing receptors (i.e., TAS2R1: owl monkey and other species, respectively; TAS2R4: marmoset and owl monkey, respectively). Interestingly, these results suggest that ancestral receptors, characterized by intermediate responsiveness, evolved lower or higher responsiveness over time, as shown in the extant study species. Molecular and environmental studies will provide the basis for elucidating these phenomena and, more broadly, the diversification of NWMs, from molecular and ecological perspectives.
We thank Dr. T. Ueda for providing the G16/gust44 expression vector.
This work was financially supported by the Leading Graduate Program in Primatology and Wildlife Sciences, Kyoto University, KAKENHI grants from the Japan Society for the Promotion of Science (Nos. 25891012 to KT, 12 J04270 to TH, 22247036 and 15H02421 to SK, 24370096, 24405018, 26117512, 15H05242 to HI) and the Yamada Science Foundation to H. I.
Availability of data and materials
Accession numbers of each haplotypes are shown as follows: CjT2R1 (LC184478), CcT2R1-1 (LC184479), CcT2R1-2 (LC184480), AaT2R1-1(LC184481), AaT2R1-2(LC184482), AgT2R1-1(LC184483), AgT2R1-2 (LC184484), ApT2R1 (LC184485), CjT2R4 (LC184486), CcT2R4 (LC184487), AaT2R4 (LC184488), AgT2R4 (LC184489), ApT2R4-1 (LC184490), and ApT2R4-2 (LC184491).
Original statistical data are shown in Additional file 3: Table S1. The fluorescence measurement data are deposited in Dryad. Data available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.6q4g6 .
KT carried out acquisition, analysis and interpretation of data, drafted the manuscript. MO and KS participated to the molecular genetic studies and acquisition of the data. NSH participated in the sequence alignment. TH participated in the sequence alignment, analysis of data, and revising the article. TM and YI contributed to acquisition of data. FA and ADM participated in the providing samples and revising the article. SK and HI conceived and designed the study, and finalized the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Blood samples of a common marmoset (Callithrix jacchus) and an Azara’s owl monkey (Aotus azarae) were taken from captive individuals housed at the Primate Research Institute, with the approval from the Animal Ethics Committee of the Primate Research Institute, Kyoto University (Permit Numbers: 2014-045 and 2015-082) under their Guidelines for the Care and Use of Laboratory Primates. The title of the project is “Genotyping of various primate species” which include the genotyping of the taste receptor genes and functional analysis of the proteins.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Bachmanov AA, Beauchamp GK. Taste receptor genes. Annu Rev Nutr. 2007;27:389–414.View ArticlePubMedPubMed CentralGoogle Scholar
- Behrens M, Meyerhof W. Gustatory and extragustatory functions of mammalian taste receptors. Physiol Behav. 2011;105(1):4–13.View ArticlePubMedGoogle Scholar
- Hayakawa T, Suzuki-Hashido N, Matsui A, Go Y. Frequent expansions of the bitter taste receptor gene repertoire during evolution of mammals in the euarchontoglires clade. Mol Biol Evol. 2014;31(8):2018–31.View ArticlePubMedGoogle Scholar
- Meyerhof W, Batram C, Kuhn C, Brockhoff A, Chudoba E, Bufe B, Appendino G, Behrens M. The molecular receptive ranges of human TAS2R bitter taste receptors. Chem Senses. 2010;35(2):157–70.View ArticlePubMedGoogle Scholar
- Imai H, Suzuki N, Ishimaru Y, Sakurai T, Yin L, Pan W, Abe K, Misaka T, Hirai H. Functional diversity of bitter taste receptor TAS2R16 in primates. Biol Lett. 2012;8(4):652–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandau MM, Goodman JR, Thomas A, Rucker JB, Rawson NE. A functional comparison of the domestic cat bitter receptors Tas2r38 and Tas2r43 with their human orthologs. BMC Neurosci. 2015;16:33.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawamura S, Hiramatsu C, Melin AD, Schaffner CM, Aureli F, Fedigan LM. Polymorphic color vision in primates: evolutionary considerations. In: Hirai H, Imai H, Go Y, editors. Post-genome biology of primates. Tokyo: Springer; 2012. p. 93–120.View ArticleGoogle Scholar
- Levenson DH, Fernandez-duque E, Evans S, Jacobs GH. Mutational changes in S-cone opsin genes common to both nocturnal and cathemeral aotus monkeys. Am J Primatol. 2007;69(7):757–65.View ArticlePubMedGoogle Scholar
- Kinzey WG. New World primates: ecology, evolution, and behavior. New York: Aldine de Gruyte; 1997.Google Scholar
- Jacobs GH. New World monkeys and color. J Primatol. 2007;28(4):729–59.Google Scholar
- Hiwatashi T, Okabe Y, Tsutsui T, Hiramatsu C, Melin AD, Oota H, Schaffner CM, Aureli F, Fedigan LM, Innan H, et al. An explicit signature of balancing selection for color-vision variation in new world monkeys. Mol Biol Evol. 2010;27(2):453–64.View ArticlePubMedGoogle Scholar
- Matsushita Y, Oota H, Welker BJ, Pavelka MS, Kawamura S. Color vision variation as evidenced by hybrid L/M opsin genes in wild populations of trichromatic New World monkeys. Int J Primatol. 2014;35:71–87.View ArticlePubMedGoogle Scholar
- Ammon C, Schafer J, Kreuzer OJ, Meyerhof W. Presence of a plasma membrane targeting sequence in the amino-terminal region of the rat somatostatin receptor 3. Arch Physiol Biochem. 2002;110(1-2):137–45.View ArticlePubMedGoogle Scholar
- Molday RS, MacKenzie D. Monoclonal antibodies to rhodopsin: characterization, cross-reactivity, and application as structural probes. Biochemistry. 1983;22(3):653–60.View ArticlePubMedGoogle Scholar
- Ueda T, Ugawa S, Yamamura H, Imaizumi Y, Shimada S. Functional interaction between T2R taste receptors and G-protein alpha subunits expressed in taste receptor cells. J Neurosci. 2003;23(19):7376–80.PubMedGoogle Scholar
- Nei M, Kumar S. Molecular Evolution and Phylogenetics. New York: Oxford University Press; 2000.Google Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Dayhoff MO, Schwartz RM, Orcutt BC. A model of evolutionary change in proteins. In: Dayhoff MO, editor. Atlas of protein sequence and structure. Washington, DC: National Biomedical Research Foundation; 1978. p. 345–52.Google Scholar
- Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Computer applications in the biosciences. CABIOS. 1992;8(3):275–82.PubMedGoogle Scholar
- Matsumoto Y, Hiramatsu C, Matsushita Y, Ozawa N, Ashino R, Nakata M, Kasagi S, Di Fiore A, Schaffner CM, Aureli F, et al. Evolutionary renovation of L/M opsin polymorphism confers a fruit discrimination advantage to ateline New World monkeys. Mol Ecol. 2014;23(7):1799–812.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider H, Sampaio I. The systematics and evolution of New World primates - a review. Mol Phylogenet Evol. 2015;82(Pt B):348–57.View ArticlePubMedGoogle Scholar
- Wang C, Wu H, Katritch V, Han GW, Huang XP, Liu W, Siu FY, Roth BL, Cherezov V, Stevens RC. Structure of the human smoothened receptor bound to an antitumour agent. Nature. 2013;497:338–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiener A, Shudler M, Levit A, Niv MY. BitterDB: a database of bitter compounds. Nucleic Acids Res. 2012;40(Database issue):D413–9.View ArticlePubMedGoogle Scholar
- Tsutsui K, Otoh M, Sakurai K, Suzuki-Hashido N, Hayakawa T, Misaka T, ISHIMARU Y, Aureli F, Melin A, Kawamura S, Imai H. Data from: Variation in ligand responses of the bitter taste receptors TAS2R1 and TAS2R4 among New World monkeys. Dryad Digital Repository. http://dx.doi.org/10.5061/dryad.6q4g6