Dietary specialization drives multiple independent losses and gains in the bitter taste gene repertoire of Laurasiatherian Mammals
- Zhijin Liu†1, 2Email author,
- Guangjian Liu†1,
- Frank Hailer†2,
- Pablo Orozco-terWengel†2,
- Xinxin Tan†1, 3,
- Jundong Tian4,
- Zhongze Yan1, 3,
- Baowei Zhang5 and
- Ming Li1Email author
© The Author(s). 2016
Received: 28 April 2016
Accepted: 22 June 2016
Published: 29 June 2016
Bitter taste perception is essential for species with selective food intake, enabling them to avoid unpalatable or toxic items. Previous studies noted a marked variation in the number of TAS2R genes among various vertebrate species, but the underlying causes are not well understood. Laurasiatherian mammals have highly diversified dietary niche, showing repeated evolution of specialized feeding preferences in multiple lineages and offering a unique chance to investigate how various feeding niches are associated with copy number variation for bitter taste receptor genes.
Here we investigated the evolutionary trajectories of TAS2Rs and their implications on bitter taste perception in whole-genome assemblies of 41 Laurasiatherian species. The number of intact TAS2Rs copies varied considerably, ranging from 0 to 52. As an extreme example of a narrow dietary niche, the Chinese pangolin possessed the lowest number of intact TAS2Rs (n = 2) among studied terrestrial vertebrates. Marine mammals (cetacea and pinnipedia), which swallow prey whole, presented a reduced copy number of TAS2Rs (n = 0-5). In contrast, independent insectivorous lineages, such as the shrew and insectivorous bats possessed a higher TAS2R diversity (n = 52 and n = 20-32, respectively), exceeding that in herbivores (n = 9-22) and omnivores (n = 18-22).
Besides herbivores, insectivores in Laurasiatheria tend to have more functional TAS2Rs in comparison to carnivores and omnivores. Furthermore, animals swallowing food whole (cetacean, pinnipedia and pangolin) have lost most functional TAS2Rs. These findings provide the most comprehensive view of the bitter taste gene repertoire in Laurasiatherian mammals to date, casting new light on the relationship between losses and gains of TAS2Rs and dietary specialization in mammals.
The ability to perceive taste provides animals with important dietary information that can be crucial for their survival. Bitter taste perception is mediated by a group of G protein-coupled receptors, which are encoded by the type 2 taste receptor gene family (TAS2Rs) . Previous studies noted a marked variation in the number of TAS2R genes among various vertebrate species, but the underlying cause is not well understood [2, 3]. A study surveying bitter taste receptor genes in 52 vertebrate species (fish, birds, amphibians, reptiles and mammals)  reported that the number of TAS2R genes per species correlates with the proportion of plants in their diet. Furthermore, carnivores were found to have fewer TAS2R genes compared to herbivores and omnivores . Because plant tissues contain more toxic compounds than animal tissues, these findings suggest that dietary toxins are a major selective force influencing the diversity of the TAS2R repertoire. Given this finding, it is reasonable to speculate that animals whose diet contains more toxins, besides herbivores, tend to have more functional TAS2Rs, while animals with less dependence on bitter taste are likely to lose a certain number of functional TAS2Rs. To test this hypothesis it is necessary to rely on a taxon with highly diversified and specialized dietary niches.
Over an estimated ~80 million years of diversification , Laurasiatherian mammals have evolved multiple lineage-specific feeding characteristics, such as the obligate ant-ingesting pangolins, the insectivorous shrews, carnivores (cats, dogs and allies), ungulate herbivores, herbivorous or insectivorous bats, and marine mammals that adopted a purely marine diet . During the process of diet differentiation in Laurasiatheria, bitter taste preferences presumably evolved adaptively to avoid ingestion of indigestible or harmful compounds. In particular, there are several independent insectivorous lineages (shrews, insectivorous bats and so on). Eating insects is similar to herbivory, where detecting toxins is important. Many insects release defensive secretions that are toxic to their predators [8–10]. It could therefore be hypothesized that insectivores may be exposed to a similar amount of toxins as herbivores. Furthermore, there are also several independent lineages (cetaceans, pinnipeds and pangolins), which have evolved specialised foraging behaviour, ingesting prey whole. These species are likely to lose a certain number of functional TAS2Rs due to less dependence on the taste perception. In summary, Laurasiatheria offers a unique chance to investigate how various feeding niches are associated with copy number variation for bitter taste receptor genes. Since the publication of Li and Zhang  additional vertebrate genomes have become available for study. Wang and Zhao  have reported a small repertoire of TAS2Rs in birds based on a genomic survey including 48 birds. Furthermore, Hayakawa et al.  have described the frequent expansions of the TAS2Rs in the Euarchontoglires clade, which is the sister clade of Laurasitherian mammals. However, less research has been carried out concerning the evolutionary history of the repertoire of TAS2Rs in Laurasitheria. Here we use a comprehensive set of 41 Laurasitherian genomes encompassing all dietary niches across mammals and test the potential impact of all dietary types (herbivory, insectivory, carnivory and omnivory) and foraging behaviours (swallowing the prey whole or not) on the evolutionary scenarios of TAS2Rs in Laurasiatherian mammals. The results showed that the diverse dietary characteristics in Laurasiatheria have had a larger impact in the TAS2Rs repertoires than previously acknowledged.
TAS2Rs repertoires in Laurasiatheria
Lineage-specific gene degeneration and expansion in Laurasiatherian TAS2Rs
Chinese pangolin (Manis pentadactyla) only contained two intact bitter taste genes (TAS2R1 and TAS2R7; Fig. 1). In Cetartiodactyla only the minke whale (Balaenoptera acutorostrata) and the Antarctic minke whale (Balaenoptera bonaerensis) retain a single TAS2R gene, although each species retains a different one, i.e. TAS2R16 and TAS2R67, respectively. The remaining whales and dolphins lost all intact TAS2Rs that our analysis indicates were present in their common ancestor. Similarly, another group of marine mammals, pinnipedia (walrus and Weddell seal, Carnivora), only present 4–5 TAS2Rs while having lost the rest of them (Fig. 2; Additional file 1: Table S2), which is consistent with the results of a previous genomic analysis on carnivores .
The increase in numbers of intact TAS2Rs in the common shrew was observed and caused by the lineage-specific duplication of TAS2R7, TAS2R8, TAS2R372 and TAS2R408. These TAS2Rs duplications are the outcome of duplications in the genomic region NW_004545902.1 of the common shrew (Fig. 3a and b; Additional file 1: Figure S1). Similarly to the shrew, we found lineage-specific duplication of TAS2R16, TAS2R18 and TAS2R41 in insectivorous bats (Fig. 3a; Additional file 1: Table S2) corresponding to a genomic duplication that increased the gene copy number of TAS2R16 and TAS2R18 in the insectivorous Brandt’s bat (Myotis brandtii). Contrastingly, the frugivorous big brown bat (Eptesicus fuscus) shows no duplication in its orthologous genomic region (Additional File 1: Figure S1).
Correlation between diet and TAS2Rs repertoire
Secondly, to test the impact of insectivory and herbivory on the TAS2Rs diversity separately, we divided the whole dataset into two sub-datasets. The first sub-dataset includes the insectivores, carnivores and omnivores, while excluding herbivores. The second dataset includes the herbivores, carnivores and omnivores, thereby excluding insectivores from the analysis. Subsequently, we conducted the PICs analysis on the two datasets separately. For the first sub-dataset including insectivores, omnivores and carnivores, the PICs in the dietary code were positively correlated with that of intact TAS2Rs numbers (ρ = 0.528, P = 0.017; Fig. 4b), supporting the hypothesis that insectivores exhibit more intact TAS2Rs than omnivores and carnivores. For the second sub-dataset including herbivores, omnivores and carnivores, we performed two PIC analyses. One included all the herbivores and another excluded the tylopods. The PICs analysis was based on the data including all the herbivores revealed a marginal significant correlation between the dietary code and intact TAS2Rs number (ρ = 0.289, P = 0.051). Furthermore, the PICs analysis excluding tylopods supported a significant correlation between the dietary code and intact TAS2Rs number (ρ = 0.505, P = 0.003; Fig. 4c). Camels live in a harsh desert environment and mainly feed on xerophilous plants such as camel thorn (Alhagi spp.) and plants from the Chenopodiaceae family . Similarly, alpaca live in areas of high altitude and extreme cold, with relatively few food sources . Consequently, in light of the available variety of plants in their environments, it is likely that tylopods adapted to a limited range of plant species [17, 18] and gradually decreased dependence on bitter taste and experienced a reduction of TAS2Rs numbers. The reduction of TAS2Rs in tylopods still reflects the changes in dietary preferences, with a correlation between the functional TAS2Rs number and abundance of potential toxins in the diet, and does not refute the general trend that herbivores evolved more copy numbers of intact TAS2Rs than omnivores/carnivores.
Thirdly, we coded a species as 0 for swallowing the food whole and 1 for species that chew their food. Then the PICs analysis was implemented again to test the correlation and regression analysis of the numbers of intact TAS2Rs and the foraging behaviour code for Laurasiatherian. The PICs in the foraging behaviour code also positively correlated with that of intact TAS2Rs numbers (ρ = 0.604, P <0.001, Fig. 4d), supporting the hypothesis that the animals less dependent on bitter taste are likely to pose a lower number of functional TAS2Rs.
TAS2Rs degeneration in pangolin and marine mammals
Pangolins are highly specialized in their diet, feeding almost exclusively on termites and ants that are located by scent . Different pangolin species typically show a clear feeding preference for particular species of ants and termites, often rejecting all other available prey species . Furthermore, pangolins catch living ants with their tongue, and all food is ingested whole and then crushed in the lower section of the stomach. This feeding behaviour together with our finding of very low TAS2R diversity in Chinese pangolin indicates that in this taxon bitter taste perception is highly reduced as a consequence of (a) a hyper-specialization on a very specific feeding resource, and (b) a lack of reliance on bitter taste perception in the identification of prey. Even with no knowledge of the tuning width of the only two TAS2Rs in pangolins, the exclusive diet and special feeding behaviour of pangolins might ultimately have rendered their bitter taste receptors redundant resulting in the loss of most TAS2Rs. On the other hand, it should be interesting to assess the detection range of the pangolins’ only two TAS2Rs on various bitter compounds.
Examples that illustrate the relationship between TAS2R gene number and special foraging behaviours are also found in marine mammals. Feng et al.  reported that nine of 10 TAS2Rs surveyed in 12 whales were pseudogenized. Consistent with Feng et al’s results but based on a larger genomic scale, we show that all but two of the analysed whales and dolphin lost all the 23 intact TAS2Rs in their most recent common ancestor. Similarly, another group of marine mammals, pinnipedia [walrus (Odobenus rosmarus) and the Weddell seal (Leptonychotes weddellii)] possess only 4–5 intact TAS2Rs (Fig. 2). The bitter taste function in these taxa (cetacea and pinnipedia) might have been relaxed during the evolutionary history of both of them following their colonisation of the marine environment, their dietary switch to consuming fish and cephalopods, and swallowing whole prey [20, 21]. Consistent with this speculation, it has been reported that cetacea have lost most of the sweet, umami and sour tastes receptors and pinnipedia have lost sweet tastes receptors [20, 21].
Purifying selection was suggested to be the major selective force shaping the evolution of taste genes . If selection is relaxed in some species, the corresponding taste genes may become pseudogenes. For example, vampire bats feed solely on a single food item (i.e. blood), and thus have less need for sweet, umami and bitter taste function. As expected, it was recently shown that vampire bats have lost the functional sweet and umami taste genes and pose a significantly greater percentage of pseudogenes of bitter taste genes than other bats [22, 23]. In giant panda the pseudogenization of the umami taste gene coincided with a dietary switch to bamboo . Similarly, penguins are carnivorous, swallow food whole and their tongue structure and function suggest they do not need taste perception, and as expected, they have also lost the functional sweet, umami, and bitter taste genes . These examples suggest that dietary composition and foraging behaviour indeed affect the number of intact taste genes.
TAS2Rs expansions in insectivorous bats and shrew
Bitter tasting compounds (toxins or not) in insects appear to have strongly impacted the diversity of the TAS2Rs repertoire in insectivores. Many insects that fly at night have evolved strategies to avoid predation by bats, of which the chemical defence is an important component: the chemical repertoire of nocturnal insects includes cardiac glycosides and pyrrolizidine alkaloids, which are effective feeding deterrents for bats [8, 9]. Diversity of the TAS2Rs repertoire might thus be of vital importance for insectivorous bats, allowing them to detect toxins and to avoid their ingestion. Similarly, shrews have an insectivorous diet consisting mostly of soil insects [26, 27]. The high metabolic rates of most shrews make them susceptible to food shortage and result in the requirement of a constant food supply . This likely contributes to the critical role of bitter taste for shrews, allowing them to detect toxins and to make prompt decisions while feeding.
The process of the influence of dietary specialization on TAS2R gene number
The influence of dietary specialization on the TAS2Rs number changes could be an interactive, time-dependent, and step by step process. On the one hand, individuals of a species, which started to lose some TAS2R genes caused by dietary specialization may not detect toxins in their ancestral diet, thereby increasing their mortality. Consequently, the loss in TAS2R genes could result in the species acquiring a narrow dietary niche. Conversely, a narrow diet might further reduce the strength of selection for a more diverse TAS2Rs repertoire, accelerating the genomic loss of TAS2Rs. This kind of evolutionary vortex might explain the process of the TAS2Rs loss in marine mammals and Chinese pangolin. On the other hand, individual of one species that acquire new TAS2R genes may be able to detect a wider range of toxins in new potential diets, thereby extending their diet range. The expanded dietary niche might explain the importance of a more diversified TAS2Rs repertoire.
The present study provides the most comprehensive understanding of the TAS2Rs in Laurasiatheria to date, going beyond the known correlation between herbivory and TAS2Rs numbers. Besides herbivores, we show that insectivores tend to have more functional TAS2Rs‚ while lineages with a narrow dietary niche (pangolins) and with less dependence on bitter taste (pangolins, cetacea and pinnipedia) have lost most functional TAS2Rs. Our study reveals that the dietary specialization of Laurasiatheria has had a complex influence on bitter taste gene evolution. However, this finding does not automatically refute the potential impact of other factors than dietary niche on the bitter taste gene evolution. The presence/absence of intact taste genes in mammals and other vertebrates is sometimes inconsistent with what is expected from organismal diet, inferring that our understanding of the physiological functions of these tastes and/or their receptor genes is not complete . Both genetic and behavioural studies, suggested by Behrens et al.  and Jiang et al. , are needed to deepen our understanding of the function and evolution of TAS2Rs.
Data resource and nomenclature
The genome assemblies of 41 Laurasiatherian mammals were downloaded (Additional file 1: Table S1), and the phylogenetic relationships of these species were assessed using TimeTree v3.0 (http://www.timetree.org/, last accessed June 30, 2016) .
Annotation of TAS2Rs
BLASTN and TBLASTN in BLAST v2.2.23  were used for sequence similarity searches in the 41 Laurasiatherian genome assemblies based on templates of nucleoid and protein sequences of the orthologs of TAS2Rs from human, mouse, dog and cow as queries respectively. BLASTN and TBLASTN were performed with filtering off (−F false) and the cutoff E value of 1e-5 (using all the intact and disrupted TAS2Rs as queries) and 1e-10 (using all the intact genes as queries), respectively. The outputs from one Laurasiatherian genome were used for reciprocal BLASTN against this genome itself. Outputs with a length of less than 100 bp were discarded. Sequences longer than 750 bp with complete open reading frame (ORF) and minimal flanking sequences were regarded as candidate complete TAS2Rs. Complete genes were verified by the TMHMM method for the presence of a seven transmembrane domains . The remaining outputs larger or equal to 750 bp were regarded as disrupted TAS2Rs reading frames if they were pseudogenised as indicated by a loss of their start position, a gain of premature stop codons, and/or the presence of indels that disrupted the coding frame. Subsequently, the remaining outputs where the contig end had ≤30 bp in the flanking sequences, which was not disrupted by stop codons or indels according to the multiple alignments with known complete TAS2Rs, were regarded as partial TAS2Rs, because the last 30 bp in the end of contig could be sequenced and assembled incorrectly .
Partial TAS2Rs may be complete genes, but they were not sequenced completely or assembled completely during genome assembly. It has been reported that the probability of overestimating the number of partial TAS2Rs increases when surveying genome assemblies with a contig N50c value of less than 10 kb . However, the N50c values of 40 herein surveyed assemblies are higher than 10 kb (Additional file 1: Table S1). Our data was thus gathered in a way to minimize the probability of such assembly errors, ensuring the reliability of the data. All the overlapping sequences of hits with the same orientations at the contig level were merged (Additional file 1: Table S3). Both complete and partial TAS2Rs were regarded as intact TAS2Rs in subsequent analysis.
Reconstruction of TAS2Rs repertoire evolution
In this study, we follow the nomenclature of TAS2Rs in Euarchontoglires (sister clade of Laurasiatheria) reported by Hayakawa et al. . The annotated TAS2Rs were merged with the TAS2R genes representing all the reported 28 Euarchontoglires TAS2Rs clades  into a multiple alignments of nucleotide sequences using E-INS-i in MAFFT v7 . After removal of gap sites, the phylogenetic trees of TAS2Rs were constructed with MEGA 5 using the neighbour-joining (NJ) method  and with RAxML v8.1.X using the maximum-likelihood (ML) method [35, 36]. The orthologous and paralogous relationships with the TAS2Rs were inferred based on the phylogenetic tree. The evolution of the TAS2R repertoires (births and deaths) in Laurasiatherian lineages was reconstructed as in Hayakawa et al. . Dot-plots were constructed using UGENE v1.10.0  to compare the structure of scaffold/chromosomal regions containing TAS2R loci to detect duplications and deletions.
Phylogenetically independent contrasts (PIC) analyses
The PIC analyses were conducted using the pic function within the R package of Analyses of Phylogenetic and Evolution (ape) in order to identify whether there was any trend in the gain/loss of TAS2Rs along the phylogenetic tree while taking into account the position of species along the phylogenetic tree [12, 38]. The nonparametric Spearman’s rank correlation coefficient (ρ) was used to assess the correlation. The information of topology and branch length of the tree was shown in Fig. 3, which are reconstructed on the basis of TimeTree v3.0  and Zhou et al. .
This project was supported by the National Natural Science Foundation of China (No. 31272301), National Key Program (2016YFC0503200) from Ministry of Science and Technology of China, Foundation of Chinese Academy of Sciences (KSCX3-IOZ-1001 and KSCX2-EW-Q-7-2) and the Youth Innovation Promotion Association of CAS under Grant No. 2011077. We thank Martin Burrows and Wenhua Xiong for assistance in data analysis and writing.
ZL and ML contributed to the design of this research. GL and ZL performed the experiments. GL, ZL, XT, JT and BZ contributed to data analysis. ZL, FH and POtW wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
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
- Dong D, Jones G, Zhang S. Dynamic evolution of bitter taste receptor genes in vertebrates. BMC Evol Biol. 2009;9:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi P, Zhang J. Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol Biol Evol. 2006;23:292–300.View ArticlePubMedGoogle Scholar
- Li D, Zhang J. Diet shapes the evolution of the vertebrate bitter taste receptor gene repertoire. Mol Biol Evol. 2014;31:303–9.View ArticlePubMedGoogle Scholar
- Hu LL, Shi P. Smallest bitter taste receptor (T2Rs) gene repertoire in carnivores. Zool Res. 2013;34:E75–81.PubMedGoogle Scholar
- Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledgebase of divergence times among organisms. Bioinformatics. 2006;22:2971–2.View ArticlePubMedGoogle Scholar
- Wilson DE, Reeder DM. Mammal species of the world, A taxonomic and geographic reference. Baltimore: Johns Hopkins University Press; 2005.Google Scholar
- Hristov N, Conner WE. Effectiveness of tiger moth (Lepidoptera, Arctiidae) chemical defenses against an insectivorous bat (Eptesicus fuscus). Chemoecology. 2005;15:105–13.View ArticleGoogle Scholar
- Ratcliffe JM, Fullard JH. The adaptive function of tiger moth clicks against echolocating bats: an experimental and synthetic approach. J Exp Biol. 2005;208:689–4698.View ArticleGoogle Scholar
- Nishida R. Sequestration of defensive substances from plants by Lepidoptera. Annu Rev Entomol. 2002;47:57–92.View ArticlePubMedGoogle Scholar
- Wang K, Zhao H. Birds generally carry a small repertoire of bitter taste receptor genes. Geno Biol Evol. 2014;7:2705–15.View ArticleGoogle 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:2018–31.View ArticlePubMedGoogle Scholar
- Lei W, Ravoninjohary A, Li X, Margolskee RF, Reed DR, Beauchamp GK, Jiang P. Functional Analyses of Bitter Taste Receptors in Domestic Cats (Felis catus). PLoS One. 2015;10(10):e0139670.View ArticlePubMedPubMed CentralGoogle Scholar
- Paradis E. Analysis of Phylogenetics and Evolution with R. 2nd ed. New York: Springer; 2012.View ArticleGoogle Scholar
- Fisher DO, Owens IP. The comparative method in conservation biology. Trends Ecol Evol. 2004;19:391–8.View ArticlePubMedGoogle Scholar
- Hutchins M, Kleiman DG, Geist V, McDade MC. Grzimek’s Animal Life Encyclopedia, Mammals V. Michigan, vol. Volumes 16. 2nd ed. Farmington Hills: Gale Group; 2003.Google Scholar
- Newman, DMR. The feeds and feeding habits of Old and New World camels. In: Cockrill WR, editor. The camelid: an all-purpose animal. Uppsala: Scandinavian Institute of African Studies; 1984. p. p. 250–92.Google Scholar
- Wackwitz B, Castellaro G, Schwartz HJ, Raggi A. Session: Sustainable Technology Development in Animal Agriculture, Botanical composition of the diets of alpaca (Lama pacos) and llama (Lama glama) in the Andean ragelands of Chile. 1999.Google Scholar
- Swart JM, Richardson PRK, Ferguson JWH. Ecological factors affecting the feeding behaviour of pangolins (Manis temminckii). J Zool. 1999;247:281–92.View ArticleGoogle Scholar
- Feng P, Zheng J, Rossiter SJ, Wang D, Zhao H. Massive Losses of Taste Receptor Genes in Toothed and Baleen Whales. Geno Biol Evol. 2014;6:1254–65.View ArticleGoogle Scholar
- Sato JJ, Wolsan M. Loss or major reduction of umami taste sensation in Pinnipedia. Naturwissenschaften. 2012;99:655–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hong W, Zhao H. Vampire bats exhibit evolutionary reduction of bitter taste receptor genes common to other bats. Proc R Soc B. 2014;281:20141079.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao H, Zhou Y, Pinto CM, Charles-Dominique P, Galindo-González J, Zhang S, Zhang J. Evolution of the sweet taste receptor gene tas1r2 in bats. Mol Biol Evol. 2010;27:2642–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao H, Yang J, Xu H, Zhang J. Pseudogenization of the umami taste receptor gene tas1r1 in the giant panda coincided with its dietary switch to bamboo. Mol Biol Evol. 2010;27:2669–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao H, Li J, Zhang J. Molecular evidence for the loss of three basic tastes in penguins. Curr Biol. 2015;25:R141–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Rychlik L, Jancewicz E. Prey size, prey nutrition and food handling by shrews of different body sizes. Behav Ecol. 2002;13:216–23.View ArticleGoogle Scholar
- Zhou X, Sun F, Xu S, Yang G, Li M. The position of tree shrews in the mammalian tree: Comparing multi-gene analyses with phylogenomic results leaves monophyly of Euarchonta doubtful. Integr Zool. 2015;10:186–98.View ArticlePubMedGoogle Scholar
- Zhao HB, Zhang JZ. Mismatches between feeding ecology and taste receptor evolution: An inconvenient truth. Proc Natl Acad Sci U S A. 2012;103:E1464.View ArticleGoogle Scholar
- Behrens M, Korsching SI, Meyerhof W. Tuning Properties of Avian and Frog Bitter Taste Receptors Dynamically Fit Gene Repertoire sizes. Mol Biol Evol. 2014;31:3216–27.View ArticlePubMedGoogle Scholar
- Jiang P, Josue J, Li X, Glaser D, Li W, Brand JG, Margolskee RF, Reed DR, Beauchamp GK. Major taste loss in carnivorous mammals. Proc Natl Acad Sci U S A. 2012;109:4956–61.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. 1997;25:3389–402.View ArticlePubMedPubMed CentralGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.View ArticlePubMedGoogle Scholar
- Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.PubMedGoogle 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:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.View ArticlePubMedGoogle Scholar
- Okonechnikov K, Golosova O, Fursov M, UGENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28:1166–7.View ArticlePubMedGoogle Scholar
- R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria (http://www.r-project.org); 2015.
- Zhou X, Xu S, Xu J, Chen B, Zhou K, Yang G. Phylogenomic Analysis Resolves the Interordinal Relationships and Rapid Diversification of the Laurasiatherian Mammals. Sys Biol. 2012;61:150–64.View ArticleGoogle Scholar
- Kishida T, Thewissen JGM, Hayakawa T, Imai H, Agata K. Aquatic adaptation and the evolution of smell and taste in whales. Zool Lett. 2015;1:9.View ArticleGoogle Scholar