Giant pandas are not an evolutionary cul-de-sac: evidence from multidisciplinary research

Conservation Genomics and Metagenomics of Giant and Red Pandas in the Wild

Giant pandas and red pandas are endangered species with similar specialized bamboo diet and partial sympatric distribution in China. Over the last two decades, the rapid development of genomics and metagenomics research on these species has enriched our knowledge of their biology, ecology, physiology, genetics, and evolution, which is crucial and useful for their conservation. We describe the evolutionary history, endangerment processes, genetic diversity, and population structure of wild giant pandas and two species of red pandas (Chinese and Himalayan red pandas). In addition, we explore how genomics and metagenomics studies have provided insight into the convergent adaptation of pandas to the specialized bamboo diet. Finally, we discuss how these findings are applied to effective conservation management of giant and red pandas in the wild and in captivity to promote the long-term persistence of these species.

Forest fragmentation causes an isolated population of the golden takin (Budorcas taxicolor bedfordi Thomas, 1911) (Artiodactyla: Bovidae) in the Qinling Mountains (China)

Budorcas taxicolor bedfordi is a rare animal uniquely distributed in the Qinling Mountains (China). Human disturbance and habitat fragmentation have directly affected the survival of B. t. bedfordi. It is urgent to clarify the genetic diversity and genetic structure of the B. t. bedfordi population and implement effective conservation measures. In this study, 20 new polymorphic microsatellite loci were isolated by Illumina sequencing. The genetic diversity and population structure of 124 B. t. bedfordi individuals from three populations (Niubeliang population, Zhouzhi population, and Foping population) were analysed according to these 20 microsatellite loci. Our results indicated that B. t. bedfordi had a low level of genetic variability and that there was inbreeding in the three populations. The population genetic structure analyses showed that the Niubeliang population had a trend of differentiation from other populations. National roads can affect population dispersal, while ecological corridors can promote population gene exchange. None of the three B. t. bedfordi populations experienced bottleneck effects. For conservation management plans, the Zhouzhi population and Foping population should be considered one management unit, and the Niubeliang population should be considered another management unit. We suggest building an ecological corridor to keep the habitat connected and formulating tourism management measures to reduce the influence of human disturbance on B. t. bedfordi.

Single-nucleus transcriptome inventory of giant panda reveals cellular basis for fitness optimization under low metabolism

BACKGROUND

Energy homeostasis is essential for the adaptation of animals to their environment and some wild animals keep low metabolism adaptive to their low-nutrient dietary supply. Giant panda is such a typical low-metabolic mammal exhibiting species specialization of extremely low daily energy expenditure. It has low levels of basal metabolic rate, thyroid hormone, and physical activities, whereas the cellular bases of its low metabolic adaptation remain rarely explored.

RESULTS

In this study, we generate a single-nucleus transcriptome atlas of 21 organs/tissues from a female giant panda. We focused on the central metabolic organ (liver) and dissected cellular metabolic status by cross-species comparison. Adaptive expression mode (i.e., AMPK related) was prominently displayed in the hepatocyte of giant panda. In the highest energy-consuming organ, the heart, we found a possibly optimized utilization of fatty acid. Detailed cell subtype annotation of endothelial cells showed the uterine-specific deficiency of blood vascular subclasses, indicating a potential adaptation for a low reproductive energy expenditure.

CONCLUSIONS

Our findings shed light on the possible cellular basis and transcriptomic regulatory clues for the low metabolism in giant pandas and helped to understand physiological adaptation response to nutrient stress.

PandaGUT provides new insights into bacterial diversity, function, and resistome landscapes with implications for conservation

BACKGROUND

The gut microbiota play important roles in host adaptation and evolution, but are understudied in natural population of wild mammals. To address host adaptive evolution and improve conservation efforts of threatened mammals from a metagenomic perspective, we established a high-quality gut microbiome catalog of the giant panda (pandaGUT) to resolve the microbiome diversity, functional, and resistome landscapes using approximately 7 Tbp of long- and short-read sequencing data from 439 stool samples.

RESULTS

The pandaGUT catalog comprises 820 metagenome-assembled genomes, including 40 complete closed genomes, and 64.5% of which belong to species that have not been previously reported, greatly expanding the coverage of most prokaryotic lineages. The catalog contains 2.37 million unique genes, with 74.8% possessing complete open read frames, facilitating future mining of microbial functional potential. We identified three microbial enterotypes across wild and captive panda populations characterized by Clostridium, Pseudomonas, and Escherichia, respectively. We found that wild pandas exhibited host genetic-specific microbial structures and functions, suggesting host-gut microbiota phylosymbiosis, while the captive cohorts encoded more multi-drug resistance genes.

CONCLUSIONS

Our study provides largely untapped resources for biochemical and biotechnological applications as well as potential intervention avenues via the rational manipulation of microbial diversity and reducing antibiotic usage for future conservation management of wildlife. Video Abstract.

Signal detection theory applied to giant pandas: Do pandas go out of their way to make sure their scent marks are found?

Inter-animal communication allows signals released by an animal to be perceived by others. Scent-marking is the primary mode of such communication in giant pandas (Ailuropoda melanoleuca). Signal detection theory propounds that animals choose the substrate and location of their scent marks so that the signals released are transmitted more widely and last longer. We believe that pandas trade-off scent-marking because they are an energetically marginal species and it is costly to generate and mark chemical signals. Existing studies only indicate where pandas mark more frequently, but their selection preferences remain unknown. This study investigates whether the marking behavior of pandas is consistent with signal detection theory. Feces count, reflecting habitat use intensity, was combined with mark count to determine the selection preference for marking. The results showed that pandas preferred to mark ridges with animal trails and that most marked tree species were locally dominant. In addition, marked plots and species were selected for lower energy consumption and a higher chance of being detected. Over 90% of the marks used were the longest-surviving anogenital gland secretion marks, and over 80% of the marks were oriented toward animal trails. Our research demonstrates that pandas go out of their way to make sure their marks are found. This study not only sheds light on the mechanisms of scent-marking by pandas but also guides us toward more precise conservation of the panda habitat.

The unique gut microbiome of giant pandas involved in protein metabolism contributes to the host's dietary adaption to bamboo

BACKGROUND

The gut microbiota of the giant panda (Ailuropoda melanoleuca), a global symbol of conservation, are believed to be involved in the host's dietary switch to a fibrous bamboo diet. However, their exact roles are still largely unknown.

RESULTS

In this study, we first comprehensively analyzed a large number of gut metagenomes giant pandas (n = 322), including 98 pandas sequenced in this study with deep sequencing (Illumina) and third-generation sequencing (nanopore). We reconstructed 408 metagenome-assembled genomes (MAGs), and 148 of which (36.27%) were near complete. The most abundant MAG was classified as Streptococcus alactolyticus. A pairwise comparison of the metagenomes and meta-transcriptomes in 14 feces revealed genes involved in carbohydrate metabolism were lower, but those involved in protein metabolism were greater in abundance and expression in giant pandas compared to those in herbivores and omnivores. Of note, S. alactolyticus was positively correlated to the KEGG modules of essential amino-acid biosynthesis. After being isolated from pandas and gavaged to mice, S. alactolyticus significantly increased the relative abundance of essential amino acids in mice jejunum.

CONCLUSIONS

The study highlights the unique protein metabolic profiles in the giant panda's gut microbiome. The findings suggest that S. alactolyticus is an important player in the gut microbiota that contributes to the giant panda's dietary adaptation by more involvement in protein rather than carbohydrate metabolism. Video Abstract.

Fecal Metabolomics Reveals the Foraging Strategies of Giant Pandas for Different Parts of Bamboo

Climate change-induced food shortages pose major threats to wildlife conservation, and the exclusive reliance of giant pandas on bamboo makes them particularly vulnerable. The aim of this study was to provide insight into the reasons for the foraging strategies of giant pandas to selectively forage for different bamboo parts (bamboo shoot, culm, and leaf) during different seasons. This study used a metabolomic approach to analyze the fecal metabolites of giant pandas and conducted a correlation analysis with their gut microbiota. The results indicate that the fecal metabolites of giant pandas differ significantly depending on the bamboo parts they forage on, with higher sugar content observed when they consume bamboo culm with high fiber content. By functional annotation, culm group metabolites were enriched in the galactose metabolic pathway, while shoot group metabolites were enriched in the phenylalanine, tyrosine and tryptophan biosynthesis pathways. Moreover, Streptococcus showed a significant positive correlation with glucose and acetic acid content. Therefore, the foraging strategy of giant pandas is based on the ability to utilize the nutrient content of different bamboo parts. Captive feeding and habitat construction should enrich bamboo species to allow them to express their natural foraging strategies and improve their welfare and reproductive status.

Feeding habits and foraging patch selection strategy of the giant panda in the Meigu Dafengding National Nature Reserve, Sichuan Province, China

Despite retaining a digestive system with carnivorous characteristics, the giant panda has now specialized to eat bamboo. Therefore, research on its feeding habits and foraging strategies has always been an important field in the study of giant panda population ecology. Located in the core distribution area of the giant panda population in the Liangshan mountains, the Meigu Dafengding Nature Reserve acts as the link between the major panda reserves in the Liangshan mountains and the key corridor zone connecting the giant panda populations in the east and west. This study was performed in 2017 and aimed to determine the feeding habits of giant pandas in this region, the heterogeneity of food resources in different seasons, and the selection strategy for foraging patches. To achieve this, we used the line transect method to investigate the giant panda population, collected fresh feces, analyzed the nutritional composition of bamboo samples, and set up a sampling grid of foraging patches in different seasonal distribution areas. It was found that wild giant pandas in the region spend most of the year at lower elevations(1939-3296 m) foraging for Yushania maculata and Yushania ailuropodina, and only move to higher elevations(2844-3770 m) in summer for Bashania fangiana. Their preferred foraging patches(n = 65) had certain topographical features, such as gentle slopes(< 25°) and closer proximity to water sources(< 500 m), and the nutritional quality of bamboo resources in foraging patches was significantly better than in those they avoided or rarely used. The food habits and foraging patch selection strategies of giant pandas in the Liangshan and Qinling mountains are relatively similar, and the results of this study can help provide a scientific basis for the conservation and management of wild giant panda habitats in different mountain regions.

The anthropogenic effect of land use on population genetics of Malcus inconspicuus

Since the beginning of the Holocene era, human activities have seriously impacted animal habitats and vegetative environments. Species that are dependent on natural habitats or with narrow niches might be more severely affected by habitat changes. Malcus inconspicuus is distributed in subtropical China and highly dependent on the mountain environment. Our study investigated the role of the mountainous landscape in the historical evolution of M. inconspicuus and the impact of Holocene human activities on it. A phylogeographical approach was implemented with integrative datasets including double-digest restriction site-associated DNA (ddRAD), mitochondrial data, and distribution data. Three obvious clades and an east-west phylogeographical pattern were found in subtropical China. Mountainous landscape has "multifaceted" effects on the evolutionary history of M. inconspicuus, it has contributed to population differentiation, provided glacial refuges, and provided population expansion corridors during the postglacial period. The effective population size (Ne) of M. inconspicuus showed a sharp decline during the Holocene era, which revealed a significantly negative correlation with the development of cropland in a hilly area at the same time and space. It supported that the species which are highly dependent on natural habitats might undergo greater impact when the habitat was damaged by agricultural activities and we should pay more attention to them, especially in the land development of their distribution areas.

Using Blood Transcriptome Analysis to Determine the Changes in Immunity and Metabolism of Giant Pandas with Age

A low reproductive rate coupled with human activities has endangered the giant panda, a species endemic to southwest China. Although giant pandas feed almost exclusively on bamboo, they retain carnivorous traits and suffer from carnivorous diseases. Additionally, their immune system is susceptible to aging, resulting in a reduced ability to respond to diseases. This study aimed to determine the genes and pathways expressed differentially with age in blood tissues. The differentially expressed genes in different age groups of giant pandas were identified by RNA-seq. The elderly giant pandas had many differentially expressed genes compared with the young group (3 years old), including 548 upregulated genes and 401 downregulated genes. Further, functional enrichment revealed that innate immune upregulation and adaptive immune downregulation were observed in the elderly giant pandas compared with the young giant pandas. Meanwhile, the immune genes in the elderly giant pandas changed considerably, including genes involved in innate immunity and adaptive immunity such as PLSCR1, CLEC7A, CCL5, CCR9, and EPAS1. Time series analysis found that giant pandas store glycogen by prioritizing fat metabolism at age 11, verifying changes in the immune system. The results reported in this study will provide a foundation for further research on disease prevention and the energy metabolism of giant pandas.

Specific Gene Duplication and Loss of Cytochrome P450 in Families 1-3 in Carnivora (Mammalia, Laurasiatheria)

Simple Summary

In this study we investigated the specific duplication and loss events of cytochrome P450 (CYP) genes in families 1-3 in Carnivora. These genes have been recognized as essential detoxification enzymes, and, using genomic data, we demonstrated a synteny analysis of the CYP coding cluster and a phylogenetic analysis of these genes. We discovered the CYP2Cs and CYP3As expansion in omnivorous species such as the badger, the brown bear, the black bear, and the dog. Furthermore, phylogenetic analysis revealed the evolution of CYP2Cs and 3As in Carnivora. These findings are essential for the appropriate estimation of pharmacokinetics or toxicokinetic in wild carnivorans.

Abstract

Cytochrome P450s are among the most important xenobiotic metabolism enzymes that catalyze the metabolism of a wide range of chemicals. Through duplication and loss events, CYPs have created their original feature of detoxification in each mammal. We performed a comprehensive genomic analysis to reveal the evolutionary features of the main xenobiotic metabolizing family: the CYP1-3 families in Carnivora. We found specific gene expansion of CYP2Cs and CYP3As in omnivorous animals, such as the brown bear, the black bear, the dog, and the badger, revealing their daily phytochemical intake as providing the causes of their evolutionary adaptation. Further phylogenetic analysis of CYP2Cs revealed Carnivora CYP2Cs were divided into CYP2C21, 2C41, and 2C23 orthologs. Additionally, CYP3As phylogeny also revealed the 3As’ evolution was completely different to that of the Caniformia and Feliformia taxa. These studies provide us with fundamental genetic and evolutionary information on CYPs in Carnivora, which is essential for the appropriate interpretation and extrapolation of pharmacokinetics or toxicokinetic data from experimental mammals to wild Carnivora.

Giant pandas' staple food bamboo phyllosphere fungal community and its influencing factors

Giant pandas have developed a series of foraging strategies to adapt to their special bamboo diets. Although bamboo is an important food resource for giant pandas in Liziping National Nature Reserve (Liziping NR), China, there are relatively few studies on their phyllosphere fungal community and its influencing factors. Herein, we used ITS1 amplification and metagenomic sequencing to analyze the phyllosphere fungi diversity and functions (KEGG, CAZyme, and antibiotic resistance gene) and explore the influencing factors for the three giant pandas foraging bamboo species (Arundinaria spanostachya, AS; Yushania lineolate, YL; and Fargesia ferax, FF) over different seasons (spring vs. autumn) in Liziping NR, China. We found that Ascomycota and Basidiomycota were the most dominant phyla in the bamboo phyllosphere. The alpha diversity (e.g., the Sobs index and Shannon index) was relatively higher in autumn samples than in spring samples, and the community structure differed significantly between the three bamboo species in spring and autumn. Some biotic and abiotic variables (e.g., the elevation and mean base diameter of bamboo) significantly influenced the abundance, diversity, and community structure of the bamboo phyllosphere fungal community. Moreover, the functional analysis showed the differences in the glycoside hydrolase community and antibiotic resistance gene (ARG) profile between spring and autumn samples. Co-occurrence network modeling suggested that AS phyllosphere fungal communities in autumn employed a much more complex network than that in spring, and the abundance of multidrug, tetracycline, and glycopeptide resistance genes was high and closely correlated with other ARGs. These results indicate that fungal community's abundance, diversity, and community structure are mainly affected by the season, host species, and elevation. The season and host species are major factors affecting the biological functions (KEGG and CAZyme), ARGs, and interactions between sympatric bacterial and fungal communities in bamboo phyllosphere. This integrated study can provide a reference basis for the seasonal management of bamboo resources foraged by wild giant pandas, and predict the risk of antibiotic resistance in bamboo phyllosphere fungal flora in Liziping NR (Xiaoxiangling mountains), China.

Adaptation of gut microbiome and host metabolic systems to lignocellulosic degradation in bamboo rats

Bamboo rats (Rhizomys pruinosus) are among the few mammals that lives on a bamboo-based diet which is mainly composed of lignocellulose. However, the mechanisms of adaptation of their gut microbiome and metabolic systems in the degradation of lignocellulose are largely unknown. Here, we conducted a multi-omics analysis on bamboo rats to investigate the interaction between their gut microbiomes and metabolic systems in the pre- and post-weaning periods, and observed significant relationships between dietary types, gut microbiome, serum metabolome and host gene expression. For comparison, published gut microbial data from the famous bamboo-eating giant panda (Ailuropoda melanoleuca) were also used for analysis. We found that the adaptation of the gut microbiome of the bamboo rat to a lignocellulose diet is related to a member switch in the order Bacteroidales from family Bacteroidaceae to family Muribaculaceae, while for the famous bamboo-eating giant panda, several aerobes and facultative anaerobes increase after weaning. The conversion of bacteria with an increased relative abundance in bamboo rats after weaning enriched diverse carbohydrate-active enzymes (CAZymes) associated with lignocellulose degradation and functionally enhanced the biosynthesis of amino acids and B vitamins. Meanwhile, the circulating concentration of short-chain fatty acids (SCFAs) derived metabolites and the metabolic capacity of linoleic acid in the host were significantly elevated. Our findings suggest that fatty acid metabolism, including linoleic acid and SCFAs, are the main energy sources for bamboo rats in response to the low-nutrient bamboo diet.

Discussions of the "Not So Fit": How Ableism Limits Diverse Thought and Investigative Potential in Evolutionary Biology

AbstractEvolutionary biology and many of its foundational concepts are grounded in a history of ableism and eugenics. The field has not made a concerted effort to divest our concepts and investigative tools from this fraught history, and as a result, an ableist investigative lens has persisted in present-day evolutionary research, limiting the scope of research and harming the ability to communicate and synthesize knowledge about evolutionary processes. This failure to divest from our eugenicist and ableist history has harmed progress in evolutionary biology and allowed principles from evolutionary biology to continue to be weaponized against marginalized communities in the modern day. To rectify this problem, scholars in evolutionary research must come to terms with how the history of the field has influenced their investigations and work to establish a new framework for defining and investigating concepts such as selection and fitness.

Earliest giant panda false thumb suggests conflicting demands for locomotion and feeding

Of the many peculiarities that enable the giant panda (Ailuropoda melanoleuca), a member of the order Carnivora, to adapt to life as a dedicated bamboo feeder, its extra “thumb” is arguably the most celebrated yet enigmatic. In addition to the normal five digits in the hands of most mammals, the giant panda has a greatly enlarged wrist bone, the radial sesamoid, that acts as a sixth digit, an opposable “thumb” for manipulating bamboo. We report the earliest enlarged radial sesamoid, already a functional opposable “thumb,” in the ancestral panda Ailurarctos from the late Miocene site of Shuitangba in Yunnan Province, China. However, since the late Miocene, the “thumb” has not enlarged further because it must be balanced with the constraints of weight bearing while walking in a plantigrade posture. This morphological adaptation in panda evolution thus reflects a dual function of the radial sesamoid for both bamboo manipulation and weight distribution. The latter constraint could be the main reason why the panda’s false thumb never evolved into a full digit. This crude “thumb” suggests that the origin of the panda’s dedicated bamboo diet goes back to as early as 6–7 Ma.

Transcriptome sequencing of black and white hair follicles in the giant panda

With the completion of the draft assembly of the giant panda genome sequence, RNA sequencing technology has been widely used in genetic research on giant pandas. We used RNA-seq to examine black and white hair follicle samples from adult pandas. By comparison with the giant panda genome, 75 963 SNP loci were labeled, 2 426 differentially expressed genes were identified, and 2 029 new genes were discovered, among which 631 were functionally annotated. A cluster analysis of the differentially expressed genes showed that they were mainly related to the Wnt signaling pathway, ECM-receptor interaction, the p53 signaling pathway and ribosome processing. The enrichment results showed that there were significant differences in the regulatory networks of hair follicles with different colors during the transitional stage of hair follicle resting growth, which may play a regulatory role in melanin synthesis during growth. In conclusion, our results provide new insights and more data support for research on the color formation in giant pandas. This article is protected by copyright. All rights reserved.

Genome of the Giant Panda Roundworm Illuminates Its Host Shift and Parasitic Adaptation

Baylisascaris schroederi, a roundworm (ascaridoid) parasite specific to the bamboo-feeding giant panda (Ailuropoda melanoleuca), represents a leading cause of mortality in wild giant panda populations. Here, we present a 293-megabase chromosome-level genome assembly of B. schroederi to infer its biology, including host adaptations. Comparative genomics revealed an evolutionary trajectory accompanied by host-shift events in ascaridoid parasite lineages after host separations, suggesting their potential for transmission and rapid adaptation to new hosts. Genomic and anatomical lines of evidence, including expansion and positive selection of genes related to the cuticle and basal metabolisms, indicate that B. schroederi undergoes specific adaptations to survive in the sharp-edged bamboo-enriched gut of giant pandas by structurally increasing its cuticle thickness and efficiently utilizing host nutrients through gut parasitism. Additionally, we characterized the secretome of B. schroederi and predicted potential drug and vaccine targets for new control strategies. Overall, this genome resource provides new insights into the host adaptation of B. schroederi to the giant panda as well as the host-shift events in ascaridoid parasite lineages. Our findings on the unique biology of B. schroederi will also aid in the development of prevention and treatment measures to protect giant panda populations from roundworm parasitism.

Seasonal shift of the gut microbiome synchronizes host peripheral circadian rhythm for physiological adaptation to a low-fat diet in the giant panda

Characteristics of the gut microbiome vary synchronously with changes in host diet. However, the underlying effects of these fluctuations remain unclear. Here, we performed fecal microbiota transplantation (FMT) of diet-specific feces from an endangered mammal (the giant panda) into a germ-free mouse model. We demonstrated that the butyrate-producing bacterium Clostridium butyricum was more abundant during shoot-eating season than during the leaf-eating season, congruent with the significant increase in host body mass. Following season-specific FMT, the microbiota of the mouse model resembled that of the donor, and mice transplanted with the microbiota from the shoot-eating season grew faster and stored more fat. Mechanistic investigations revealed that butyrate extended the upregulation of hepatic circadian gene Per2, subsequently increasing phospholipid biosynthesis. Validation experiments further confirmed this causal relationship. This study demonstrated that seasonal shifts in the gut microbiome affect growth performance, facilitating a deeper understanding of host-microbe interactions in wild mammals.

Contrasting Vaginal Bacterial Communities Between Estrus and Non-estrus of Giant Pandas (Ailuropoda melanoleuca)

Bacterial infection and imbalance of bacterial community in the genitourinary system of giant panda could affect the reproductive health. In severe cases, it can also lead to abortion. In this study, 13 of vaginal secretions in the estrue (E) group and seven of vaginal secretions in the non-estrue (NE) group were used to study the composition and diversity of vaginal bacterial communities between estrus and non-estrus by 16S rRNA gene sequencing analysis. The results showed that the vaginal microbiome in giant pandas shared the same top five abundant species between estrus and non-estrus at the phylum level. However, the vaginal microbiome changed significantly during estrus at the genus level. In top 10 genera, the abundance of Escherichia, Streptococcus, and Bacteroides in the E group was significantly higher than that in the NE group (p<0.05); Azomonas, Porphyromonas, Prevotella, Campylobacter, and Peptoniphilus in the NE group was significantly higher than that in the E group (p<0.05). The richness and diversity of vaginal microbiome in giant panda on estrus were significantly lower than those on non-estrus (p<0.05). It is noteworthy that the abundance of Streptococcus, Escherichia, and Bacteroides of vagina in giant pandas maintained low abundance in the daily. Whereas, they increased significantly during estrus period, which may play an important role in female giant pandas during estrus period. It was hypothesized that hormones may be responsible for the changes in the vaginal microbiome of giant pandas between estrus and no-estrus stages.

Multi-omics reveals the positive leverage of plant secondary metabolites on the gut microbiota in a non-model mammal

BACKGROUND

Flavonoids are important plant secondary metabolites (PSMs) that have been widely used for their health-promoting effects. However, little is known about overall flavonoid metabolism and the interactive effects between flavonoids and the gut microbiota. The flavonoid-rich bamboo and the giant panda provide an ideal system to bridge this gap.

RESULTS

Here, integrating metabolomic and metagenomic approaches, and in vitro culture experiment, we identified 97 flavonoids in bamboo and most of them have not been identified previously; the utilization of more than 70% flavonoid monomers was attributed to gut microbiota; the variation of flavonoid in bamboo leaves and shoots shaped the seasonal microbial fluctuation. The greater the flavonoid content in the diet was, the lower microbial diversity and virulence factor, but the more cellulose-degrading species.

CONCLUSIONS

Our study shows an unprecedented landscape of beneficial PSMs in a non-model mammal and reveals that PSMs remodel the gut microbiota conferring host adaptation to diet transition in an ecological context, providing a novel insight into host-microbe interaction. Video abstract.

Preliminary investigation on iodine nutrition in captive giant pandas

BACKGROUND/OBJECTIVE

The giant panda belongs to the family Ursidae and, as a species of bear, still retains the simple digestive system of a Carnivoran. However, under the pressure of a specific habitat they had to adapt to a plant mono-diet consisting of bamboo with different species and growth stages around the year. A plant-based diet has relatively low iodine content with risk of iodine deficiency. Furthermore, bamboo contains cyanogenic glycosides releasing cyanide whose detoxification metabolite the thiocyanate acts as antagonist against iodine uptake and storage in the thyroid. To date very little is known about the iodine nutritional status of the giant panda, thus this study was conducted to receive the first information about the iodine nutrition of captive giant panda.

SUBJECTS/METHODS

Here we investigated the iodine content of bamboo with different plant parts/vegetation stage and species and further compounds of the captive giant panda diet. Next, the urinary iodine (UI) and urinary thiocyanate (UT) levels of infant, sub-adult, adult and geriatric captive giant pandas was measured during the periods when the pandas consume both bamboo leaves- and culm (bamboo leaf-culm stage). Afterwards, the UI of 19 adult giant pandas was measured again for the different iodine intake during bamboo shoot stage. Finally, in this study part also the fecal iodine concentration was analyzed for calculation of total iodine excretion in relation to the iodine intake.

RESULTS

Bamboo leaves had the highest iodine content (453 μg/kg dry matter (DM)), followed by the shoots (84 μg/kg DM, p <  0.05), while bamboo culm had the lowest value (12 μg/kg DM, p <  0.05). During bamboo leaf-culm stage, giant pandas of different age groups had different UI and UT levels (p <  0.05). Furthermore, UI and UT were positively correlated among sub-adult, adult and geriatric giant pandas (p <  0.05). In adult giant pandas during bamboo shoot stage, the iodine excretion in feces was not different from that in urine while their total iodine excretion was less than their iodine intake (p <  0.05). Moreover, during bamboo shoot stage, the UI level of adult giant pandas was much lower than noted during bamboo leaf-culm stage (p <  0.05).

CONCLUSIONS

Our results indicate that UI of captive giant pandas was related to their age as well as to the vegetation stage/part of bamboo they consumed reflecting a different periodic iodine supply. Thiocyanate and fecal excretion should be emphasized when considering the iodine nutrition of giant pandas.

A single nucleotide mutation in the dual-oxidase 2 (DUOX2) gene causes some of the panda's unique metabolic phenotypes

The giant panda (Ailuropoda melanoleuca) is an iconic bear native to China, famous for eating almost exclusively bamboo. This unusual dietary behavior for a carnivore is enabled by several key adaptations including low physical activity, reduced organ sizes and hypothyroidism leading to lowered energy expenditure. These adaptive phenotypes have been hypothesized to arise from a panda-unique single-nucleotide mutation in the dual-oxidase 2 (DUOX2) gene, involved in thyroid hormone synthesis. To test this hypothesis, we created genome-edited mice carrying the same point mutation as the panda and investigated its effect on metabolic phenotype. Homozygous mice were 27% smaller than heterozygous and wild-type ones, had 13% lower body mass-adjusted food intake, 55% decreased physical activity, lower mass of kidneys (11%) and brain (5%), lower serum thyroxine (T4: 36%), decreased absolute (12%) and mass-adjusted (5%) daily energy expenditure, and altered gut microbiota. Supplementation with T4 reversed the effects of the mutation. This work uses a state-of-the-art genome editing approach to demonstrate the link between a single-nucleotide mutation in a key endocrine-related gene and profound adaptive changes in the metabolic phenotype, with great importance in ecology and evolution.

Temporomandibular joint and Giant Panda's (Ailuropoda melanoleuca) adaptation to bamboo diet

Here, we present new evidence that evolutionary adaptation of the Ailuripodinae lineage to bamboo diet has taken place by morphological adaptations in the masticatory system. The giant panda in the wild and in captivity removes without an exception the outer skin of all bamboo shoots, rich in abrasive and toxic compounds, by the highly adapted premolars P3 and P4. The temporomandibular joint (TMJ) allows sidewise movement of the jaw and the premolars can, in a cusp-to-cusp position, remove the poorly digestible outer skin of the bamboo before crushing the bamboo with molars. Based on the evidence presented here, we suggest that adaptation of TMJ to lateral movement for enabling cusp-to-cusp contact of premolars is the crucial evolutionary factor as which we consider the key to understand the Ailuropodinae lineage adaptive pathway to utilize the bamboo resource.

Genomic Signatures of Coevolution between Nonmodel Mammals and Parasitic Roundworms

Antagonistic coevolution between host and parasite drives species evolution. However, most of the studies only focus on parasitism adaptation and do not explore the coevolution mechanisms from the perspective of both host and parasite. Here, through the de novo sequencing and assembly of the genomes of giant panda roundworm, red panda roundworm, and lion roundworm parasitic on tiger, we investigated the genomic mechanisms of coevolution between nonmodel mammals and their parasitic roundworms and those of roundworm parasitism in general. The genome-wide phylogeny revealed that these parasitic roundworms have not phylogenetically coevolved with their hosts. The CTSZ and prolyl 4-hydroxylase subunit beta (P4HB) immunoregulatory proteins played a central role in protein interaction between mammals and parasitic roundworms. The gene tree comparison identified that seven pairs of interactive proteins had consistent phylogenetic topology, suggesting their coevolution during host–parasite interaction. These coevolutionary proteins were particularly relevant to immune response. In addition, we found that the roundworms of both pandas exhibited higher proportions of metallopeptidase genes, and some positively selected genes were highly related to their larvae’s fast development. Our findings provide novel insights into the genetic mechanisms of coevolution between nonmodel mammals and parasites and offer the valuable genomic resources for scientific ascariasis prevention in both pandas.

Diet, habitat environment and lifestyle conversion affect the gut microbiomes of giant pandas

Gut microbiota (GM) are important for the health of giant pandas (GPs), in addition to the utilization of bamboo in their diets. However, it is not fully understood how diet, habitat environment and lifestyle contribute to the composition of GM in GP. Consequently, we evaluated how dietary changes, habitat environment conversions and lifestyle shifts influence the GM of GPs using high-throughput sequencing and genome-resolved metagenomics. The GM of GPs were more similar when their hosts exhibited the same diet. High fiber diets significantly increased the diversity and decreased the richness of gut bacterial communities alone or interacted with the age factor (p < 0.05). The abundances of Streptococcus, Pseudomonas, Enterococcus, Lactococcus, Acinetobacter, and Clostridium significantly increased during diet conversion process (Non-parametric factorial Kruskal-Wallis sum-rank test, LDA > 4). Reconstruction of 60 metagenome-assembled-genomes (MAGs) indicated that these bacteria were likely responsible for bamboo digestion via gene complements involved in cellulose, hemicellulose, and lignin degradation. While habitat environment may play a more important role in shaping the GM of GP, lifestyle can also greatly affect bacterial communities. The GM structure in reintroduced GPs notably converged to that of wild pandas. Importantly, the main bacterial genera of wild GPs could aid in lignin degradation, while those of reintroduced GPs were related to cellulose and hemicellulose digestion. Streptococcus, Pseudomonas, Enterococcus, Lactococcus, Acinetobacter, and Clostridium may contribute to lignocellulose digestion in GP. The results revealed that diet conversion, habitat environment and lifestyle could remarkably influence the GM of GP. In addition, results suggested that increasing the ability of lignin degradation with GM may aid to change the GM of reintroduced pandas to resemble those of wild pandas.

Laterality in Responses to Acoustic Stimuli in Giant Pandas

Simple Summary

Functional lateralization in the auditory system has been widely studied. Accordingly, behavioral laterality responses affected by acoustic stimuli have been observed in many vertebrate species. In this study, we assessed giant pandas’ behavioral responses to different acoustic stimuli in order to examine cerebral lateralization. We concluded that adult giant pandas showed a left-hemisphere bias in response to positive acoustic stimuli. Furthermore, we found the specific valence of cerebral lateralization for different categories of acoustic stimuli, of which some were relevant to the lateralization while others were not relevant. Our findings support an evolutionary strategy that giant pandas process auditory signals similar to other mammals.

Abstract

Cerebral lateralization is a common feature present in many vertebrates and is often observed in response to various sensory stimuli. Numerous studies have proposed that some vertebrate species have a right hemisphere or left hemisphere dominance in response to specific types of acoustic stimuli. We investigated lateralization of eight giant pandas (Ailuropoda melanoleuca) by using a head turning paradigm and twenty-eight acoustic stimuli with different emotional valences which included twenty-four conspecific and four non-conspecific acoustic stimuli (white noise, thunder, and vocalization of a predator). There was no significant difference in auditory laterality in responses to conspecific or non-conspecific sounds. However, the left cerebral hemisphere processed the positive stimuli, whereas neither of the two hemispheres exhibited a preference for processing the negative stimuli. Furthermore, the right hemisphere was faster than the left hemisphere in processing emotional stimuli and conspecific stimuli. These findings demonstrate that giant pandas exhibit lateralization in response to different acoustic stimuli, which provides evidence of hemispheric asymmetry in this species.

Diversity and Composition of Gut Bacterial Community in Giant Panda with Anorexia

The giant panda (GP) is the most precious animal in China. Gastrointestinal tract disease, especially associated with dysbiosis of gut microbiota, is the leading cause of death in GPs. Here, we performed 16S rRNA high-throughput sequencing to investigate the gut microbiota of GPs having symptoms of anorexia. Results showed that gut microbiota of GP with anorexia had lower richness (Chao1 index) than the healthy GP. However, no significant differences in alpha diversity were observed. There is a significance in the microbial structure between anorexia and healthy GPs. The abundance of phylum Firmicutes (99.23% ± 7.1%), unidentified genus Clostridiales (24.75% ± 2.5%), was significantly higher in the subadult anorexia group (P < 0.01), and that of the unidentified genus Clostridiales (4.53% ± 1.2%) was also significantly higher in the adult anorexia group (P < 0.01). Weissella and Streptococcus were found to be decreased in both anorexia groups. The decreased abundance of Weissella (0.02% ± 0.0%, 0.08% ± 0.0%) and Streptococcus (73.89% ± 4.3%, 91.15% ± 7.6%) and increase in Clostridium may cause symptoms of anorexia in giant pandas. The correlation analysis indicated that there is a symbiotic relationship among Streptococcus, Leuconostoc, Weissella, and Bacillus which are classified as probiotics (r > 0.6, P < 0.05). Importantly, a negative correlation has been found between Streptococcus and unidentified_Clostridium in two groups (r > 0.6, P < 0.05). Our results suggested that Streptococcus might be used as probiotics to control the growth of Clostridium causing the anorexia.

Immune profiles of male giant panda (Ailuropoda melanoleuca) during the breeding season

Background

The giant panda (Ailuropoda melanoleuca) is a threatened endemic Chinese species and a flagship species of national and global conservation concern. Life history theory proposes that reproduction and immunity can be mutually constraining and interrelated. Knowledge of immunity changes of male giant pandas during the breeding season is limited.

Results

Here, we researched peripheral blood gene expression profiles associated with immunity. Thirteen captive giant pandas, ranging from 9 to 11 years old, were divided into two groups based on their reproductive status. We identified 318 up-regulated DEGs and 43 down-regulated DEGs, which were enriched in 87 GO terms and 6 KEGG pathways. Additionally, we obtained 45 immune-related genes with altered expression, mostly up-regulated, and identified four hub genes HSPA4, SUGT1, SOD1, and IL1B in PPI analysis. These 45 genes were related to pattern recognition receptors, autophagy, peroxisome, proteasome, natural killer cell, antigen processing and presentation. SUGT1 and IL1B were related to pattern recognition receptors. HSP90AA1 was the most up-regulated gene and is a member of heat shock protein 90 family. HSP90 contributes to the translocation of extracellular antigen. KLRD1 encodes CD94, whose complex is an inhibitor of the cytotoxic activity of NK cells, was down-regulated. IGIP, which has the capability of inducing IgA production by B cells, was down-regulated, suggesting low concentration of IgA in male giant pandas. Our results suggest that most immune-related genes were up-regulated and more related to innate immune than adaptive immune.

Conclusions

Our results indicated that breeding male giant pandas presented an immunoenhancement in innate immunity, enhanced antigen presentation and processing in cellular immunity compared to non-breeding males. The humoral immunity of male giant pandas may show a tendency to decrease during the breeding season. This study will provide a foundation for further studies of immunity and reproduction in male giant pandas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07456-x.

Metagenomic analysis revealed a wide distribution of antibiotic resistance genes and biosynthesis of antibiotics in the gut of giant pandas

BACKGROUND

The gut microbiome is essential for the host's health and serves as an essential reservoir of antibiotic resistance genes (ARGs). We investigated the effects of different factors, including the dietary shifts and age, on the functional characteristics of the giant panda's gut microbiome (GPs) through shotgun metagenome sequencing. We explored the association between gut bacterial genera and ARGs within the gut based on network analysis.

RESULTS

Fecal samples (n=60) from captive juvenile, adult, and geriatric GPs were processed, and variations were identified in the gut microbiome according to different ages, the abundance of novel ARGs and the biosynthesis of antibiotics. Among 667 ARGs identified, nine from the top ten ARGs had a higher abundance in juveniles. For 102 ARGs against bacteria, a co-occurrence pattern revealed a positive association for predominant ARGs with Streptococcus. A comparative KEGG pathways analysis revealed an abundant biosynthesis of antibiotics among three different groups of GPs, where it was more significantly observed in the juvenile group. A co-occurrence pattern further revealed a positive association for the top ten ARGs, biosynthesis of antibiotics, and metabolic pathways.

CONCLUSION

Gut of GPs serve as a reservoir for novel ARGs and biosynthesis of antibiotics. Dietary changes and age may influence the gut microbiome's functional characteristics; however, it needs further studies to ascertain the study outcomes.

Impact of sympatric carnivores on den selection of wild giant pandas

Interspecific killing is a primary reason for the low survival rates of some animal species. The giant panda (Ailuropoda melanoleuca) is an altricial eutherian mammal and thus, in comparison to other infants, panda cubs are highly vulnerable, which may significantly influence the selection of breeding sites by females. Here, we used infrared camera traps to monitor giant panda dens for 5.5 years in Foping National Nature Reserve (FNNR) to determine how interspecific factors affect den selection by wild female pandas. Results indicated that Asian black bears (Ursusthibetanus), yellow-throated martens (Martes flavigula), leopard cats (Prionailurus bengalensis), and masked palm civets (Paguma larvata) visited the dens frequently, and the presence of these species negatively influenced den selection by female pandas. Interestingly, the presence of rodents and terrestrial birds appeared to indicate den safety, and female giant pandas were not averse and even preferred dens with a high abundance index of rodents and terrestrial birds. The den suitability index (DSI) was a reliable tool for evaluating whether dens were suitable for female giant pandas to give birth to and rear cubs, with preference for dens with high DSI values. This study increases our understanding of the den selection criteria of female giant pandas and the main threats to the survival of their cubs, thus providing important guidance for the conservation and management of this species.

A New Alignment-Free Whole Metagenome Comparison Tool and Its Application on Gut Microbiomes of Wild Giant Pandas

The comparison of metagenomes is crucial for studying the relationship between microbial communities and environmental factors. One recently published alignment-free whole metagenome comparison method based on k-mer frequencies, Libra, showed higher resolutions than the present fastest method, Mash, on whole metagenomic sequencing reads, but it did not perform as well on the assembled contigs. Here, we developed a new alignment-free tool, KmerFreqCalc, for the comparison of the whole metagenomic data, which first calculated the frequencies of both forward and reverse complementary sequences of k-mers like Mash and then computed the cosine distance between the samples based on k-mer frequency vectors like Libra. We applied KmerFreqCalc on the assembled contigs of the gut microbiomes of wild giant pandas and compared the results to Libra and Mash. The results indicated that KmerFreqCalc was able to detect the subtle difference between giant panda samples caused by seasonal diet change, showing better clustering than Libra and Mash. Therefore, KmerFreqCalc has high resolution and accuracy for assembled contigs, being very suitable for comparison of samples with low dissimilarity.

Nutrient contents predict the bamboo-leaf-based diet of Assamese macaques living in limestone forests of southwest Guangxi, China

Determining the nutrient factors influencing food choice provides important insight into the feeding strategy of animals, which is crucial for understanding their behavioral response to environmental changes. A bamboo-leaf-based diet is rare among mammals. Animals' food choice and nutritional goals have been explained by several frameworks; however, the influence of nutrients on food choice in bamboo-leaf-based macaques is not yet available. Assamese macaques (Macaca assamensis) inhabiting limestone forests are characterized by such a bamboo-leaf-based diet, predominantly consuming young leaves of Bonia saxatilis, a shrubby, karst-endemic bamboo. We studied the feeding behavior of one group of Assamese macaques using instantaneous scan sampling in limestone forests of the Guangxi Nonggang National Nature Reserve in southwest Guangxi, China. We compared the nutrient content of staple food and nonfood items and examine the role of key nutrients in the food selection of macaques. Our results showed that young leaves of bamboo B. saxatilis contained more water, crude protein, phosphorus, and less tannin than nonfood items. Furthermore, staple foods contained a higher content of water and less content of calcium than nonfood items. More specifically, quantities of water, crude protein, calcium, and phosphorus in food items were critical factors affecting feeding time on a specific plant item. Our results suggest that young bamboo leaves could meet macaques' required protein and water intake, while enabling them to maintain their mineral balance, consequently facilitating to maintain the primates' bamboo-leaf-diet in the limestone forest. Our findings confirm the effects of nutrient contents in food choice of Assamese macaques, highlighting the importance of the nutrient contents in maintaining their bamboo-based diet and the need to increase the knowledge on their nutritional strategy adapted to the bamboo-dominated diet inhabiting the unique limestone habitat.

TAS2R20 variants confer dietary adaptation to high-quercitrin bamboo leaves in Qinling giant pandas

Sensitivity to bitter tastes provides animals with an important means of interacting with their environment and thus, influences their dietary preferences. Genetic variants encoding functionally distinct receptor types contribute to variation in bitter taste sensitivity. Our previous study showed that two nonsynonymous sites, A52V and Q296H, in the TAS2R20 gene are directionally selected in giant pandas from the Qinling Mountains, which are speculated to be the causative base-pair changes of Qinling pandas for the higher preference for bamboo leaves in comparison with other pandas. Here, we used functional expression in engineered cells to identify agonists of pTAS2R20 (i.e., giant panda's TAS2R20) and interrogated the differences in perception in the in vitro responses of pTAS2R20 variants to the agonists. Our results show that pTAS2R20 is specifically activated by quercitrin and that pTAS2R20 variants exhibit differences in the sensitivity of their response to the agonist. Compared with pTAS2R20 in pandas from other areas, the receptor variant with A52V and Q296H, which is most commonly found in Qinling pandas, confers a significantly decreased sensitivity to quercitrin. We subsequently quantified the quercitrin content of the leaves of bamboo distributed in the Qinling Mountains, which was found to be significantly higher than that of the leaves of bamboo from panda habitats in other areas. Our results suggest that the decreased sensitivity to quercitrin in Qinling pandas results in higher-quercitrin-containing bamboo leaves to be tasting less bitter to them and thus, influences their dietary preference. This study illustrates the genetic adaptation of Qinling pandas to their environments and provides a fine example of the functional effects of directional selection in the giant panda.

Comparative genomics reveals bamboo feeding adaptability in the giant panda (Ailuropoda melanoleuca)

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.

Transcriptome analysis reveals immune-related gene expression changes with age in giant panda (Ailuropoda melanoleuca) blood

The giant panda (Ailuropoda melanoleuca), an endangered species endemic to western China, has long been threatened with extinction that is exacerbated by highly contagious and fatal diseases. Aging is the most well-defined risk factor for diseases and is associated with a decline in immune function leading to increased susceptibility to infection and reduced response to vaccination. Therefore, this study aimed to determine which genes and pathways show differential expression with age in blood tissues. We obtained 210 differentially expressed genes by RNA-seq, including 146 up-regulated and 64 down-regulated genes in old pandas (18-21yrs) compared to young pandas (2-6yrs). We identified ISG15, STAT1, IRF7 and DDX58 as the hub genes in the protein-protein interaction network. All of these genes were up-regulated with age and played important roles in response to pathogen invasion. Functional enrichment analysis indicated that up-regulated genes were mainly involved in innate immune response, while the down-regulated genes were mainly related to B cell activation. These may suggest that the innate immunity is relatively well preserved to compensate for the decline in the adaptive immune function. In conclusion, our findings will provide a foundation for future studies on the molecular mechanisms underlying immune changes associated with ageing.

The carnivorous digestive system and bamboo diet of giant pandas may shape their low gut bacterial diversity

The gut microbiota diversity of eight panda cubs was assessed during a dietary switch.Gut microbiota diversity of panda cubs significantly decreased after bamboo consumption.Carnivorous species living on a plant-based diet possess low microbial diversity.Mice were fed a bamboo diet but did not display low gut microbiota diversity. Giant pandas have an exclusive diet of bamboo; however, their gut microbiotas are more similar to carnivores than herbivores in terms of bacterial composition and their functional potential. This is inconsistent with observations that typical herbivores possess highly diverse gut microbiotas. It is unclear why the gut bacterial diversity of giant pandas is so low. Herein, the dynamic variations in the gut microbiota of eight giant panda cubs were measured using 16S rRNA gene paired-end sequencing during a dietary switch. Similar data from red panda (an herbivorous carnivore) and carnivorous species were compared with that of giant pandas. In addition, mice were fed a high-bamboo diet (80% bamboo and 20% rat feed) to determine whether a bamboo diet could lower the gut bacterial diversity in a non-carnivorous digestive tract. The diversity of giant panda gut microbiotas decreased significantly after switching from milk and complementary food to bamboo diet. Carnivorous species living on a plant-based diet, including giant and red pandas, possess a lower microbial diversity than other carnivore species. Mouse gut microbiota diversity significantly increased after adding high-fibre bamboo to their diet. Findings suggest that a very restricted diet (bamboo) within a carnivorous digestive system might be critical for shaping a low gut bacterial diversity in giant pandas.

Steep topography buffers threatened gymnosperm species against anthropogenic pressures in China

China is one of the most species-rich countries in the world, harboring many rare gymnosperms. Following recent human-led loss of forests, China is now experiencing increases in forest cover resulting from efforts of reforestation schemes. As anthropogenic activities have previously been found to interact with topography in shaping forest cover in China and considering the large human population and the ongoing population increase of the country, it is important to understand the role of anthropogenic pressures relative to environmental drivers for shaping species distributions here. Based on the well-established relationship between human population density and topography, we propose a hypothesis for explaining species distributions in a country dominated by human activities, predicting that species are more likely to occur in areas of steep topography under medium human population densities compared to low and high human population densities. Using species occurrence data from the Chinese Vascular Plant Distribution Database along with a common SDM method (maximum entropy modeling), we tested this hypothesis. Our results show that steep topography has the highest importance for predicting Chinese gymnosperm species occurrences in general, and threatened species specifically, in areas of medium human population densities. Consequently, these species are more often found in areas of steep terrain, supporting the proposed hypothesis. Results from this study highlight the need to include topographically heterogeneous habitats when planning new protected areas for species conservation.

Comparative skeletal anatomy of neonatal ursids and the extreme altriciality of the giant panda

Mammalian neonates are born at a wide range of maturity levels. Altricial newborns are born with limited sensory agency and require extensive parental care. In contrast, precocial neonates are relatively mature physically and often capable of independent function shortly after birth. In extant mammals, placental newborns vary from altricial to precocial, while marsupials and monotremes are all extremely altricial at birth. Bears (family Ursidae) have one of the lowest neonatal-maternal mass ratios in placental mammals, and are thought to also have the most altricial placental newborns. In particular, giant pandas (Ailuropoda melanoleuca) are thought to be exceptionally altricial at birth, and possibly marsupial-like. Here we used micro-computer (micro-computed) tomography scanning to visualize the skeletal anatomy of ursid neonates and compare their skeletal maturity with the neonates of other caniform outgroups. Specifically, we asked whether ursid neonates have exceptionally altricial skeletons at birth compared with other caniform neonates. We found that most bear neonates are similar to outgroup neonates in levels of skeletal ossification, with little variation in degree of ossification between ursine bears neonates (i.e. bears of the subfamily Ursinae). Perinatal giant pandas, however, have skeletal maturity levels most similar to a 42-45-day-old beagle fetus (~70% of total beagle gestation period). No bear exhibits the skeletal heterochronies seen in marsupial development. With regards to skeletal development, ursine bears are not exceptionally altricial relative to other caniform outgroups, but characterized largely by the drastic difference between newborn and adult body sizes. A review on the existing hypotheses for ursids' unique reproductive strategy suggests that the extremely small neonatal-maternal mass ratio of ursids may be related to the recent evolution of large adult body size, while life history characteristics retained an ancestral condition. A relatively short post-implantation gestation time may be the proximal mechanism behind the giant panda neonates' small size relative to maternal size and altricial skeletal development at birth.

Single-base-resolution methylome of giant panda's brain, liver and pancreatic tissue

The giant panda (Ailuropoda melanoleuca) is one of the most endangered mammals, and its conservation has significant ecosystem and cultural service value. Cytosine DNA methylation (5mC) is a stable epigenetic modification to the genome and has multiple functions such as gene regulation. However, DNA methylome of giant panda and its function have not been reported as of yet. Bisulfite sequencing was performed on a 4-day-old male giant panda's brain, liver and pancreatic tissues. We found that the whole genome methylation level was about 0.05% based on reads normalization and mitochondrial DNA was not methylated. Three tissues showed similar methylation tendency in the protein-coding genes of their genomes, but the brain genome had a higher count of methylated genes. We obtained 467 and 1,013 different methylation regions (DMR) genes in brain vs. pancreas and liver, while only 260 DMR genes were obtained in liver vs pancreas. Some lncRNA were also DMR genes, indicating that methylation may affect biological processes by regulating other epigenetic factors. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis indicated that low methylated promoter, high methylated promoter and DMR genes were enriched at some important and tissue-specific items and pathways, like neurogenesis, metabolism and immunity. DNA methylation may drive or maintain tissue specificity and organic functions and it could be a crucial regulating factor for the development of newborn cubs. Our study offers the first insight into giant panda's DNA methylome, laying a foundation for further exploration of the giant panda's epigenetics.

Are the gut microbial systems of giant pandas unstable?

Animals have stable dominant gut microbiomes under similar diets. Similar diets can also lead to similar gut microbial communities within host species levels. Giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens) have had long-term and stable bamboo diets, and seem well adapted to this highly fibrous diet. When compared to the gut microbiomes of Père David's deer (Elaphurus davidianus), humans, cheetah (Acinonyx jubatus), black-backed jackal (Canis-mesomelas), and black bear (Ursus thibetanus), giant panda gut microbiomes have high variation in the abundance of Pseudomonadaceae and Clostridiaceae, and are somewhat unstable. This high instability and dissimilarity may reflect an unstable gut environment, perturbation or selective pressure because of their carnivorous gastrointestinal system. A short digestive tract, brief digestion time and fast intestinal peristalsis may result in higher oxygen concentrations that select for the growth of aerobes and facultative anaerobes in giant pandas. Potential selection of high proportion of Pseudomonadaceae in giant panda (GP-HP) and red panda gut microbiomes may arise because of their postulated ability to degrade secondary compounds (e.g., cyanide compounds and aromatic compounds). However, high proportion of Clostridiaceae (GP-HF) may focus on cellulose and hemicellulose digestion. Thus, GP-HP and GP-HF groups have high dissimilarity on the functional level. These findings show that long-term similarities in diet do not always lead to similar or stable gut microbial system within the same host species and that other factors can drive the selection of gut taxa.

Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil

As Charles Darwin anticipated, living fossils provide excellent opportunities to study evolutionary questions related to extinction, competition, and adaptation. Ginkgo (Ginkgo biloba L.) is one of the oldest living plants and a fascinating example of how people have saved a species from extinction and assisted its resurgence. By resequencing 545 genomes of ginkgo trees sampled from 51 populations across the world, we identify three refugia in China and detect multiple cycles of population expansion and reduction along with glacial admixture between relict populations in the southwestern and southern refugia. We demonstrate multiple anthropogenic introductions of ginkgo from eastern China into different continents. Further analyses reveal bioclimatic variables that have affected the geographic distribution of ginkgo and the role of natural selection in ginkgo’s adaptation and resilience. These investigations provide insights into the evolutionary history of ginkgo trees and valuable genomic resources for further addressing various questions involving living fossil species.

Ginkgo is one of the living fossils from the plant kingdom. Here, authors conduct population genomics analyses to reveal its refugia and demographic history, and provide evidence of multiple anthropogenic introductions of ginkgo from eastern China into different continents.

Comparative Analysis of Gut Microbiota Among the Male, Female and Pregnant Giant Pandas (Ailuropoda Melanoleuca)

The giant panda (GP) was the most endangered species in China, and gut microbiota plays a vital role in host health. To determine the differences of the gut microbiota among the male, female and pregnant GPs, a comparative analysis of gut microbiota in GPs was carried out by 16S rRNA and ITS high-throughput sequencing. In 16S rRNA sequencing, 435 OTUs, 17 phyla and 182 genera were totally detected. Firmicutes (53.6%) was the predominant phylum followed by Proteobacteria (37.8%) and Fusobacteria (7.1%). Escherichia/Shigella (35.9%) was the most prevalent genus followed by Streptococcus (25.9%) and Clostridium (11.1%). In ITS sequencing, 920 OTUs, 6 phyla and 322 genera were also detected. Ascomycota (71.3%) was the predominant phylum followed by Basidiomycota (28.4%) and Zygomycota (0.15%). Purpureocillium (4.4%) was the most prevalent genus followed by Cladosporium (2.5%) and Pezicula (2.4%). Comparative analysis indicated that the male GPs harbor a higher abundance of phylum Firmicutes than female GPs with the contribution from genus Streptococcus. Meanwhile, the female GPs harbor a higher abundance of phylum Proteobacteria than male GPs with the contribution from genus Escherichia/ Shigella. In addition, the shift in bacteria from female to pregnant GPs indicated that phylum Firmicutes increased significantly with the contribution from Clostridium in the gut, which may provide an opportunity to study possible associations with low reproduction of the GPs.

Giant Pandas Are Macronutritional Carnivores

Giant pandas are unusual in belonging to a primarily carnivorous clade and yet being extremely specialized herbivores that feed almost exclusively on highly fibrous bamboo [1]. Paradoxically, they appear inconsistently adapted to their plant diet, bearing a mix of herbivore and carnivore traits. Herbivore traits include a skull, jaw musculature, and dentition that are adapted for fibrous diets and a specialized "pseudo-thumb" used for handling bamboo [2, 3]. They have lost functional versions of the T1R1 gene codes for umami taste receptors, which are often associated with meat eating [3]. They also have an herbivore-like subcellular distribution of the metabolic enzyme alanine: glyoxylate aminotransferase [4]. But meanwhile, giant pandas have a digestive tract [5], digestive enzymes [6], and a gut microbiota composition that resemble those of carnivores and not of herbivores [6, 7]. We draw on recent developments in multi-dimensional niche theory [8] to examine this apparent paradox. We show that the pandas' diet clustered in a macronutrient space among carnivores and was distinct from that of herbivores. The similarity with carnivore diets applied not only to the ingested diet but also to the absorbed diet, with the absorbed macronutrient ratios similar to those of the ingested foods. Comparison of the macronutrient composition of pandas' milk with those of other species shows that the carnivore-like dietary macronutrient composition extends across the life cycle. These results cast new light on the seemingly incongruous constellation of dietary adaptations in pandas, suggesting that the transition from carnivorous and omnivorous ancestry to specialized herbivory might be less abrupt than it might otherwise appear.