Research Article |
Corresponding author: Hairui Wang ( wang200108143@aliyun.com ) Academic editor: Jesus Maldonado
© 2020 Xin He, Walter H. Hsu, Rong Hou, Ying Yao, Qin Xu, Dandan Jiang, Longqiong Wang, Hairui Wang.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
He X, Hsu WH, Hou R, Yao Y, Xu Q, Jiang D, Wang L, Wang H (2020) Comparative genomics reveals bamboo feeding adaptability in the giant panda (Ailuropoda melanoleuca). ZooKeys 923: 141-156. https://doi.org/10.3897/zookeys.923.39665
|
The giant panda (Ailuropoda melanoleuca) is one of the world’s most endangered mammals and remains threatened as a result of intense environmental and anthropogenic pressure. The transformation and specialization of the giant panda’s diet into a herbivorous diet have resulted in unique adaptabilities in many aspects of their biology, physiology and behavior. However, little is known about their adaptability at the molecular level. Through comparative analysis of the giant panda’s genome with those of nine other mammalian species, we found some genetic characteristics of the giant panda that can be associated with adaptive changes for effective digestion of plant material. We also found that giant pandas have similar genetic characteristics to carnivores in terms of olfactory perception but have similar genetic characteristics to herbivores in terms of immunity and hydrolytic enzyme activity. Through the analysis of gene family expansion, 3752 gene families were found, which were enriched in functions such as digestion. A total of 93 genes under positive selection were screened out and gene enrichment identified these genes for the following processes: negative regulation of cellular metabolic process, negative regulation of nitrogen compound metabolic process, negative regulation of macromolecule metabolic process and negative regulation of metabolic process. Combined with the KEGG pathway, it was found that genes such as CREB3L1, CYP450 2S1, HSD11B2, LRPAP1 play a key role in digestion. These genes may have played a key role in the pandas’ adaptation to its bamboo diet.
adaptation, bamboo diet, dietary transition, digestion, feeding habits
Diet may be the most important selective force in animal evolution (
Bamboo is a low nutrition/energy food comprising of 70–80% cellulose, hemicellulose, and lignin and 20–30% protein, soluble carbohydrate, and fat (
The genome sequences of wild animal species are rapidly being accumulated, providing rich resources for the study of adaptation, trait evolution, species divergence, and population structure analyses (
Myr: million years; BLAST: basic local alignment search tool; GO: gene ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; NCBI: National Center for Biotechnology Information; dN/dS: non-synonymous/synonymous rate ratio; ML: maximum likelihood; CREB3L1: cAMP-responsive element binding protein 3 like 1; CYP450 2S1: cytochrome P450, family 2, subfamily s, polypeptide 1; HSD11B2: corticosteroid 11-beta-dehydrogenase isozyme 2; LRPAP1: low density lipoprotein receptor related protein associated protein 1
We performed comparative genomic analysis of 10 mammalian species with different diets including the giant panda (Ailuropoda melanoleuca); four species of mammals with carnivorous diets: wolf (Canis lupus familiaris), tiger (Panthera tiger) polar bear (Ursus maritimus) and sperm whale (Physeter catodon)); two species with herbivorous diets: white rhinoceros (Ceratotherium simum simum) and gorilla (Gorilla gorilla); and three species of omnivores: macaque (Macaca mulatta), human (Homo sapiens) and mouse (Mus musculus). For the present study, we downloaded the genome sequences, protein sequences and annotation files of the 10 species from the NCBI website (2018/3/7) (Table
OrthoFinder (OrthoFinder-2.2.7) software was used for the gene orthology analysis (
The single-copy gene protein sequences were compared using the MAFFT (v7.158b (2014/06/27)) software (
The giant panda gene was chosen as the foreground and the genes of the remaining species as the background. The positive selection analysis was used to determine whether there was a significant difference between the non-synonymous replacement rate and the proportion of synonymous replacement rates (dN/dS) between foreground and background branches. The null hypothesis parameters were: model = 2, NSsites = 2, fix_omega = 1, omega = 1. Alternative hypothesis parameters were: model = 2, NSsites = 2, fix_omega = 0. Chi programs using Paml software (
A total of 517,058 protein-coding genes from 10 species were used for the gene family analysis; 481,081 genes were identified in 24,788 gene families, including 911 single-copy true orthologous genes across all 10 species (Table
The enrichment analysis of shared genes between the giant panda and mammalian species with different feeding habits. Giant pandas have the characteristics of both carnivores and herbivores. Studies show that it is close to carnivores in perception and close to herbivores in physiological functions. The abscissa is the pair value of the corrected p value, and the corrected p < 0.05 is taken as the threshold value. a shared genes between the giant panda and other mammalian species with different feeding habits b gene enrichment analysis of the giant panda and carnivores c gene enrichment analysis of the giant panda and herbivores.
Property | Value |
---|---|
Number of genes | 517,058 |
Number of genes in orthogroups | 481,081 |
Number of unassigned genes | 35,977 |
Percentage of genes in orthogroups | 93.00% |
Percentage of unassigned genes | 7.00% |
Number of orthogroups | 24,788 |
Number of species-specific orthogroups | 169 |
Number of genes in species-specific orthogroups | 1108 |
Percentage of genes in species-specific orthogroups | 0.20% |
Mean orthogroup size | 19.4 |
Median orthogroup size | 15 |
Number of orthogroups with all species present | 14,680 |
Number of single-copy orthogroups | 911 |
The phylogenetic tree constructed using all 10 mammalian species based on single-copy gene family data estimated the time of divergence between giant panda and polar bear to be 21 Myr (Fig.
A total of 93 genes were identified to be under positive selection in the giant panda branch. In the GO enrichment analysis, it was found that the positive gene enrichment was selected for the negative regulation of the cellular metabolic process, negative regulation of the nitrogen compound metabolic process, negative regulation of the macromolecule metabolic process and negative regulation of the metabolic process (Fig.
The carnivorous diet of the ancestors of the giant panda has gradually evolved into a strict bamboo diet since the late Pliocene or early Pleistocene. In order to adapt to this change in diet, giant pandas developed adaptive evolution genetically, physiologically and morphologically. By comparing the genes that the giant panda shares with several carnivorous and herbivorous species, we found that they share genes with similar smell and perception functions (
Giant pandas have several families of amplified genes, among which salivary secretion, pancreatic secretion, insulin secretion and parathyroid hormone synthesis, secretion and function pathways may be important for the digestion and adaptation of giant pandas to bamboo. The sialaden and pancreas are important glands for the digestive system; the increased saliva secretion can help giant pandas lubricate the gut for digesting the starch in bamboo (
In the present study, 93 genes were found under positive selection in the genome of the giant panda. The GO enrichment analysis detected that these genes were concentrated in the negative regulation mechanism of metabolism. Activation of this mechanism could reduce their own metabolism and energy needed for adaptation to ecological changes. The results of the present study are consistent with the report showing that in order to cope with the low energy and nutrition content of the bamboo diet, giant pandas reduce their own metabolism, reduce behavioral activities, and eat voraciously to meet their own energy demand (
Other genes hint at adaptations in other directions. For example, the HSD11B2 gene can encode type 2 11β-hydroxysteroid dehydrogenase and participate in intracellular homeostasis, and convert cortisol to cortisone. It is an inactive corticosteroid that can prevent obesity and high blood pressure (
Changes in diet composition during animal evolution likely created strong selective pressures in multiple biological processes. The giant pandas’ dietary habits have shifted dramatically from a carnivorous diet to a strict bamboo diet. The reasons for the changes in the diet of giant pandas are still largely unknown. With the development of high throughput, next generation sequencing technology and the open availability of genome-wide data, we explored the adaptability of giant pandas from a genomic perspective.
In this study, the comparative genomics of species with different feeding habits was used to elucidate the genes responsible for the adaptation of giant pandas to herbivory. The giant panda retains the characteristics of a carnivore in its sensory abilities, but is similar to a herbivore in terms of physiology, such as digestion. These adaptations also help the panda digest bamboo and metabolize nutrients by regulating its digestive and endocrine secretions. Here, we elucidated the molecular/genomic mechanism of giant pandas’ adaptation to dietary changes and the adoption of dietary specializations in the digestive system. However, the function of these genes in giant pandas need to be further verified. In response to changes in diet, the giant panda’s genome has undergone important evolutionary adaptations that help it better digest bamboo and metabolize nutrients.
This study also provides new perspectives and insights for the adaptive evolution of giant pandas. In addition, these findings also provide the molecular theoretical basis for the adaptation mechanism of giant panda’s dietary changes, and present additional information that can be used for the management and protection of this endangered species.
The authors wish to thank the Shanghai Personal Biotechnology Co., Ltd for the bioinformatics analysis support provided.
This work was supported by Sichuan Science and Technology Program (2018JY0349) and the Program of the Giant Panda Breeding Research Foundation (CPF2017-10).