The Bitter Taste Receptor Genes of the Raccoon Dog (Nyctereutes procyonoides).
Mammals have five basic tastes, including sweet, bitter, umami, sour and salty. Among these tastes, the bitter sense is believed to help animals identify poisonous compounds. Bitter taste receptors coded by bitter taste receptor (Tas2r) genes has been identified in many species, but the evolution and repertoire of Tas2r genes in raccoon dog (Nyctereutes procyonoides) was still unknown. As a true omnivorous mammal, the Tas2r genes may play a critical role in food selection of raccoon dogs. Thus, we explored the evolution of Tas2r genes in raccoon dogs. In this study, Tas2r genes including eleven intact genes and four pseudogenes from raccoon dogs were identified for the first time. We also obtained Tas2r genes from fourteen additional species based on their genome sequences and our previous study, including thirteen carnivorous species and one omnivorous species.
Then, we constructed a phylogenetic tree using Tas2r genes sequences from these fifteen species. Phylogenetic analyses showed that most Tas2r genes of raccoon dogs were closed to other canid-species, indicating that the Tas2r genes were very conservative. Besides, we also detected positive selection in 11 intact genes in raccoon dogs. The results of the positive selection analyses indicated that Tas2r10 and Tas2r67 had positively selected sites and other intact genes in raccoon dogs were under purifying selection. Overall, the phylogenetic and evolutionary relationship of Tas2r genes in raccoon dogs was first studied in this study. Most of the raccoon dog' Tas2r genes were under purifying selection, indicating that its true omnivorous diet was an important role in the Tas2r evolution. Our study provided valuable data of Tas2r genes in raccoon dogs, and provided novel insights into conservation concern of natural populations.
Raccoon dog, Tas2r genes, Phylogenetic analysis, Positive selection, Diet.
Mammals have five basic tastes, including sweet, bitter, umami, sour and salty (Li and Zhang, 2014). Bitter taste receptors which were coded by bitter taste receptor (Tas2r) genes play an important role in identifying the toxins and choosing food in mammals (Chandrashekar et al., 2000). The Tas2r gene repertoires of vertebrates diverge tremendously in size, from 10 pseudogenes in the dolphin (Tursiops truncatus) (Lei et al., 2015) to 80 in the Latimeria chalumnae (Syed and Korsching, 2014). The relationship between the number of Tas2r genes and diets in 54 vertebrate species suggested that dietary toxins were major selective force for shaping the diversity of Tas2r repertoire (Li and Zhang, 2014). By analyzing the evolutionary trajectories of Tas2r genes in the 41 Laurasiatherian species, Liu found that insectivores in Laurasiatherian tended to have more functional Tas2r genes in comparison to carnivoresand omnivores (Liu et al., 2016).
However, no positive correlation were detected between diet composition and Tas2r genes in carnivores (Shang et al., 2017).
Although Tas2r genes have been identified in many species (He et al., 2017), including seven canids, the evolution and repertoire of Tas2r genes in raccoon dog (Nyctereutes procyonoides) were still unknown. As a nocturnal canid, the raccoon dog originally came from East Asia (Ward and Wurster-Hill, 1990). Raccoon dogs have spread over eastern, central and north Europe since its introduction from Asia (Kowalczyk and Bunevich, 2009). Previous studies mainly focused on the diet, habitat use, seasonal physiology, hair and reproduction of the raccoon dog (Kauhala et al., 1993; Hirasawa et al., 2006; Matsuyama et al., 2006; Kowalczyk and Bunevich, 2009). As a true omnivore (Hirasawa et al., 2006), the raccoon dog may face poisonous compounds in food selection, but the research of Tas2r genes in the raccoon dog was obscure. In the present study, we first identified and explored the evolution of Tas2r genes in raccoon dogs.
MATERIALS AND METHODS
Samples, primers and DNA extraction
The sample was obtained from Dalai Lake. Genomic DNA of the raccoon dog was isolated from the muscle tissue using QIAGEN DNeasy kits. The primers were designed by Primer premier 5.0 (Lalitha, 2000), which can be found in Supplementary Material S1. PCR condition was as follow: 95AdegC 30s, 52-58AdegC 40s, 72AdegC 45s, and 72AdegC elongation for 10 min. Then polymerase chain reaction (PCR) productions were sequenced in both directions. The sequences of the newly generated were submitted in GenBank under the accession number (KT426765-KT426779), and additional Tas2r gene sequences used in this study was present in Supplementary Material S2.
This research was approved by the Qufu Normal University Institutional Animal Care and Use Committee with the permit number of QFNU2015-004. This study has no harm for the canid species. Meanwhile, it may help us protect these animals better.
Sequencing data and phylogenetic analysis
The DNA sequences of intact genes were translated into protein sequences by the MEGA 6 (Tamura et al., 2013). Then the protein sequences (Supplementary Material S2) were used for multiple sequence alignment by Multiple Sequence Comparison by Log-Expectation (MUSCLE), free available online at http://www.ebi.ac.uk/Tools/msa/muscle (version 3.8.31) (Edgar, 2004). The software was also used for analyzing the percent identity matrix. The aligned protein data were used for phylogenetic analysis using the software MEGA 6. The neighbor-joining statistical method of analysis and the bootstrap consensus tree was inferred from 1000 replicates (Felsenstein, 1985; Saitou and Nei, 1987). And the zebra fish V1R gene was used as outgroup. The best fitting model was calculated by ProtTest 3.4 (Abascal, 2005).
Selective pressure on Tas2r genes
In order to detect the selective pressure on the intact bitter taste receptor genes in raccoon dogs, we used the site model in the codeml program in PAML (Yang, 2007). The I ratios of nonsynonymous to synonymous substitution (dN/dS) being less than 1, equal to 1 and more than 1 indicate purifying selection, neutral evolution and positive selection, respectively. Twice the log likelihood differences (2Il) were detected in the Chi-square test. The LRT suggest positive selection if there is a significant difference between the two models. Then we use three distinct different models, including random effect likelihood (REL), fixed-effect likelihood (FEL) and single likelihood ancestor counting (SLAC) to analyze the potential selective pressure of 11 intact genes (Kosakovsky-Pond and Frost, 2005).
Tas2r gene amplification and sequencing
The dog Tas2r genes were used to design the primers of the raccoon dogs. The PCR products of Tas2r genes were obtained from the raccoon dogs, including 11 genes (at least 270 codons, start codon, stop codon and seven transmembrane domains) and 4 pseudogenes (Fig. 1). The similarity of the Tas2r protein sequence between dog and raccoon dog ranged from 97.6% to 99.25%, thus the newly Tas2r genes of the raccoon dog were reliable (Supplementary Material S3). Meanwhile, we also used Tas2r genes from 14 additional species (Fig. 1). The results showed that the intact and pseudogenes were different among 15 species. For example, several Tas2r intact genes in raccoon dog were pseudogenized in other species.
In the present study, we carried out a phylogenetic analysis based on the intact Tas2r genes from 15 mammals. Intact Tas2r genes which were homologous to 11 intact Tas2r genes from the raccoon dog were used in the analysis. The pseudogenes and other intact genes were not used in the phylogenetic analysis. Total 145 Tas2r genes and one zebra fish V1r gene were used in the phylogenetic analysis. The result showed that most Tas2r genes of raccoon dog were closed to other canid species (except for Tas2r40 and Tas2r39) (Fig. 2). And Tas2r genes from the Feliodea species tended to cluster together. Tas2r genes of canid species were near to the carnivorous species and far from the omnivorous species.
Table I.- Positive selection of 11 intact Tas2r genes.
Gene###Species###Lnl M7###Lnl M8###P-values###M8###SLAC###FEL###REL
Tas2r5###8###-1459.87###-1457.99###0.15###26 73 76###0###0###0
Tas2r7###9###-1425.02###-1423.17###0.15###15 26 30 31 180 213 308###0###0###0
Tas2r10###8###-1494.8###-1490.98###0.02###6 7 9 40 164 200 318###79
Tas2r38###9###-1486.22###-1483.36###0.057###46 65 105 299 300 312###0###0###0
Tas2r40###9###-1524.8###-1523.44###0.26###84 112 229###0###0###0
Tas2r42###8###-1548.74###-1546.37###0.093###39 63 82 216 279 291 310###0###0###315 321
###315 319 320 321 322
Tas2r67###9###-1573.53###-1565.84###0.00045###86 308 310 312###0###0###86 308 312
Selective pressure of Tas2r genes
To explore the evolutionary forces on the Tas2r genes, we conducted an analysis of site models by using the Tas2r2, Tas2r5, Tas2r7, Tas2r10, Tas2r12, Tas2r38, Tas2r39, Tas2r40, Tas2r41, Tas2r42 and Tas2r67 sequences as input data. The likelihood values, estimated parameters, and predicted positively selected sites were listed in Table I. The LRT was significantly different in Tas2r10 and Tas2r67, indicating that the positive selection was reliable in these genes. The results showed that two and three common positively selected sites were identified in Tas2r10 and Tas2r67, respectively (Table I).
The two common positively selected sites (codon 7and 9) in Tas2r10 were located on the transmembrane region (Fig. 3). And the common positively selected sites in Tas2r67 were located on the transmembrane region (codon 86) and an intracellular domain (codon 308 and 312) (Fig. 3).
The Tas2r genes sequences were identified for the first time in the raccoon dog. The protein sequence identity between dog and raccoon dog ranged from 97.6% to 99.25% (Supplementary Material S3), suggesting that the dog gene-based primer sets can be used for the raccoon dog. Previous researches showed that the carnivorous species obtained 22 or more Tas2r genes (Hu and Shi, 2013; Li and Zhang, 2014; Liu et al., 2016; Shang et al., 2017), while in our research, 11 intact genes and 4 pseudogenes were obtained in the raccoon dog. The incomplete repertoire of Tas2r genes in the raccoon dog may be due to the high similarity but containing gaps in the template sequence at the annealing point of the primer (Monteiro-Ferreira et al., 2015).
Several Tas2r intact genes in raccoon dog were pseudogenized in other species, for example, the Tas2r67 and Tas2r2 were not existed in human, but intact in raccoon dogs and other canid species; the Tas2r39, Tas2r40 and Tas2r41 was pseudogenized in ferret, but was intact in canid species. The pseudogenes and intact genes were different species, indicating that the Tas2r genes were species-specific.
The phylogenetic tree showed that most Tas2r genes of the raccoon dog were clustered together with canid species. We speculated that Tas2r genes were very conservative among species. Thus they could be used as a molecular marker to study the evolutionary relationships among animals (Shi and Zhang, 2006; Hu and Shi, 2013). The phylogenetic results also showed that 11 intact genes from raccoon dogs were intermingled across the tree which indicated that gene duplication events were not existent, which was similar to other carnivorous species (Hu and Shi, 2013). Due to lack of genome sequences in raccoon dog, the incomplete repertoire of Tas2r genes may not reflect the real events. More comparative data are needed to critically test it.
Previous study showed that the repertoire and evolution of Tas2r genes were shaped by many factors, including diet and feeding behaviors (Hong and Zhao, 2014; Li and Zhang, 2014; Liu et al., 2016; Shang et al., 2017; Zhong et al., 2017). Our results showed that the raccoon dog had a small number of Tas2r genes, which might suggest Tas2r genes playing a less important role. But minimal number of functional Tas2r genes did not imply a reduced importance of bitter taste, as it could be compensated by large tuning width (Behrens et al., 2014).
In human, the Tas2r10 was reported to identify 20 natural bitter compounds (Meyerhof et al., 2010). The results of positive selection of Tas2r genes in the raccoon dogs showed that two intact genes had positively selected sites, the Tas2r10 and Tas2r67. Meanwhile, except the Tas2r10 and Tas2r67, the remaining intact genes in raccoon dogs had no positively selected sites, suggesting that most genes were under purifying selection. Raccoon dogs are omnivores that feed on meat (insects, rodents birds, fish, carrion and reptiles) and plants (fruits, nuts and berries) (Sasaki and Kawabata, 1994; Sutor et al., 2010). They faced the potential poisonous compounds, including the bile acids, venom, and skin secretion from carrion, reptiles, insects and toxic plants.
We speculated that Tas2r genes might be functionally important during their evolution due to its true omnivorous diet (Ward and Wurster-Hill, 1990). These taste capabilities could protect raccoon dogs from harmful substances and thus be helpful for their survival.
Overall, the phylogenetic and evolutionary relationship of Tas2r genes in raccoon dogs was first studied in this study. Most of the raccoon dog's Tas2r genes were under purifying selection, indicating that its true omnivorous diet was an important role in the Tas2r evolution. Our study provided valuable data of Tas2r genes in raccoon dogs, and provided novel insights into conservation concern of natural populations.
This work was supported by the Special Fund for Forest Scientific Research in the Public Welfare (No. 201404420) and the National Natural Science Fund of China (No. 31372220, 31672313).
There is supplementary material associated with this article. Access the material online at: http://dx.doi.org/10.17582/journal.pjz/2018.50.4.1361.1366
Statement of conflict of interest
Authors have declared no conflict of interest.
Abascal, F., 2005. Prottest: Selection of best-fit models of protein evolution. Bioinformatics, 21: 2104-2105. https://doi.org/10.1093/bioinformatics/bti263
Behrens, M., Korsching, S. and Meyerhof, W., 2014. Tuning properties of avian and frog bitter taste receptors dynamically fit gene repertoire sizes. Mol. Biol. Evolut., 31: 3216-3227. https://doi.org/10.1093/molbev/msu254
Chandrashekar, J., Mueller, K.L., Hoon, M.A., Adler, E., Feng, L., Guo, W., Zuker, C.S. and Ryba, N.J., 2000. T2rs function as bitter taste receptors. Cell, 100: 703-711. https://doi.org/10.1016/S0092-8674(00)80706-0
Edgar, R.C., 2004. Muscle: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics, 5: 1-19. https://doi.org/10.1186/1471-2105-5-1
Felsenstein, J., 1985. Phylogenies and the comparative method. Am. Natural., 125: 1-15. https://doi.org/10.1086/284325
He, Y., Ma, S., Liu, B., Xue, T., Zhou, Q., Jin, W. and Chen, K., 2017. Molecular cloning and expression of taste receptor gene T2R1 of obscure puffer, takifugu fasciatus. Pakistan J. Zool, 49: 1683- 1691. http://dx.doi.org/10.17582/journal.pjz/2017.49.5.1683.1691
Hirasawa, M., Kanda, E. and Takatsuki, S., 2006. Seasonal food habits of the raccoon dog at a western suburb of Tokyo. Mammal Study, 31: 9-14. https://doi.org/10.3106/1348-6160(2006)31[9:SFHOTR]2.0.CO;2
Hong, W. and Zhao, H., 2014. Vampire bats exhibit evolutionary reduction of bitter taste receptor genes common to other bats. Proc. R. Soc. B: biol. Sci., 281: 20141079. https://doi.org/10.1098/rspb.2014.1079
Hu, L.L. and Shi, P., 2013. Smallest bitter taste receptor (t2rs) gene repertoire in carnivores. Zool. Res., 34: E75.
Kauhala, K., Helle, E. and Taskinen, K., 1993. Home range of the raccoon dog (Nyctereutes procyonoides) in Southern Finland. J. Zool., 231: 95-106. https://doi.org/10.1111/j.1469-7998.1993.tb05355.x
Kosakovsky-Pond, S.L. and Frost, S.D., 2005. Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evolut., 22: 1208. https://doi.org/10.1093/molbev/msi105
Kowalczyk, R. and Bunevich, A.N., 2009. Reproduction and mortality of invasive raccoon dogs (Nyctereutes procyonoides) in the Bialowieza primeval forest (Eastern Poland). Annls. Zool. Fenn., 46: 291-301. https://doi.org/10.5735/086.046.0406
Lalitha, S., 2000. Primer premier 5. Biotech Softw. Internet Rep., 1: 270-272. https://doi.org/10.1089/152791600459894
Lei, W., Ravoninjohary, A., Li, X., Margolskee, R.F., Reed, D.R., Beauchamp, G.K. and Jiang, P., 2015. Functional analyses of bitter taste receptors in domestic cats (Felis catus). PLoS One, 10: e0139670. https://doi.org/10.1371/journal.pone.0139670
Li, D. and Zhang, J., 2014. Diet shapes the evolution of the vertebrate bitter taste receptor gene repertoire. Mol. Biol. Evolut., 31: 303-309. https://doi.org/10.1093/molbev/mst219
Liu, Z., Liu, G., Hailer, F., Orozco-Terwengel, P., Tan, X., Tian, J., Yan, Z., Zhang, B. and Ming, L., 2016. Dietary specialization drives multiple independent losses and gains in the bitter taste gene repertoire of laurasiatherian mammals. Front. Zool., 13: 1-10. https://doi.org/10.1186/s12983-016-0161-1
Matsuyama, J., Hata, K. and Sone, K., 2006. Food habits of raccoon dogs, nyctereutes procyonoides gray, in kagoshima prefecture. Res. Bull. Kagoshima Univ. Forests, 34: 75-80.
Meyerhof, W., Batram, C., Kuhn, C., Brockhoff, A., Chudoba, E., Bufe, B., Appendino, G. and Behrens, M., 2010. The molecular receptive ranges of human tas2r bitter taste receptors. Chem. Senses, 35: 157-170. https://doi.org/10.1093/chemse/bjp092
Monteiro-Ferreira, A., Tomas-Marques, A., Bhide, M., Cubric-Curik, V., Hollung, K., Knight, C.H., Raundrup, K., Lippolis, J., Palmer, M., Sales-Baptista, E., Araujo, S.S. and de Almeida, A.M., 2015. Sequence analysis of bitter taste receptor gene repertoires in different ruminant species. PLoS One, 10: e0124933. https://doi.org/10.1371/journal.pone.0124933
Saitou, N. and Nei, M., 1987. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evolut., 4: 406-425.
Sasaki, H. and Kawabata, M., 1994. Food habits of the raccoon dog Nyctereutes procyonoides viverrinus in a mountainous area of Japan. J. Mammal. Soc. Japan, 19: 1-8.
Shang, S., Wu, X., Chen, J., Zhang, H., Zhong, H., Wei, Q., Yan, J., Li, H., Liu, G. and Sha, W., 2017. The repertoire of bitter taste receptor genes in canids. Amino Acids, 49: 1159-1167. https://doi.org/10.1007/s00726-017-2422-5
Shang, S., Zhang, H., Wu, X., Chen, J., Zhong, H., Wei, Q., Zhao, C., Yan, J., Chen, Y. and Tang, X., 2017. The repertoire of bitter taste receptor genes in ovalentaria fish. Environ. Biol. Fishes, 100: 1489-1492. https://doi.org/10.1007/s10641-017-0659-1
Shi, P. and Zhang, J., 2006. Contrasting modes of evolution between vertebrate sweet/umami receptor genes and bitter receptor genes. Mol. Biol. Evolut., 23: 292-300. https://doi.org/10.1093/molbev/msj028
Sutor, A., Kauhala, K. and Ansorge, H., 2010. Diet of the raccoon dog Nyctereutes procyonoides - A canid with an opportunistic foraging strategy. Acta Theriol., 55: 165-176. https://doi.org/10.4098/j.at.0001-7051.035.2009
Syed, A.S. and Korsching, S.I., 2014. Positive darwinian selection in the singularly large taste receptor gene family of an 'ancient'fish, latimeria chalumnae. BMC Genomics, 15: 1. https://doi.org/10.1186/1471-2164-15-650
Tamura, K., Stecher, G., Peterson, D., Filipski, A. and Kumar, S., 2013. Mega6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evolut., 30: 2725-2729. https://doi.org/10.1093/molbev/mst197
Ward, O.G. and Wurster-Hill, D.H., 1990. Nyctereutes procyonoides. Mammal. Sp., 103: 472-474. https://doi.org/10.2307/3504213
Yang, Z., 2007. Paml 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evolut., 24: 1586-1591. https://doi.org/10.1093/molbev/msm088
Zhong, H., Shang, S., Wu, X., Chen, J., Zhu, W., Yan, J., Li, H. and Zhang, H., 2017. Genomic evidence of bitter taste in snakes and phylogenetic analysis of bitter taste receptor genes in reptiles. Peer J., 5: e3708. https://doi.org/10.7717/peerj.3708
|Printer friendly Cite/link Email Feedback|
|Author:||Shang, Shuai; Zhang, Honghai; Wu, Xiaoyang; Wei, Qinguo; Chen, Jun; Zhang, Huanxin; Zhong, Huaming;|
|Publication:||Pakistan Journal of Zoology|
|Date:||Aug 31, 2018|
|Previous Article:||Therapeutic Effects of Curcuma longa against Rotenone-Induced Gross Motor Skill Deficits in Rats.|
|Next Article:||Gene Expression in Escherichia coli and Purification of Recombinant Type II Pullulanase from a Hyperthermophilic Archaeon, Pyrobaculum calidifontis.|