Why Does the Giant Panda Eat Bamboo? A Comparative Analysis of Appetite-Reward-Related Genes among Mammals

Abstract

Background

The giant panda has an interesting bamboo diet unlike the other species in the order of Carnivora. The umami taste receptor gene T1R1 has been identified as a pseudogene during its genome sequencing project and confirmed using a different giant panda sample. The estimated mutation time for this gene is about 4.2 Myr. Such mutation coincided with the giant panda's dietary change and also reinforced its herbivorous life style. However, as this gene is preserved in herbivores such as cow and horse, we need to look for other reasons behind the giant panda's diet switch.

Methodology/Principal Findings

Since taste is part of the reward properties of food related to its energy and nutrition contents, we did a systematic analysis on those genes involved in the appetite-reward system for the giant panda. We extracted the giant panda sequence information for those genes and compared with the human sequence first and then with seven other species including chimpanzee, mouse, rat, dog, cat, horse, and cow. Orthologs in panda were further analyzed based on the coding region, Kozak consensus sequence, and potential microRNA binding of those genes.

Conclusions/Significance

Our results revealed an interesting dopamine metabolic involvement in the panda's food choice. This finding suggests a new direction for molecular evolution studies behind the panda's dietary switch.

About the Authors:

Ke Jin

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Chenyi Xue

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Xiaoli Wu

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Jinyi Qian

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Yong Zhu

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Zhen Yang

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Takahiro Yonezawa

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

M. James C. Crabbe

Affiliation: Faculty of Creative Arts, Technologies and Science, Institute of Applied Natural Sciences, University of Bedfordshire, Luton, United Kingdom

Ying Cao

Affiliation: Institute of Statistical Mathematics, Tokyo, Japan

Masami Hasegawa

Affiliations School of Life Sciences, Fudan University, Shanghai, China, Institute of Statistical Mathematics, Tokyo, Japan

Yang Zhong

* E-mail: [email protected] (Y. Zheng); [email protected] (Y. Zhong)

Affiliations School of Life Sciences, Fudan University, Shanghai, China, Institute of Biodiversity Science and Institute of High Altitude Medicine, Tibet University, Lhasa, China

Yufang Zheng

* E-mail: [email protected] (Y. Zheng); [email protected] (Y. Zhong)

Affiliation: School of Life Sciences, Fudan University, Shanghai, China

Introduction

Biodiversity is showing up not only how different animals appear, but also on how different their diets are. Some species have unique diets and it is interesting to look for the evolutionary reasons behind such natural selection. One such example is the giant panda. As one species of Ursidae (the bear family), Carnivora, the giant panda (Ailuropoda melanoleuca) has a very special bamboo diet. The bears are normally carnivores to omnivores. For example, the polar bears survive solely on meat and fat, while the black bears are more or less omnivores. About 7 Myr ago, the ancient giant panda was still omnivorous [1]. About 2–2.4 Myr ago, they become herbivores as soft bamboo shoots, stems, and leaves became their major food source [1]. The giant panda also developed an enlarged radial sesamoid that functions as a thumb to grab bamboo [2], [3], [4]. The bamboos have relatively low energy and nutrition compare to meat and fruits. Therefore, the average giant panda needs to consume up to 6% of its body weight of bamboo in dry matter per day to survive [5], [6], [7]. However, the giant panda's digestive system is still more fit for a meat diet than bamboo as they can use less than 20% of the bamboo they ingest [7], [8], and both its gut anatomy and microbiome have not yet adapted to degrade those bamboo fibers [9]. Therefore, what is the driving force behind the panda's diet switch is still a question for evolutionary biologists.

Recently, two papers reported that the T1R1 gene in the giant panda turned into a pseudogene due to two frame-shifting mutations in exon 3 and 6, respectively [8], [10]. T1R1 is part of the T1R1/T1R3 heterodimer receptor that mediates umami taste. Thus, this lost-of-function on the T1R1 gene in the giant panda may contribute to the panda's food choice [8], [10]. This is a breakthrough on how the giant panda has become the species it is. However, there are some missing links between the lost meat taste and the bamboo diet of the giant panda. First of all, taste is not the only environmental cue that can affect animals' eating behaviours. To survive, energy and nutrition properties of food can highly influence animals' food choices, while smell and taste are associated with those properties [11], [12]. Therefore, even without the ability to taste meat, giant pandas can still choose meat as their main diet since meat contains much higher energy and nutrition than bamboo. Meat is also available in the giant panda's habitats as other carnivores such as wolf and dog share the same area. Secondly, the estimated mutation time for the T1R1 gene is probably 4.2 Myr for the giant panda [10]. The fossil evidence showed that the giant panda started eating bamboo at least 7 Myr ago, and at about 2.0–2.4 Myr ago they probably had already completed their dietary switch [1]. Therefore, the pseudogenization of T1R1 is probably the result of, not the reason for, its dietary change. Thirdly, the T1R1 gene is intact in some herbivores such as cow and horse [10], which indicates that the taste is probably not the only reason for an animal's food choices. Therefore, we need to look beyond taste to understand the driving force for the panda's diet.

To answer such a question, we also need to look on how and why taste can trigger eating behaviour. In the wild, animals are normally attracted by sweet and umami (savory) stimuli. The umami taste is excited by L-glutamate, which is abundant in meat [13], [14]. The sweet taste is excited by saccharides [13], [14]. Both chemicals indicate the food is high in nutrition and energy, which is crucial for an animal's survival in the wild. Therefore, such a taste is related with the physical rewarding properties of food that in turn further affect the feeding behaviour of the animal. In other words, both the cues (such as umami taste) and properties (such as rich in nutrition and energy) of the food are stored in memory to guide future behaviour, such as to orient the animal back to the source of food [11], [12]. To the giant panda, without being able to taste the meat, the nutrition and energy properties of food should play a more critical role in driving its feeding behaviour. However, the bamboo diet does not fit to this hypothesis, suggesting that there might be something special in the appetite-reward system for the giant panda. To understand such bias, we looked deep in the appetite-reward circuitry and did a comprehensive analysis on genes involved in this appetite-reward system based on the giant panda's draft genome published last year [8]. Our results revealed a complex genetic background and an interesting dopamine metabolic involvement behind the taste for the giant panda's bamboo diet.

Results and Discussion

The tendency to engage in or maintain feeding behaviour is potently influenced by the flavour of food, the gut reaction to the components of food, and the reward pathways in the brain [11], [15], [16], [17], [18]. Animal experiments have shown that both opioid and dopamine are related to appetite-reward circuitry in food intake behaviour [17], [18]. Based on those studies in human and rodents [11], [15], [16], [17], [18], we undertook a comprehensive analysis on about 166 major genes (see a list of gene full name and Ensembl number in table S1) in the panda that are involved in the appetite and food reward system for their coding regions and untranslated regions, especially the 3′UTR and ATG region. We also constructed an online database for analysis on those genes, which can be accessed at http://idm.fudan.edu.cn/Apanda/ (Fig S3).

1. Analysis on coding region

We first screened out the coding region of those selected genes. To look for possible unique sites for the giant panda, we first compared the sequences of giant panda with human (Homo sapiens). Interestingly, most of those genes are highly conserved which is probably due to the biological importance of food intake for survival. Therefore, to look for the possible structural changes of those proteins based on the panda's sequence, we then screened those genes based on their similarity to human orthologs. The screening was performed following two criteria: (i)>75% consensus sequence between panda and human; and (ii) with human protein structures reported in PDB (sequence coverage >75%). Six genes of the panda were screened out (Table 1, COMT, MAOA, MAOB, LEP, ALDH2, PNMT) and structural simulations for those six proteins were done based on the structure of human orthologs by Modeller 9v8 software [19]. Our results revealed that most of them had no major changes in the panda (Figure S1A for MAOA), except one, COMT (catechol-o-methyltransferase), which had a significant change in its substrate/cofactor SAM (S-adenosyl-L-methionine) binding domain (Figure 1A). COMT is one of enzymes that inactive catecholamine neurotransmitters, such as dopamine, epinephrine, and norepinephrine [20], [21]. The enzyme introduces a methyl group to the catecholamine, which is donated by the CE methyl group on the methionine portion of SAM [22], [23]. The methionine portion of SAM is fixed to proper position by hydrogen bonds with several crucial residues in α3 (V42), α4 (S72), and β4 (D141) in human COMT [23], [24]. Although those crucial amino acids were conserved in panda, the α4 helix is lost due to a four amino acids deletion mutation in panda COMT and this part became a more flexible loop in the panda's COMT structure (Figure 1A, blue square & Figure 2 for alignment). This flexible loop may very likely affect the conformation of the methionine portion of SAM and the methyl transfer between the CE methyl group of SAM and catecholamines. Therefore, such structural variation of the panda's COMT suggests a possible relatively slow metabolic turnover rate of catecholamine neurotransmitters by COMT in panda. It has been shown that COMT-deficient male mice had elevated dopamine levels in the frontal cortex [25], and elevated dopamine could enhance motivation for food in mice [16]. In humans, several mutations on COMT were associated with eating disorders [26], [27] and obsessive-compulsive disorder in men [28]. Therefore, it is very likely that the potential lower activity of COMT in panda is related to its special feeding behavior.

[Figure omitted. See PDF.]

Figure 1. Panda and dog COMT protein structural simulation.

The simulated panda or dog proteins were indicated with pink color for backbone and red color for special amino acid residues. Human proteins were indicated with blue color for backbone and yellow color for special amino acid residues. The side chain of COMT catecholamine substrate binding sites (Lys144, Asn170, Glu199) and SAM binding sites (Val42, Ser72, Glu90, Asp141) were shown in ball and stick model. (A) The simulated panda COMT structure compared with human COMT. The α4 helix in blue square has turned into a loop in panda. (B) The simulated dog COMT structure compared with human COMT.

https://doi.org/10.1371/journal.pone.0022602.g001

[Figure omitted. See PDF.]

Figure 2. Sequence alignment of COMT from all nine species.

The number for amino acids is based on human soluble COMT. The conserved amino acids were highlighted. The regions for alpha helix and beta sheet were marked at the bottom of the alignment (wave for alpha helix and arrow for beta sheet). The region with missing amino acids in panda is highlighted with purple square box.

https://doi.org/10.1371/journal.pone.0022602.g002

[Figure omitted. See PDF.]

Table 1. Six genes were selected out based on protein structure and sequence consensus percentages.

https://doi.org/10.1371/journal.pone.0022602.t001

To further check the evolutional effect on those six genes, we then checked the orthologs of seven other species from Ensembl, including chimpanzee (Pan troglodytes), mouse (Mus musculus), rat (Rattus norvegicus), dog (Canis lupus familiaris), cat (Felis catus), horse (Equus caballus), and cow (Bos taurus). Those species have various diets as indicated in Li's paper [8]. For example, cat and dog are carnivores, horse and cow are herbivores, and the others are omnivores. We performed the structure simulation for the six proteins from those species and compared to human orthologs. In general, the structures are highly conserved and very similar to human orthologs. As for COMT, the substrate binding domains of other species were very similar to human COMT even though there is some small shift on the whole structure (Figure 1B for dog COMT and others in our database). The sequence for the α4 helix is totally conserved in all species we checked except panda (Figure 2).

2. Analysis of Kozak motif at the ATG starting site

Besides the coding region, a gene's expression level is highly regulated by multiple ways, including transcription regulation and post-transcription regulation. The 5′UTR region is critical in transcription regulation as there are several mechanisms involved, such as transcription factors binding to the promoter region to modify the efficiency of transcription, the start of a protein-coding sequence that affects the polymerase complex binding to DNA, etc. Since the transcription factor binding sites are complex, as there are both enhancing and depressing effects by transcription factors, it is inconclusive to analyze them based solely on DNA sequence information. Therefore, we decided to focus on the starting site of those genes, which is the Kozak consensus sequence.

The Kozak consensus sequence is a sequence that occurs on eukaryotic mRNA and has the consensus gccRccAUGG in vertebrates [29], [30], [31]. Some nucleotides in this sequence are more important than others: the AUG is essential since it is the actual initiation codon. For a ‘strong’ consensus, the nucleotides at positions +4 and -3 must both match the consensus. An ‘adequate’ consensus has only 1 of these sites, while a ‘weak’ consensus has neither. Although a G in the -6 position and CC at -1 and -2 contribute to the overall strength, they are less essential in initiating translation [29], [30]. Therefore, to compare the strength between human and panda Kozak sequence, we used AUG, -3, and +4 positions as criteria (Table 2 and Table S2). To eliminate the effects of sequencing errors, we first removed the genes without the AUG initiation codon from the 166 genes targeted. There are 78 genes removed and the 88 genes left were divided into two parts, (i) 7 genes present the different pattern at -3 or +4 position or both positions (Table 2, MC4R, OPRD1, COMT, ADRA1D, GRIA3, HTR3E, GRM7), and (ii) 81 genes have the same pattern at -3 and +4 position (Table S2). Four of those seven genes (COMT, MC4R, OPRD1, GRM7) have weaker Kozak consensus sequence in the giant panda that probably can also cause a lower expression of the respective gene (Table 2, gene name in bold). Interestingly, COMT is also included in this list and such a weaker Kozak motif probably can further enhance the effect of its coding region.

[Figure omitted. See PDF.]

Table 2. Seven genes with different Kozak sequence between panda and human.

https://doi.org/10.1371/journal.pone.0022602.t002

We also checked the Kozak motif of those seven genes in other species and all species have its own expression pattern for those seven genes. Such variations and diversity may somehow relate to the complex genomic background and specific characters of each species (Table S3).

3. Analysis on potential MicroRNA binding on 3′UTR

Post-transcription regulation has also been recognized more and more importantly in regulating gene expression level. After being produced, the stability and distribution of the different transcripts can be regulated by several important mechanisms such as RNA interference (RNAi) and RNA binding protein. Two types of small RNA molecules, microRNA (miRNA) and small interfering RNA (siRNA) are central to RNAi. MiRNAs are post-transcriptional regulators that bind to complementary sequences on target messenger RNA transcripts (mRNAs), usually resulting in translational repression and gene silencing [32]. MiRNAs are well conserved in eukaryotic organisms and are thought to be a vital and evolutionary ancient component of genetic regulation [33], [34]. Therefore, we focused on those genes that have been screened out by the first two methods to predict potential miRNA binding. There are 6 genes screened out by the first method and 7 genes screened out by the second method, and COMT is screened out by both methods. Therefore, there are total 12 genes checked to predict potential miRNA binding. 9 out of 12 genes have at least one computational predicted miRNA (Table S4). Two of those genes, COMT and ADRA1D (adrenergic receptor, alpha-1D), were predicted to be regulated by miRNA in panda but not in human. MiR-30C was computationally predicted to bind the 3′UTR of panda ADRA1D gene and miR-199a-5p to the 3′UTR of panda COMT gene (Figure 3 & Figure S2). The adrenergic receptors are G protein-coupled receptors that are targets of the catecholamine, especially norepinephrine and epinephrine. The predicted miRNA regulation may down-regulate the COMT and ADRA1D's protein expression level in the giant panda and enhance the effects of deficient catecholamine metabolism.

[Figure omitted. See PDF.]

Figure 3. Part of the panda COMT 3′UTR and predicted miRNA-199a-5p binding.

COMT 3′UTR secondary structure was predicted by Mfold and part of it was shown on the left. The possible miRNA-199a-5p target with panda COMT gene and the calculated free energy was shown on the right side.

https://doi.org/10.1371/journal.pone.0022602.g003

Conclusion

The giant panda's dietary switch to bamboo is unique and interesting. The pseudogenization of its umami receptor gene T1R1 is coincident with such bamboo dietary switch. But there are probably other factors involved in such a switch. To look beyond taste and the microbiome in panda that has the potential affect on this dietary switch, we studied the giant panda's appetite-reward system.

Our comprehensive sequence analysis on the giant panda's appetite-reward systems indicated that the panda probably has some defects in its catecholamine metabolic pathways, which thus affect its food choices. This prediction provides a new insight for better understanding the giant panda's specific characteristics. However, eating behaviour is an extremely complex process. Besides the homeostatic regulation by the hypothalamus, the cortical-limbic system is also involved in regulating response to environmental conditions and stimuli such as the smell and taste of food [1], [3], [8]. In addition, too much or too little dopamine both has profound effects on feeding behaviours [2]. It would be interesting to carry out detailed biochemical assay on those enzymes, such as COMT, that are involved in appetite-reward circuitry and compare them to other family of bears and other species, e.g. dog.

Based on our analysis, the dopamine metabolic system is probably not that competent in the giant panda. Recently, it has been shown that dopamine is essential for stimulus-reward learning behaviour [35]. Therefore, such deficient dopamine metabolism in the giant panda may have some profound effect on its reward-cue directed behaviour. It is possible that some ingredient in the bamboo may able to help the catecholamine, especially dopamine, metabolism in panda. Such an ingredient may able to stimulate the appetite-reward circuitry in the giant panda and play certain role in setting up its bamboo diet. Further analysis on bamboo chemical components and their effects on nerve system will be needed.

Even then, the initial evolutionary reason for the giant panda to choose such a diet is still a mystery. Further investigation on the carnivore catecholamine metabolism and its effect on food choice will be essential. In addition, comparative genomics of the bear family, Ursidae, as well as Ailurus fulgens (red panda), which has a similar diet to the giant panda, will be helpful to understand carnivore food choices.

Materials and Methods

The genome sequence of the Giant Panda was downloaded from BGI (Beijing Genomics Institute, http://panda.genomics.org.cn/page/panda/index.jsp) and compared with the data from Ensembl (Ensembl Genome Browser). The high quality genome and reference sequence of eight other species were also retrieved from Ensembl, which include human (Homo sapiens), chimpanzee (Pan troglodytes), mouse (Mus musculus), rat (Rattus norvegicus), dog (Canis lupus familiaris), cat (Felis catus), horse (Equus caballus), and cow (Bos taurus). The 166 targeted genes are selected based on the studies about appetite and food intake behaviour [11], [15], [16], [17], [18]. Orthologs were determined and mapped to Ensembl gene identifiers using BioMart [36], [37]. Please see a list of gene name, full name, and Ensembl number in our database and supplementary table S1.

Panda protein 3D-structure simulation

To see how the different amino acids in a panda protein would affect its structure, we carried out the following simulation analysis. First, we filtered the datasets containing 166 proteins by two criteria (i) the sequence consensus between the protein in panda and its ortholog in human is above 75%; and (ii) the gene's ortholog in human must have 3D-structure with sequence coverage above 75%. With these two criteria, we filtered out six proteins (Table 1) and simulated panda proteins using Modeller 9v8 [19] to compare with their human orthologs. Then we also performed structure simulation of those six proteins for the other seven species.

For the COMT protein, we performed structure simulation based on the 3D-structure of human COMT protein (PDB ID: 3BWM). We obtained the simulated 3D-structure of COMT protein in panda and aligned with human COMT to visualize the difference (Figure 1A). The dog and human COMT structure alignment was presented in Figure 1B.

By using the strategy described above, we performed structure simulation in panda on the 3D-structure of human MAOA protein (PDB ID: 2BXR) and aligned with its human ortholog. These two structures were presented in Figure S1A. The simulated dog MAOA aligned with human ortholog were presented in Figure S1B.The simulated structures for other species are stored in our database.

For the other four proteins (LEP, ALDH2, PNMT, MAOB), the sequence between human and panda are highly conserved (Table 1). Similar methods were applied and all simulation structures are stored in our database (http://idm.fudan.edu.cn/Apanda/) (Figure S3).

Sequence alignment of COMT

The alignment of COMT from nine species were performed by MAFFT 6.850 [38], which has an auto mode allowing the program to choose proper parameters for each sequence sets.

Kozak motif analysis

For all previous described 166 genes, we carried out Kozak motif analysis to evaluate the Kozak consensus of each gene within nine species. For all genes, we downloaded the DNA genome sequences from Ensembl. Based on the annotation of the released genome sequences on Ensembl, we identified where the coding sequence starts. To eliminate the effects of sequencing errors, we first removed 78 panda genes without the AUG initiation codon from the 166 genes targeted. Then we extracted the 30 bp nucleotide sequences of each gene at position −25∼+5 based on the annotation for those 88 panda genes left. If the sequence included the nucleotides without masked characters, we extracted the Kozak sequence directly from them. Otherwise, we manually identified a 30 bp sequence between position −25∼+5 by using the NCBI sequence viewer.

Prediction of 3′ UTR in Panda

The genes screened out by the first two methods were combined here to check for potential miRNA binding. Since COMT was screened out by both methods, there are 12 genes (6+7−1 = 12) checked here for potential miRNA binding. As the 5′- and 3′- untranslated regions (UTR) for genes in the panda genome are without annotation, 4 kb sequences downstream of the coding sequence (CDS) of the above 12 genes were extracted for 3′UTR analysis. We used the Polyadq program [39] with default settings to predict the polyadenylation site (polyA) (http://rulai.cshl.org/tools/polyadq/polyadq_form.html), which is a major signal for transcription termination. This sequence was also compared with the multiple sequence alignment file of the 3′UTR of each human gene obtained from the TargetScanS web server [40], [41], [42] (http://www.targetscan.org/).

For COMT, 604 bps downstream of the CDS (details can be accessed at http://idm.fudan.edu.cn/Apanda/) was determined as the 3′UTR of the COMT gene. The secondary structure of this sequence was predicted by Mfold with default settings [43] and the best predicted structure was selected (Figure S2).

For MAOA, the length of 3′UTR of human MAOA gene is 2,307 bps in the annotation of TargetScanS and the coding region of MAOA gene is on the negative strand. Therefore, a 4 kb sequence on positive strand upstream of the MAOA gene was selected, and the reverse complimentary sequences were analyzed by the same method described above to predict the polyA site of MAOA gene in panda. This sequence was also compared with the multiple sequences alignment file of the 3′UTR of human MAOA gene obtained from the TargetScanS. Finally, a 2320 bps sequence (details can be accessed at http://idm.fudan.edu.cn/Apanda/) was determined as the 3′UTR of the MAOA gene in panda.

For the other 10 genes (LEP, ALDH2, PNMT, MAOB, ADRA1D, MC4R, OPRD1, GRIA3, HTR3E, GRM7), detailed information can be seen in our database (http://idm.fudan.edu.cn/Apanda/).

MicroRNA Target Prediction in Panda

For the above 12 genes, the same protocol was used to predict potential panda miRNA target. In detail, we first obtained all the miRNAs that have at least one target site in the 3′UTR of each human gene by TargetScanS. Then, those highly conserved miRNAs were filtered out by the following criteria: (i) miRNA should be present in at least six species; and (ii) only allow one mismatch in 22 nucleotides in all species. MiRNA sequences were downloaded from miRBase [44] (February, 2010). For those 12 genes, potential target sites on the 3′UTR of the filtered panda miRNAs were predicted by RNA22 server [45]. To ensure accuracy and sensitivity, we changed two parameter settings: (i) seed/nucleus is changed from 6 (default setting) to 7; and (ii) maximum folding energy for heteroduplex (Kcal/mol) is changed from −25 (default setting) to −20. The positive results for 12 genes were marked with “+” in table S3.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

The simulated panda and dog MAOA protein structure compared with human MAOA. The panda or dog proteins were indicated with pink color for backbone and red color for special amino acid residues. Human proteins were indicated with blue color for backbone and yellow color for special amino acid residues. (A) The simulated panda MAOA compared with human MAOA. (B) The simulated dog MAOA compared with human MAOA.

https://doi.org/10.1371/journal.pone.0022602.s001

(TIF)

Figure S2.

The secondary structure of 3′- UTR of gene COMT predicted by Mfold. The part in red box was enlarged and put into Figure 3.

https://doi.org/10.1371/journal.pone.0022602.s002

(TIF)

Figure S3.

A screen-shot of our giant panda genome analysis database. (http://idm.fudan.edu.cn/Apanda/).

https://doi.org/10.1371/journal.pone.0022602.s003

(TIF)

Table S1.

List of the 166 genes' full name and Ensembl number. Empty box indicated that there is no ortholog found in that species.

https://doi.org/10.1371/journal.pone.0022602.s004

(XLS)

Table S2.

Kozak sequence pattern of 81 genes in panda and human. The genes have the same pattern at −3 and +4 position. A “y” is marked to indicate them matching the consensus sequence.

https://doi.org/10.1371/journal.pone.0022602.s005

(DOC)

Table S3.

Kozak sequence pattern of 7 genes in all nine species. The genes have different pattern at -3 and/or +4 position between human and panda. A “Y” is marked to indicate them matching the consensus sequence. An “X” is marked to indicate no ortholog identified in that species. A “-” is marked to indicate no “AUG” start codon in that species.

https://doi.org/10.1371/journal.pone.0022602.s006

(DOC)

Table S4.

Predicted miRNAs for those 12 genes screened out by first two methods. MiRNAs for those genes were predicted by RNA22 and there were no predicted miRNA for PNMT, GRIA3 and MC4R.

https://doi.org/10.1371/journal.pone.0022602.s007

(DOC)

Acknowledgments

We would like to thank the Associate Editor and two referees for their constructive comments which have significantly improved the quality of this article.

Author Contributions

Conceived and designed the experiments: Y. Zhong Y. Zheng. Performed the experiments: KJ CX XW JQ Y. Zhu ZY. Analyzed the data: KJ CX XW JQ Y. Zhu ZY. Wrote the paper: KJ TY YC MJCC MH Y. Zhong Y. Zheng.

Citation: Jin K, Xue C, Wu X, Qian J, Zhu Y, Yang Z, et al. (2011) Why Does the Giant Panda Eat Bamboo? A Comparative Analysis of Appetite-Reward-Related Genes among Mammals. PLoS ONE6(7): e22602. https://doi.org/10.1371/journal.pone.0022602

References

1. Jin C, Ciochon RL, Dong W, Hunt RM Jr, Liu J, et al. (2007) The first skull of the earliest giant panda. Proc Natl Acad Sci U S A 104: 10932–10937.C. JinRL CiochonW. DongRM Hunt JrJ. Liu2007The first skull of the earliest giant panda.Proc Natl Acad Sci U S A1041093210937

2. Endo H, Sasaki N, Yamagiwa D, Uetake Y, Kurohmaru M, et al. (1996) Functional anatomy of the radial sesamoid bone in the giant panda (Ailuropoda melanoleuca). J Anat 189 ( Pt 3): 587–592.H. EndoN. SasakiD. YamagiwaY. UetakeM. Kurohmaru1996Functional anatomy of the radial sesamoid bone in the giant panda (Ailuropoda melanoleuca).J Anat189 ( Pt 3)587592

3. Endo H, Yamagiwa D, Hayashi Y, Koie H, Yamaya Y, et al. (1999) Role of the giant panda's ‘pseudo-thumb’. Nature 397: 309–310.H. EndoD. YamagiwaY. HayashiH. KoieY. Yamaya1999Role of the giant panda's ‘pseudo-thumb’.Nature397309310

4. Salesa MJ, Anton M, Peigne S, Morales J (2006) Evidence of a false thumb in a fossil carnivore clarifies the evolution of pandas. Proc Natl Acad Sci U S A 103: 379–382.MJ SalesaM. AntonS. PeigneJ. Morales2006Evidence of a false thumb in a fossil carnivore clarifies the evolution of pandas.Proc Natl Acad Sci U S A103379382

5. Schaller GB (1985) The Giant pandas of Wolong. Chicago: University of Chicago Press. xix, 298.GB Schaller1985The Giant pandas of Wolong.Chicago: University of Chicago Press. xix,298

6. Gittleman JL (1989) Carnivore behavior, ecology, and evolution. Ithaca: Comstock Pub. Associates. JL Gittleman1989Carnivore behavior, ecology, and evolution.IthacaComstock Pub. Associates

7. Dierenfeld ES, Hintz HF, Robertson JB, Van Soest PJ, Oftedal OT (1982) Utilization of bamboo by the giant panda. J Nutr 112: 636–641.ES DierenfeldHF HintzJB RobertsonPJ Van SoestOT Oftedal1982Utilization of bamboo by the giant panda.J Nutr112636641

8. Li R, Fan W, Tian G, Zhu H, He L, et al. (2010) The sequence and de novo assembly of the giant panda genome. Nature 463: 311–317.R. LiW. FanG. TianH. ZhuL. He2010The sequence and de novo assembly of the giant panda genome.Nature463311317

9. Tuohy KM, Hinton DJ, Davies SJ, Crabbe MJ, Gibson GR, et al. (2006) Metabolism of Maillard reaction products by the human gut microbiota--implications for health. Mol Nutr Food Res 50: 847–857.KM TuohyDJ HintonSJ DaviesMJ CrabbeGR Gibson2006Metabolism of Maillard reaction products by the human gut microbiota--implications for health.Mol Nutr Food Res50847857

10. Zhao H, Yang JR, Xu H, Zhang J (2010) Pseudogenization of the umami taste receptor gene Tas1r1 in the giant panda coincided with its dietary switch to bamboo. Mol Biol Evol 27: 2669–2673.H. ZhaoJR YangH. XuJ. Zhang2010Pseudogenization of the umami taste receptor gene Tas1r1 in the giant panda coincided with its dietary switch to bamboo.Mol Biol Evol2726692673

11. Fulton S (2010) Appetite and reward. Front Neuroendocrinol 31: 85–103.S. Fulton2010Appetite and reward.Front Neuroendocrinol3185103

12. Berthoud HR, Morrison C (2008) The brain, appetite, and obesity. Annu Rev Psychol 59: 55–92.HR BerthoudC. Morrison2008The brain, appetite, and obesity.Annu Rev Psychol595592

13. Zhao GQ, Zhang Y, Hoon MA, Chandrashekar J, Erlenbach I, et al. (2003) The receptors for mammalian sweet and umami taste. Cell 115: 255–266.GQ ZhaoY. ZhangMA HoonJ. ChandrashekarI. Erlenbach2003The receptors for mammalian sweet and umami taste.Cell115255266

14. Chandrashekar J, Hoon MA, Ryba NJ, Zuker CS (2006) The receptors and cells for mammalian taste. Nature 444: 288–294.J. ChandrashekarMA HoonNJ RybaCS Zuker2006The receptors and cells for mammalian taste.Nature444288294

15. Palmiter RD (2007) Is dopamine a physiologically relevant mediator of feeding behavior? Trends Neurosci 30: 375–381.RD Palmiter2007Is dopamine a physiologically relevant mediator of feeding behavior?Trends Neurosci30375381

16. Wise RA (2004) Dopamine, learning and motivation. Nat Rev Neurosci 5: 483–494.RA Wise2004Dopamine, learning and motivation.Nat Rev Neurosci5483494

17. Lenard NR, Berthoud HR (2008) Central and peripheral regulation of food intake and physical activity: pathways and genes. Obesity (Silver Spring) 16: Suppl 3S11–22.NR LenardHR Berthoud2008Central and peripheral regulation of food intake and physical activity: pathways and genes.Obesity (Silver Spring)16Suppl 3S1122

18. Rask-Andersen M, Olszewski PK, Levine AS, Schioth HB (2009) Molecular mechanisms underlying anorexia nervosa: focus on human gene association studies and systems controlling food intake. Brain Res Rev 62: 147–164.M. Rask-AndersenPK OlszewskiAS LevineHB Schioth2009Molecular mechanisms underlying anorexia nervosa: focus on human gene association studies and systems controlling food intake.Brain Res Rev62147164

19. Eswar N, Webb B, Marti-Renom MA, Madhusudhan MS, Eramian D, et al. (2006) Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics Chapter 5: Unit 5 6. N. EswarB. WebbMA Marti-RenomMS MadhusudhanD. Eramian2006Comparative protein structure modeling using Modeller.Curr Protoc Bioinformatics Chapter 5: Unit 5 6

20. Haavik J, Blau N, Thony B (2008) Mutations in human monoamine-related neurotransmitter pathway genes. Hum Mutat 29: 891–902.J. HaavikN. BlauB. Thony2008Mutations in human monoamine-related neurotransmitter pathway genes.Hum Mutat29891902

21. Roth JA (1992) Membrane-bound catechol-O-methyltransferase: a reevaluation of its role in the O-methylation of the catecholamine neurotransmitters. Rev Physiol Biochem Pharmacol 120: 1–29.JA Roth1992Membrane-bound catechol-O-methyltransferase: a reevaluation of its role in the O-methylation of the catecholamine neurotransmitters.Rev Physiol Biochem Pharmacol120129

22. Mannisto PT, Kaakkola S (1999) Catechol-O-methyltransferase (COMT): biochemistry, molecular biology, pharmacology, and clinical efficacy of the new selective COMT inhibitors. Pharmacol Rev 51: 593–628.PT MannistoS. Kaakkola1999Catechol-O-methyltransferase (COMT): biochemistry, molecular biology, pharmacology, and clinical efficacy of the new selective COMT inhibitors.Pharmacol Rev51593628

23. Rutherford K, Le Trong I, Stenkamp RE, Parson WW (2008) Crystal structures of human 108V and 108M catechol O-methyltransferase. J Mol Biol 380: 120–130.K. RutherfordI. Le TrongRE StenkampWW Parson2008Crystal structures of human 108V and 108M catechol O-methyltransferase.J Mol Biol380120130

24. Vidgren J, Svensson LA, Liljas A (1994) Crystal structure of catechol O-methyltransferase. Nature 368: 354–358.J. VidgrenLA SvenssonA. Liljas1994Crystal structure of catechol O-methyltransferase.Nature368354358

25. Huotari M, Gogos JA, Karayiorgou M, Koponen O, Forsberg M, et al. (2002) Brain catecholamine metabolism in catechol-O-methyltransferase (COMT)-deficient mice. Eur J Neurosci 15: 246–256.M. HuotariJA GogosM. KarayiorgouO. KoponenM. Forsberg2002Brain catecholamine metabolism in catechol-O-methyltransferase (COMT)-deficient mice.Eur J Neurosci15246256

26. Mikolajczyk E, Grzywacz A, Samochowiec J (2010) The association of catechol-O-methyltransferase genotype with the phenotype of women with eating disorders. Brain Res 1307: 142–148.E. MikolajczykA. GrzywaczJ. Samochowiec2010The association of catechol-O-methyltransferase genotype with the phenotype of women with eating disorders.Brain Res1307142148

27. Frieling H, Romer KD, Wilhelm J, Hillemacher T, Kornhuber J, et al. (2006) Association of catecholamine-O-methyltransferase and 5-HTTLPR genotype with eating disorder-related behavior and attitudes in females with eating disorders. Psychiatr Genet 16: 205–208.H. FrielingKD RomerJ. WilhelmT. HillemacherJ. Kornhuber2006Association of catecholamine-O-methyltransferase and 5-HTTLPR genotype with eating disorder-related behavior and attitudes in females with eating disorders.Psychiatr Genet16205208

28. Pooley EC, Fineberg N, Harrison PJ (2007) The met(158) allele of catechol-O-methyltransferase (COMT) is associated with obsessive-compulsive disorder in men: case-control study and meta-analysis. Mol Psychiatry 12: 556–561.EC PooleyN. FinebergPJ Harrison2007The met(158) allele of catechol-O-methyltransferase (COMT) is associated with obsessive-compulsive disorder in men: case-control study and meta-analysis.Mol Psychiatry12556561

29. Kozak M (1984) Point mutations close to the AUG initiator codon affect the efficiency of translation of rat preproinsulin in vivo. Nature 308: 241–246.M. Kozak1984Point mutations close to the AUG initiator codon affect the efficiency of translation of rat preproinsulin in vivo.Nature308241246

30. Kozak M (1986) Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell 44: 283–292.M. Kozak1986Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes.Cell44283292

31. Kozak M (1987) An analysis of 5′-noncoding sequences from 699 vertebrate messenger RNAs. Nucleic Acids Res 15: 8125–8148.M. Kozak1987An analysis of 5′-noncoding sequences from 699 vertebrate messenger RNAs.Nucleic Acids Res1581258148

32. Bartel DP (2009) MicroRNAs: target recognition and regulatory functions. Cell 136: 215–233.DP Bartel2009MicroRNAs: target recognition and regulatory functions.Cell136215233

33. Tanzer A, Stadler PF (2004) Molecular evolution of a microRNA cluster. J Mol Biol 339: 327–335.A. TanzerPF Stadler2004Molecular evolution of a microRNA cluster.J Mol Biol339327335

34. Lee CT, Risom T, Strauss WM (2007) Evolutionary conservation of microRNA regulatory circuits: an examination of microRNA gene complexity and conserved microRNA-target interactions through metazoan phylogeny. DNA Cell Biol 26: 209–218.CT LeeT. RisomWM Strauss2007Evolutionary conservation of microRNA regulatory circuits: an examination of microRNA gene complexity and conserved microRNA-target interactions through metazoan phylogeny.DNA Cell Biol26209218

35. Flagel SB, Clark JJ, Robinson TE, Mayo L, Czuj A, et al. (2010) A selective role for dopamine in stimulus-reward learning. Nature 469: 53–57.SB FlagelJJ ClarkTE RobinsonL. MayoA. Czuj2010A selective role for dopamine in stimulus-reward learning.Nature4695357

36. Haider S, Ballester B, Smedley D, Zhang J, Rice P, et al. (2009) BioMart Central Portal--unified access to biological data. Nucleic Acids Res 37: W23–27.S. HaiderB. BallesterD. SmedleyJ. ZhangP. Rice2009BioMart Central Portal--unified access to biological data.Nucleic Acids Res37W2327

37. Smedley D, Haider S, Ballester B, Holland R, London D, et al. (2009) BioMart--biological queries made easy. BMC Genomics 10: 22.D. SmedleyS. HaiderB. BallesterR. HollandD. London2009BioMart--biological queries made easy.BMC Genomics1022

38. Katoh K, Toh H (2010) Parallelization of the MAFFT multiple sequence alignment program. Bioinformatics 26: 1899–1900.K. KatohH. Toh2010Parallelization of the MAFFT multiple sequence alignment program.Bioinformatics2618991900

39. Tabaska JE, Zhang MQ (1999) Detection of polyadenylation signals in human DNA sequences. Gene 231: 77–86.JE TabaskaMQ Zhang1999Detection of polyadenylation signals in human DNA sequences.Gene2317786

40. Friedman RC, Farh KK, Burge CB, Bartel DP (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105.RC FriedmanKK FarhCB BurgeDP Bartel2009Most mammalian mRNAs are conserved targets of microRNAs.Genome Res1992105

41. Lewis BP, Burge CB, Bartel DP (2005) Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120: 15–20.BP LewisCB BurgeDP Bartel2005Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.Cell1201520

42. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB (2003) Prediction of mammalian microRNA targets. Cell 115: 787–798.BP LewisIH ShihMW Jones-RhoadesDP BartelCB Burge2003Prediction of mammalian microRNA targets.Cell115787798

43. Zuker M (2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 31: 3406–3415.M. Zuker2003Mfold web server for nucleic acid folding and hybridization prediction.Nucleic Acids Res3134063415

44. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ (2008) miRBase: tools for microRNA genomics. Nucleic Acids Res 36: D154–158.S. Griffiths-JonesHK SainiS. van DongenAJ Enright2008miRBase: tools for microRNA genomics.Nucleic Acids Res36D154158

45. Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, et al. (2006) A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell 126: 1203–1217.KC MirandaT. HuynhY. TayYS AngWL Tam2006A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes.Cell12612031217