Reference gene catalog and metagenome-assembled genomes from the gut microbiome reveal the microbial composition, antibiotic resistome, and adaptability of a lignocellulose diet in the giant panda

The giant panda, a strict herbivore that feeds on bamboo, still retains a typical carnivorous digestive system. Reference catalogs of microbial genes and genomes are lacking, largely limiting the antibiotic resistome and functional exploration of the giant panda gut microbiome. Here, we integrated 177 fecal metagenomes of captive and wild giant pandas to construct a giant panda integrated gene catalog (GPIGC) comprised of approximately 4.5 million non-redundant genes and reconstruct 393 metagenome-assembled genomes (MAGs). Taxonomic and functional characterization of genes revealed that the captivity of the giant panda significantly changed the core microbial composition and the distribution of microbial genes. Higher abundance and prevalence of antibiotic resistance genes (ARGs) were detected in the guts of captive giant pandas, and ARG distribution was influenced by geography, for both captive and wild individuals. Escherichia, as the prevalent genus in the guts of captive giant pandas, was the main carrier of ARGs, meaning there is a high risk of ARG transmission by Escherichia. We also found that multiple mcr gene variants, conferring plasmid-mediated mobile colistin resistance, were widespread in the guts of captive and wild giant pandas. There were low proportions of carbohydrate-active enzyme (CAZyme) genes in GPIGC and MAGs compared with several omnivorous and herbivorous mammals. Many members of Clostridium MAGs were significantly enriched in the guts of adult, old and wild giant pandas. The genomes of isolates and MAGs of Clostridiaceae harbored key genes or enzymes in complete pathways for degrading lignocellulose and producing short-chain fatty acids (SCFAs), indicating the potential of these bacteria to utilize the low-nutrient bamboo diet. Overall, our data presented an exhaustive reference gene catalog and MAGs in giant panda gut and provided a comprehensive understanding of the antibiotic resistome and microbial adaptability for a high-lignocellulose diet.

Giant pandas in captivity undergo short-term adaptation in nerve-related pathways

BACKGROUND

Behaviors in captive animals, including changes in appetite, activity level, and social interaction, are often seen as adaptive responses. However, these behaviors may become progressively maladaptive, leading to stress, anxiety, depression, and other negative reactions in animals.

RESULTS

In this study, we investigated the whole-genome sequencing data of 39 giant panda individuals, including 11 in captivity and 28 in the wild. To eliminate the mountain range effect and focus on the factor of captivity only, we first performed a principal component analysis. We then enumerated the 21,474,180 combinations of wild giant pandas (11 chosen from 28) and calculated their distances from the 11 captive individuals. The 11 wild individuals with the closest distances were used for the subsequent analysis. The linkage disequilibrium (LD) patterns demonstrated that the population was almost eliminated. We identified 505 robust selected genomic regions harboring at least one SNP, and the absolute frequency difference was greater than 0.6 between the two populations. GO analysis revealed that genes in these regions were mainly involved in nerve-related pathways. Furthermore, we identified 22 GO terms for which the selection strength significantly differed between the two populations, and there were 10 nerve-related pathways among them. Genes in the differentially abundant regions were involved in nerve-related pathways, indicating that giant pandas in captivity underwent minor genomic selection. Additionally, we investigated the relationship between genetic variation and chromatin conformation structures. We found that nucleotide diversity (θπ) in the captive population was correlated with chromatin conformation structures, which included A/B compartments, topologically associated domains (TADs) and TAD-cliques. For each GO term, we then compared the expression level of genes regulated by the above four factors (AB index, TAD intactness, TAD clique and PEI) with the corresponding genomic background. The retained 10 GO terms were all coordinately regulated by the four factors, and three of them were associated with nerve-related pathways.

CONCLUSIONS

This study revealed that giant pandas in captivity undergo short-term adaptation in nerve-related pathways. Furthermore, it provides new insights into the molecular mechanism of gene expression regulation under short-term adaptation to environmental change.

The effectiveness of using giant panda as a surrogate for protecting sympatric species

The use of umbrella species to promote biodiversity conservation is practiced worldwide. The giant panda (Ailuropoda melanoleuca) an iconic species for world wildlife conservation, that inhabits regions with significant biodiversity. Given that the functions at wildlife of different trophic levels and in different body size groups are different within the ecosystem, it is unknown whether those groups of wildlife co-occurring with giant pandas are each likewise protected. To examine the umbrella effect of giant pandas on sympatric species, we used an extensive dataset of wildlife from more than 78% of giant panda habitats. We analysed the changes in distribution for four wildlife categories (large carnivores, large herbivores, medium carnivores and medium herbivores) using a generalized linear mixed model, and the underlying driving factors using binomial logistic regression models. Changes in forests in giant panda habitats were evaluated using Fragstats. The results have shown that the counts of herbivores and medium carnivores increased significantly during the decade. However, those of large carnivores significantly declined. Forest cover and nature reserves showed significant and positive effects on wildlife in 2001 and 2011, while the human population had significant and negative impacts on the herbivores and carnivores. Our results have also suggested that there has been a slight alleviation in forest fragmentation in areas unaffected by earthquakes. We concluded that the umbrella strategy of using the giant panda as an umbrella species achieved partial success by promoting the recovery of herbivores and medium carnivores. Meanwhile, this has indicated that the strategy was not sufficient for large carnivores, and therefore not enough for local ecosystems, given the critical role of large carnivores. We have suggested integrating habitat patches, controlling human disturbance, and preparing for potential human-wildlife conflict management in the Giant Panda National Park to restore large carnivore populations and maintain ecosystem functioning.

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.

Effects of environmental disinfection on microbial population and resistance genes: A case study of the microecology within a panda enclosure

Widespread use of disinfectants raises concerns over their involvement in altering microbial communities and promoting antimicrobial resistance. This study explores the influence of disinfection protocols on microbial populations and resistance genes within an isolated enclosure environment and in the gut of giant pandas (GPs) held within. Samples of panda feces, air conditioning ducts, soil and bamboo were collected before and after disinfection. High-throughput sequencing characterized the microbial flora of GP gut and environmental microbes inside the artificial habitat. Microbial cultures showed that Escherichia coli (34.6%), Enterococcus (15.4%) and other pathogenic bacteria deposited in feces and the enclosure. Isolates exhibit a consistent resistance to disinfectant, with the greatest resistance shown to cyanuric acid, and the lowest to glutaraldehyde-dodecyl dimethyl ammonium bromide (GD-DDAB) and dodecyl dimethyl ammonium bromide (DDAB). The total number of the culturable bacteria in soil and bamboo were significantly diminished after disinfection but increased in the gut. After disinfection, the richness (Chao1 index) of environment samples increased significantly (P < 0.05), while the richness in gut decreased significantly (P < 0.05). Ten genera showed significant change in feces after disinfection. Metagenome sequencing showed that 126 types of virulence genes were present in feces before disinfection and 37 in soil. After disinfection, 110 virulence genes localized in feces and 53 in soil. Eleven virulence genes including ECP and T2SS increased in feces. A total of 182 antibiotic resistance genes (ARGs) subtypes, potentially conferring resistance to 20 classes of drugs, were detected in the soils and feces, with most belonging to efflux pump protein pathways. After disinfection, the number of resistance genes increased both in gut and soil, which suggests disinfection protocols increase the number of resistance pathways. Our study shows that the use of disinfectants helps to shape the microbial community of GPs and their habitat, and increases populations of resistant strain bacteria.

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.

Identifying priority protection areas of key food resources of the giant panda

Animals that live in seasonal environments adjust their reproduction cycle to optimize seasonal forage quality. Giant pandas ( Ailuropoda melanoleuca) are seasonal altitudinal migrants that feed on bamboo shoots and leaves with different nutritional quality. However, the importance of bamboo shoots to giant pandas, especially small and isolated populations, is not fully appreciated. Here, we explored whether mating time of giant pandas is shaped by bamboo shoot phenology. We also assessed the intensity of ongoing bamboo shoot harvesting by local communities in 42 giant panda reserves based on questionnaire surveys. Varying intensity and protection levels of bamboo shoot harvesting were found. From these data, we developed a priority ranking scheme of protection areas for this key food resource. Our study showed that pandas time their mating behavior to coincide with bamboo shoot phenology due to the high nutritional demands associated with mating and pregnancy. We also found that bamboo shoots were not well protected in many places. Liangshan, Daxiangling, and Xiaoxiangling, containing the most isolated panda populations, were identified as the areas with the most urgent need of protection. Furthermore, equal attention should be paid to the QiongL-B population, as this region holds considerable potential to serve as a corridor between the Minshan and Qionglai populations. To address the challenges posed by bamboo shoot harvesting, we recommend establishing more practical bamboo shoot management policies, increasing public awareness of bamboo shoot protection, and providing alternative sources of income for local communities.

Identification of novel anelloviruses in the blood of giant panda (Ailuropoda melanoleuca)

In recent years, the continuous development of metagenomics has revealed that in addition to the digestive tract, some viruses are also common in mammalian blood. To explore and monitor potential novel viruses, in April 2015, a blood sample was collected from a healthy captive giant panda at the Chengdu Research Base of Giant Panda Breeding in Sichuan Province, China. The genomes of 25 different anelloviruses containing the complete ORF1 region have been identified. The BLASTp results showed that the amino acid sequence identity of these viruses with the best match in GenBank ranged from 27.15% to 41.29%. Based on phylogenetic analysis and SDT (Species Demarcation Tool) analysis of the complete ORF1 regions of these 25 viruses, these sequences were deduced to represent one or several novel virus genera or species. This virological study has increased our understanding of the diversity of anelloviruses in the blood of giant pandas, but further laboratory analysis is needed to verify its possible pathogenicity.

Chronic kidney disease in a giant panda (Ailuropoda melanoleuca): a case report

BACKGROUND

Chronic kidney disease (CKD) is a common cause of morbidity and mortality in captive wildlife species. However, CKD has been rarely documented in giant pandas.

CASE PRESENTATION

The following report describes a case of an eight-year-old female giant panda showing clinical signs of epistaxis, bloody diarrhea, polyuria, azotemia and anemia. The animal died despite of supportive treatments. Necropsy was performed. Grossly, both kidneys were shrunken and scarred with pallor. Subcutis edema and petechia on the epicardium of the heart were observed. The tissue samples were made into paraffin sections and stained by H.E and special staining including Periodic Acid-Schiff (PAS), von Kossa, Masson's trichrome, Phosphotungstic acid-hematoxylin (PTAH), and Congo red. Histopathology examination revealed severe chronic tubulointerstitial nephritis with marked interstitial fibrosis, glomerulosclerosis, tubular atrophy and calcification in kidneys, and acute necrotizing hemorrhagic myocarditis with calcification in heart. Other lesions included intestinal hemorrhage, hepatic fatty degeneration and necrosis with hemosiderin, and splenic hemosiderin.

CONCLUSIONS

In summary, chronic kidney disease was finally diagnosed based on the association of clinical, gross, and histopathological findings. Heart failure secondary to CKD is the leading cause of death in this giant panda. The potential cause of CKD in this animal is possibly due to long term and uncontrolled hypertension. Blood pressure monitoring is essential in establishing the diagnosis and management of hypertension in giant panda.

Characteristics of mesenchymal stem cells and their exosomes derived from giant panda (Ailuropoda melanoleuca) endometrium

Conservation of genetic resources is an important way to protect endangered species. At present, mesenchymal stem cells (MSCs) have been isolated from the bone marrow and umbilical cords of giant pandas. However, the types and quantities of preserved cell resources were rare and limited, and none of MSCs was derived from female reproductive organs. Here, we first isolated MSCs from the endometrium of giant panda. These cells showed fibroblast morphology and expressed Sox2, Klf4, Thy1, CD73, CD105, CD44, CD49f, and CD105. Endometrium mesenchymal stem cells (eMSCs) of giant panda could induce differentiation into three germ layers in vitro. RNA-seq analysis showed that 833 genes were upregulated and 716 genes were downregulated in eMSCs compared with skin fibroblast cells. The results of GO and the KEGG analysis of differentially expressed genes (DEGs) were mainly focused on transporter activity, signal transducer activity, pathways regulating pluripotency of stem cells, MAPK signaling pathway, and PI3K-Akt signaling pathway. The genes PLCG2, FRK, JAK3, LYN, PIK3CB, JAK2, CBLB, and MET were identified as hub genes by PPI network analysis. In addition, the exosomes of eMSCs were also isolated and identified. The average diameter of exosomes was 74.26 ± 13.75 nm and highly expressed TSG101 and CD9 but did not express CALNEXIN. A total of 277 miRNAs were detected in the exosomes; the highest expression of miRNA was the has-miR-21-5p. A total of 14461 target genes of the whole miRNAs were predicted and proceeded with functional analysis. In conclusion, we successfully isolated and characterized the giant panda eMSCs and their exosomes, and analyzed their functions through bioinformatics techniques. It not only enriched the conservation types of giant panda cell resources and promoted the protection of genetic diversity, but also laid a foundation for the application of eMSCs and exosomes in the disease treatment of giant pandas.

Evaluation of the immunoprotective effects of eight recombinant proteins from Baylisascaris schroederi in mice model

BACKGROUND

Baylisascaris schroederi is the most common and harmful intestinal parasitic nematode of giant pandas, causing ascariasis. Although drug deworming is the main measure to control ascariasis in captive giant pandas, prolonged and repeated use of deworming drugs might induce resistance in nematodes and drug residues in giant pandas. Therefore, developing a safe and effective vaccine might provide a novel strategy to prevent ascariasis in captive giant pandas.

METHODS

Four highly expressed secretome genes encoding excretory and secretory proteins of B. schroederi, including transthyretin-like protein 46 (BsTLP), uncharacterized protein (BsUP), hypothetical protein 1 (BsHP1), and hypothetical protein 2 (BsHP2) and four functional genes [(encoding Galectin (BsGAL), glutathione S-transferase (BsGST), fatty acid-binding protein (BsFABP), and thioredoxin peroxidase (BsTPX)] were identified based on genome and transcriptome databases of B. schroederi and used to construct recombinant proteins via prokaryotic expression. Kunming mice were vaccinated subcutaneously twice with the recombinant proteins (50 μg/mouse) mixed with Quil A adjuvant with a 2-week interval and then orally challenged with 3000 infective eggs. The immunoprotective effects of the eight recombinant proteins on mice were assessed comprehensively using surface lesion histology scores of the mouse liver and lung, larval worm reduction, serum antibody levels (IgG, IgE, IgA, IgG1, and IgG2a), and cytokine production [interferon gamma (IFN-γ), interleukin (IL)-2, IL-4, IL-5, and IL-10].

RESULTS

Mice vaccinated with recombinant (r)BsUP (76.5%), rBsGAL (74.7%), and rBsHP2 (71.5%) showed a significant (P < 0.001) reduction in the larval worm rate compared with that in the adjuvant control. Besides, the surface lesions in the liver and lung of the vaccinated mice were alleviated. Serum levels of total IgG, IgE, IgA, IgG1, IgG2a, and cytokines, including IL-10, IL-5, and IFN-γ, were significantly higher (P < 0.001) than those in the control group.

CONCLUSIONS

The results showed that candidate three vaccines (rBsUP, rBsGAL, and rBsHP2) could provide effective protection against egg infection in mice associated with a mixed Th1/2-type immune response.

Assessing the effectiveness of protected areas for panda conservation under future climate and land use change scenarios

While the relatively stable land use and land cover (LULC) patterns is an important feature of protected areas (PAs), the influence of this feature on future species distribution and the effectiveness of the PAs has rarely been explored. Here, we assessed the role of land use patterns within PAs on the projected range of the giant panda (Ailuropoda melanoleuca) by comparing projections inside and outside of PAs for four model configurations: (1) only climate covariates, (2) climate and dynamic land use covariates, (3) climate and static land use covariates and (4) climate and hybrid dynamic-static land use covariates. Our objectives were twofold: to understand the role of protected status on projected panda habitat suitability and evaluate the relative efficacy of different climate modeling approaches. The climate and land use change scenarios used in the models include two shared socio-economic pathways (SSPs) scenarios: SSP126 [an optimistic scenario] and SSP585 [a pessimistic scenario]. We found that models including land-use covariates performed significantly better than climate-only models and that these projected more suitable habitat than climate-only models. Static land-use models projected more suitable habitat than both the dynamic and hybrid models under SSP126, while these models did not differ under SSP585. China's panda reserve system was projected to effectively maintain suitable habitat inside PAs. Panda dispersal ability also significantly impacted outcomes, with most models assuming unlimited dispersal forecasting range expansion and models assuming zero dispersal consistently forecasting range contraction. Our findings highlight that policies targeting improved land-use practices should be an effective means for offsetting some of the negative effects of climate change on pandas. As the effectiveness of PAs is projected to be maintained, we recommend the judicious management and expansion of the PA system to ensure the resilience of panda populations into the future.

Diagnosis and surgical management of testicular seminoma in captive giant panda (Ailuropoda melanoleuca)

A captive adult male giant panda (Guangzhou Chimelong Safari Park, CHINA) presented with azoospermia and enlarged left testicle. A tentative diagnosis of testicular neoplasia was made, and confirmed as testicular seminoma cases by testicular ultrasound, computed tomography (CT), testicular biopsy, and tumor marker examination findings. Based on the diagnostic results, the treatment of choice was surgical resection of the testicular tumor under general anesthesia. And the histopathological findings of the excised neoplasm were consistent with those of testicular seminoma. In addition, no tumor recurrence was observed after surgery, which indicates that our surgical and post-operative treatments were effective. The surgical treatment adopted in this case report is safe for patients and provides the best solution for the diagnosis and treatment of giant panda testicular seminoma. To our knowledge this is the first detailed report of surgical resection of testicular seminoma in the giant panda.

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.

Stable seasonal migration patterns in giant pandas

A critical function of animal movement is to maximize access to essential resources in temporally fluctuating and spatially heterogeneous environments. Seasonally mediated resource fluctuations may influence animal movements, enabling them to track changing resource distributions, resulting in annual migration patterns. The conservation-dependent giant panda ( Ailuropoda melanoleuca) displays seasonal movement patterns; however, the key factor driving these seasonal migration patterns remains poorly understood. Here, we used GPS tracking collars to monitor the movements of six giant pandas over a 12-year period across different elevations, and performed statistical analysis of seasonal migration directions, routes, habitat revisitation, home range overlap, first arrival events, and stability. Our results revealed a compelling pattern of seasonal migrations that facilitated the ability of the pandas to forage at the appropriate time and place to maximize nutritional intake. Our results indicated that pandas utilize spatial memory to locate reliable food resources, as evidenced by their annual return to the same or similar winter and summer home ranges and the consistently maintained percentage of home range overlap. These novel insights into giant panda foraging and movement ecology not only enhance our understanding of its ability to adapt to nutritionally poor dietary resources but also provide important information for the development of resource utilization-based protection and management strategies.

Complete genome sequence of a blaNDM-5-producing Escherichia coli DC71 assigned as ST410-O8:H9 and recovered from a captive giant panda (Ailuropoda melanoleuca) in China

OBJECTIVES

In this study, we report the complete genome sequence of a multidrug-resistant Escherichia coli strain recovered from a fecal sample from a captive giant panda in China.

METHODS

Antimicrobial susceptibility testing was performed. Genomic DNA from E. coli DC71 was sequenced using a Nanopore PromethION sequencer instrument (Oxford Nanopore Technologies, UK) and MGI High-throughput Sequencing MGISEQ-2000 platforms. The clean reads were de novo assembled using SPAdes v3.11. The complete genome was annotated and analyzed using multilocus sequence typing, serotyping, plasmid replicons, fimH typing, chromosomal point mutations, acquired antimicrobial resistance, and virulence genes with web tools available at the Center for Genomic Epidemiology.

RESULTS

The complete genome, 4 991 906 bp in length and comprising 4677 protein-coding sequences, was generated. In silico analysis revealed that E. coli DC71 belonged to the ST410-O8:H9 subclone. A carbapenem resistance gene, bla_NDM-5, was located on the pDC71-2 plasmid, coproducing bla_TEM-1. Many other resistance determinants encoded by chromosomes and pDC71-3 were found. The virulence related genes carried by chromosomes were mostly related to enterohemorrhagic E. coli (EHEC) O157:H7.

CONCLUSIONS

To our knowledge, this is the first complete genome of an E. coli ST410-O8:H9 strain recovered from captive giant panda in China. This multidrug-resistant E. coli subclone may pose potential risks to human and animal health. The genome sequence will be helpful to understand the genomic structure, its diversity, and the molecular mechanism allowing bacteria to disseminate the resistance gene.

Virulence gene detection and antimicrobial resistance analysis of Enterococcus faecium in captive giant pandas (Ailuropoda melanoleuca) in China

BACKGROUND

The emergence of multidrug resistance among enterococci makes effective treatment of enterococcal infections more challenging. Giant pandas (Ailuropoda melanoleuca) are vulnerable to oral trauma and lesions as they feast on bamboo. Enterococci may contaminate such oral lesions and cause infection necessitating treatment with antibiotics. However, few studies have focused on the virulence and drug resistance of oral-derived enterococci, including Enterococcus faecium, in giant pandas. In this study, we analyzed the prevalence of 8 virulence genes and 14 drug resistance genes in E. faecium isolates isolated from saliva samples of giant pandas held in captivity in China and examined the antimicrobial drug susceptibility patterns of the E. faecium isolates.

RESULTS

Twenty-eight isolates of E. faecium were successfully isolated from the saliva samples. Four virulence genes were detected, with the acm gene showing the highest prevalence (89%). The cylA, cpd, esp, and hyl genes were not detected. The isolated E. faecium isolates possessed strong resistance to a variety of drugs; however, they were sensitive to high concentrations of aminoglycosides. The resistance rates to vancomycin, linezolid, and nitrofurantoin were higher than those previously revealed by similar studies in China and other countries.

CONCLUSIONS

The findings of the present study indicate the drugs of choice for treatment of oral E. faecium infection in the giant panda.

Gene expressions between obligate bamboo-eating pandas and non-herbivorous mammals reveal converged specialized bamboo diet adaptation

BACKGROUND

It is inevitable to change the function or expression of genes during the environmental adaption of species. Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to Carnivora and have developed similar adaptations to the same dietary switch to bamboos at the morphological and genomic levels. However, the genetic adaptation at the gene expression level is unclear. Therefore, we aimed to examine the gene expression patterns of giant and red panda convergent specialized bamboo-diets. We examined differences in liver and pancreas transcriptomes between the two panda species and other non-herbivorous species.

RESULTS

The clustering and PCA plots suggested that the specialized bamboo diet may drive similar expression shifts in these two species of pandas. Therefore, we focused on shared liver and pancreas DEGs (differentially expressed genes) in the giant and red panda relative to other non-herbivorous species. Genetic convergence occurred at multiple levels spanning carbohydrate metabolism, lipid metabolism, and lysine degradation. The shared adaptive convergence DEGs in both organs probably be an evolutionary response to the high carbohydrate, low lipid and lysine bamboo diet. Convergent expression of those nutrient metabolism-related genes in both pandas was an intricate process and subjected to multi-level regulation, including DNA methylation and transcription factor. A large number of lysine degradation and lipid metabolism related genes were hypermethylated in promoter regions in the red panda. Most genes related to carbohydrate metabolism had reduced DNA methylation with increased mRNA expression in giant pandas. Unlike the red panda, the core gene of the lysine degradation pathway (AASS) doesn't exhibit hypermethylation modification in the giant panda, and dual-luciferase reporter assay showed that transcription factor, NR3C1, functions as a transcriptional activator in AASS transcription through the binding to AASS promoter region.

CONCLUSIONS

Our results revealed the adaptive expressions and regulations of the metabolism-related genes responding to the unique nutrients in bamboo food and provided data accumulation and research hints for the future revelation of complex mechanism of two pandas underlying convergent adaptation to a specialized bamboo diet.

Community structure of the solitary giant pandas is maintained by indirect social connections

BACKGROUND

Indirect interactions between individual solitary mammals, such as the giant panda, are often overlooked because of their nature, yet are important for maintaining the necessary sociality in solitary species.

METHODS AND RESULTS

Here, we determined the genetic identity of all giant panda individuals in a local population and matched these identities with their associations to determine social network of this solitary animal. Total thirty-five giant panda individuals were found in our field survey, and we constructed genetic and social networks for thirty-three individuals who successfully obtained genetic, age and sex information. The results showed that sex had great impact on both social network and genetic network, and age may have the potential to influence the social network of the giant pandas. Adult males, mostly in the central of the social network, which appeared significantly larger network connections than adult females. Due to the female-biased dispersal pattern of wild giant pandas, male-male pairs showed higher relatedness than female-female ones and multi-generational patrilinear assemblages are expected in the study area.

CONCLUSIONS

The relatedness of individuals has an influence on the formation of community social structure of giant pandas, and indirect interactions among solitary giant pandas potentially function to reduce competition for resources and inbreeding.

An improved, chromosome-level genome of the giant panda (Ailuropoda melanoleuca)

BACKGROUND

The iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species.

FINDINGS

We generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named "Jing Jing" with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for "Jing Jing", which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI).

CONCLUSIONS

We presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.

Identification of extended-spectrum beta-lactamase (CTX-M)-producing Klebsiella pneumoniae belonging to ST37, ST290, and ST2640 in captive giant pandas

Background

Extended-spectrum β-lactamases (ESBL)-producing strains of Klebsiella pneumoniae remain a worldwide, critical clinical concern. However, limited information was available concerning ESBL-producing Klebsiella pneumoniae in giant pandas. The objective of this study was to characterize ESBL-producing Klebsiella pneumoniae isolates from captive giant pandas. A total of 211 Klebsiella pneumoniae isolates were collected from 108 giant pandas housed at the Chengdu Research Base of Giant Panda Breeding (CRBGP), China. Samples were screened for the ESBL-producing phenotype via the double-disk synergy test.

Result

A total of three (1.42%, n = 3/211) ESBL-producing Klebsiella pneumoniae strains were identified, and characterization of ESBL-producing Klebsiella pneumoniae isolates were studied by the detection of ESBL genes and mobile genetic elements (MGEs), evaluation of antimicrobial susceptibility and detection of associated resistance genes. Clonal analysis was performed by multi-locus sequencing type (MLST). Among the three ESBL-producing isolates, different ESBL-encoding genes, including bla_CTX-M, and bla_TEM, were detected. These three isolates were found to carry MGEs genes (i.e., IS903 and tnpU) and antimicrobial resistance genes (i.e., aac(6')-Ib, aac(6')-I, qnrA, and qnrB). Furthermore, it was found that the three isolates were not hypermucoviscosity, resistant to at least 13 antibiotics and belonged to different ST types (ST37, ST290, and ST2640).

Conclusion

Effective surveillance and strict infection control strategies should be implemented to prevent outbreaks of ESBL-producing Klebsiella pneumoniae in giant pandas.

Viral metagenomics unveiled extensive communications of viruses within giant pandas and their associated organisms in the same ecosystem

Cross-species transmission events were commonplace, with numerous cases of host-switching during the viral evolutionary history, but relatively little evidence for onward transmission in different species living in the same ecosystem. For understanding the communications of viruses in giant pandas (Ailuropoda melanoleuca) and their associated organisms, based on a large size of samples (N = 2305) collected between 2015 and 2020 from giant panda (N = 776) and other four giant panda-associated organisms in the same ecosystem, red pandas (N = 700), stray cats (N = 32), wild rats (N = 42), and mosquitoes (N = 755), viromics was used for the virus identification and subsequent virus traceability. The results showed that a feline panleukopenia virus (FPV) was found in giant pandas with clinical signs of vomiting and mild diarrhea. Meanwhile, the same FPV strain was also prevalent in the healthy red panda (Ailurus fulgens) population. From the viromes of the five different organisms, 250 virus genomes were determined. Our data revealed that besides FPV, other putative pathogenic viruses, such as red panda amdoparvoviruses (RPAVs) and Getah viruses (GETVs) were responsible for previous disease or death of some red pandas. We also demonstrated that a number of viruses were involved in potential interspecies jumping events between giant pandas and their associated species. Collectively, our results shed light on the genetic diversity and relationship of diverse viral pathogens in 'Giant pandas-Associated animals-Arthropods' and report some cases of possible viral host-switching among these host species living in the same ecosystem.

Comparison of genetic characteristics between captive and wild giant pandas based on 13 mitochondrial coding genes

BACKGROUND

Research on genetic diversity based on mitochondrial DNA of giant pandas mainly focused on a single marker or a few genes.

OBJECTIVE

To provide a more comprehensive assessment of the genetic diversity on giant pandas based on 13 mitochondrial protein coding genes.

METHODS

We assembled 13 protein coding genes in the mitochondrial genome of the giant panda based on the whole genome sequencing data, including ND1, ND2, COX1, COX2, ATP8, ATP6, COX3, ND3, ND4L, ND4, ND5, ND6 and Cyt b.

RESULTS

We successfully obtained long sequence of 11,416 base pairs with all 13 genes for 110 giant panda individual, accounting for 67.93% in length of the mitochondrial reference genome. Haplotype diversity was 0.9518 ± 0.009 and nucleotide diversity (π) was 0.00157 ± 0.00014. We detected three new haplotypes, including GPC10 and GPC21 for the CR sequence and GPB12 for the Cyt b gene.

CONCLUSION

These multi-gene sequences provided more genetic variable information to compare captive and wild giant panda population.

Epigenomic profiling indicates a role for DNA methylation in the postnatal liver and pancreas development of giant pandas

Giant pandas are unique within Carnivora with a strict bamboo diet. Here, the epigenomic profiles of giant panda liver and pancreas tissues collected from three important feeding stages were investigated using BS-seq. Few differences in DNA methylation profiles were exhibited between no feeding and suckling groups in both tissues. However, we observed a tendency toward a global loss of DNA methylation in the gene-body and promoter region of metabolism-related genes from newborn to adult. Correlation analysis revealed a significant negative correlation between the changes in methylation levels within gene promoters and gene expression. The majority of genes related to nutrition metabolism had lost DNA methylation with increased mRNA expression in adult giant pandas. The few galactose metabolism and unsaturated fatty acid metabolism related genes that were hypomethylated and highly-expressed at early stages of giant panda development may meet the nutritional requirement of this species' highly altricial neonates.

Whole blood transcriptome profiling identifies candidate genes associated with alopecia in male giant pandas (Ailuropoda melanoleuca)

Background

The giant panda (Ailuropoda melanoleuca) is a threatened species endemic to China. Alopecia, characterized by thinning and broken hair, mostly occurs in breeding males. Alopecia significantly affects the health and public image of the giant panda and the cause of alopecia is unclear.

Results

Here, we researched gene expression profiles of four alopecia giant pandas and seven healthy giant pandas. All pandas were approximately ten years old and their blood samples collected during the breeding season. A total of 458 up-regulated DEGs and 211 down-regulated DEGs were identified. KEGG pathway enrichment identified that upregulated genes were enriched in the Notch signaling pathway and downregulated genes were enriched in ribosome, oxidative phosphorylation, and thermogenesis pathways. We obtained 28 hair growth-related DEGs, and identified three hub genes NOTCH1, SMAD3, and TGFB1 in PPI analysis. Five hair growth-related signaling pathways were identified with abnormal expression, these were Notch, Wnt, TGF-β, Mapk, and PI3K-Akt. The overexpression of NOTCH1 delays inner root sheath differentiation and results in hair shaft abnormalities. The delayed hair regression was associated with a significant decrease in the expression levels of TGFB1.

Conclusions

Our data confirmed the abnormal expression of several hair-related genes and pathways and identified alopecia candidate genes in the giant panda. Results of this study provide theoretical basis for the establishment of prevention and treatment strategies for giant pandas with alopecia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08501-z.

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.

Multi-functional Potential of Five Lactic Acid Bacteria Strains Derived from Giant Panda (Ailuropoda melanoleuca)

The multi-functional properties of lactic acid bacteria (LAB) on host health have been a popular research topic. The aim of present study was to assess the multi-functional potential of five LAB strains isolated from giant panda. In this study, we analyzed five giant panda LAB strains (Weissella confuse WJ202003 (W3), WJ202009 (W9), WJ202021 (W21), BSP201703 (X₃); Lactiplantibacillus plantarum BSGP201683 (G83)) and found that they exhibited rapid growth as well as strong acid production capacity. The five LAB strains possessed high cell surface hydrophobicity to the four tested solvents (xylene, hexadecane, chloroform, ethyl acetate; except strain W9), auto-aggregation ability, co-aggregation ability with three pathogens (Escherichia coli, Enterotoxigenic Escherichia coli, Salmonella), adhesion ability to Caco-2 cell line, and strongly biofilm formation ability, suggesting an adhesion property. As investigated for their antioxidative potential, all the strains showed good tolerance to H₂O₂, high scavenging ability against 1, 1-diphenyl-2-picrylhydrazyl (DPPH), and hydroxyl (OH^-), and reduction ability. Furthermore, the five LAB strains could produce multiple probiotic substances, including exopolysaccharide (EPS), gamma-aminobutyric acid (GABA), bile salt hydrolase (BSH), cellulase (only strain G83), and protease (except strain X₃), which was the first to report the production of EPS, GABA, BSH, cellulase, and protease in giant panda-derived LAB strain. These results demonstrated that strains W3, W9, W21, X₃, and G83 had multi-functional potential and could be utilized as potential probiotics for giant panda.

Blood transcriptome analysis revealed the immune changes and immunological adaptation of wildness training giant pandas

The giant panda (Ailuropoda melanoleuca) is a global flagship species for biodiversity conservation. As the time for captive giant pandas to be released into the wild matures, wildness training is provided to allow adaptation to their natural environment. It is assumed that changes in the immune system would be integral in this adaptation from captive to wild, where many more pathogens would be encountered in their natural habitats. Therefore, this study aims to determine the expression changes of immune-related genes and their potential as immunoassay markers for adaptation monitoring in wildness training giant pandas, and then to understand the adaptation strategy of wildness training giant pandas to the wild environment, thereby improving the success rate of panda reintroduction. We obtained 300 differentially expressed genes (DEGs) by RNA-seq, with 239 up-regulated and 61 down-regulated DEGs in wildness training giant pandas compared to captive pandas. Functional enrichment analysis indicated that up-regulated DEGs were enriched in several immune-related terms and pathways. There were 21 immune-related DEGs, in which most of them were up-regulated in wildness training giant pandas, including several critical innate and cellular immune genes. IL1R2 was the most significantly up-regulated gene and is a signature of homeostasis within the immune system. In the protein-protein interaction (PPI) analysis, CXCL8, CXCL10, and CCL5 were identified as the hub immune genes. Our results suggested that wildness training giant pandas have stronger innate and cellular immunity than captive giant pandas, and we proposed that a gene set of CXCL8, CXCL10, CCL5, CD3D, NFKBIA, TBX21, IL12RB2, and IL1R2 may serve as potential immunoassay markers to monitor and assess the immune status of wildness training giant pandas. Our study offers the first insight into immune alterations of wildness training giant pandas, paving the way for monitoring and evaluating the immune status of giant pandas when reintroducing them into the wild.

The potential risks and exposure of Qinling giant pandas to polycyclic aromatic hydrocarbon (PAH) pollution

Rapid industrialization and urbanization have created a substantial urban-rural gradient for various pollutants. The Qinling Mountains are highly important in terms of biodiversity, providing habitat for giant pandas, which are endemic to China and are a widely recognized symbol for conservation. Whether polycyclic aromatic hydrocarbon (PAH) exposure risks regarding in situ animal conservation zones are affected by environmental pollution or even enhanced by the mountain-trapping effect requires further research. Our group carried out a large-scale investigation on the area ranging from Xi'an to Hanzhong across the giant panda habitat in the Qinling Mountains by collecting atmosphere, soil, bamboo, and fecal samples from different sites over a two-year period. The total toxicity of atmospheric PAHs and the frequencies of soil PAHs above effect range low (ERL) values showed a decreasing trend from urban areas to the central mountains, suggesting a distance effect from the city. The proportions of total 5- and 6-ring PAHs in the atmosphere were higher in the central mountainous areas than in the urban areas, while this difference was reversed in the soil. Health risk assessments showed that the incremental lifetime carcinogenic risks (ILCR) of PAH exposure by bamboo ingestion ranged from 2.16 × 10^-4 to 3.11 × 10^-4, above the critical level of 10^-4. Bamboo ingestion was the main driver of the PAH exposure risks. The concentration difference between bamboo and fecal samples provided a reference for the level of PAHs absorbed by the panda digestive system. Since the Qinling Mountains possess the highest density of giant pandas and provide habitats to many other endangered animal and plant species, we should not ignore the probability of health risks posed by PAHs. Monitoring the pollution level and reducing the atmospheric emissions of toxic pollutants are recommended actions. Further detailed research should also be implemented on pandas' health effects of contaminant exposure.

Nanoparticles of conformation-stabilized canine distemper virus hemagglutinin are highly immunogenic and induce robust immunity

Background

Canine distemper virus (CDV) infection of ferrets, dogs, and giant pandas causes an acute systemic disease involving multiple organ systems, including the respiratory tract, lymphoid system, and central nervous system. In this study, we tested a new candidate CDV vaccine-CDV nanoparticles-based on hemagglutinin protein.

Methods

The nanoparticles were generated from conformation-stabilized CDV hemagglutinin tetramers. Immune responses against CDV were evaluated in mice. Immunization was initiated 6 weeks after birth and boosted two times with 4-week intervals. The blood and mucosal samples were collected 2 weeks after each immunization.

Results

Vaccination with CDV nanoparticles elicited high levels of IgG antibody titers in mice (approximately sevenfold to eightfold higher than that obtained with soluble CDV H protein) and mucosal immune responses and developed increased CDV-specific neutralizing antibody. The mice that received nanoparticles showed significantly higher IFN-γ- and IL-4-secreting cell population in the spleen and lymph node compared with mice immunized with soluble H protein. The co-stimulatory molecular expression of CD80 and CD86 on the surface of DCs was also upregulated.

Conclusion

The results demonstrate that self-assembly into nanoparticles can increase the immunogenicity of vaccine antigens, and nanoparticles assembled from conformation-stabilized CDV H protein can serve as a new CDV vaccine.

Chromosome-scale genomes provide new insights into subspecies divergence and evolutionary characteristics of the giant panda

Extant giant pandas are divided into Sichuan and Qinling subspecies. The giant panda has many species-specific characteristics, including comparatively small organs for body size, small genitalia of male individuals, and low reproduction. Here, we report the most contiguous, high-quality chromosome-level genomes of two extant giant panda subspecies to date, with the first genome assembly of the Qinling subspecies. Compared with the previously assembled giant panda genomes based on short reads, our two assembled genomes increased contiguity over 200-fold at the contig level. Additional sequencing of 25 individuals dated the divergence of the Sichuan and Qinling subspecies into two distinct clusters from 10,000 to 12,000 years ago. Comparative genomic analyses identified the loss of regulatory elements in the dachshund family transcription factor 2 (DACH2) gene and specific changes in the synaptotagmin 6 (SYT6) gene, which may be responsible for the reduced fertility of the giant panda. Positive selection analysis between the two subspecies indicated that the reproduction-associated IQ motif containing D (IQCD) gene may at least partly explain the different reproduction rates of the two subspecies. Furthermore, several genes in the Hippo pathway exhibited signs of rapid evolution with giant panda-specific variants and divergent regulatory elements, which may contribute to the reduced inner organ sizes of the giant panda.

Identification of a novel circovirus in blood sample of giant pandas (Ailuropoda melanoleuca)

The members of the family Circoviridae are considered to be one of the smallest autonomously replicating viruses that are classified into two genera, Circovirus and Cyclovirus. Circoviruses have been found in a variety of vertebrates, but whether they infect endangered protected animals has not been studied in much detail. Here, viral metagenomics and PCR methods were used to detect and verify viral nucleic acid in the blood sample from giant pandas. According to these methods, the complete genome sequence of a novel circovirus, the giant panda associated circovirus (GPCV) from the blood sample of three giant pandas was identified. The GPCV genome is 2090 bp in size and reveals two putative ambisense open-reading frames, encoding the major structural capsid protein and the replication associated protein, respectively, the latter having two predicted introns. Pairwise sequence comparison and phylogenetic analyses indicated GPCV was a putative new species within genus Circovirus based on the species demarcation criteria of the International Committee on the Taxonomy of Viruses. It is the first time that circovirus has been identified from blood sample of giant pandas. These efforts will contribute to future analyses to illuminate the evolutionary relationships between classified and newly identified members of the family Circoviridae.

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.

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.

First report of fatal baylisascariasis-induced acute pancreatitis in a giant panda

A wild adult male giant panda that was rescued from a nature reserve in Sichuan Province, China, has died. The panda had been in poor physical condition: it was wheezing and had increased serum amylase. A pathological examination was performed in order to determine the cause of death. Gross examination revealed 1380 mL of yellowish fluid in the abdominal cavity, 356 nematodes in the digestive tract and one filling the pancreatic duct, contractions and variably-sized dark purple areas in the spleen, a collapsed right lung and consolidation of the left lung. Acute pancreatitis was confirmed histopathologically via edema, focal necrosis and hemorrhage with inflammatory cell infiltration. Other major histopathological changes included serous-hemorrhagic pneumonia, lymphocytic necrosis and depletion in the spleen, and degeneration and necrosis of renal tubular epithelial cells. The nematodes were identified as Baylisascaris schroederi via molecular assays. In conclusion, the cause of death of the giant panda was determined to be multiple organ dysfunction syndrome caused by baylisascariasis-induced acute pancreatitis. To our knowledge, this is the first report of fatal baylisascariasis-induced acute pancreatitis in the giant panda.

Development of recombinase polymerase amplification assays for rapid and visual detection of canine distemper virus infecting giant panda

Background

Canine distemper virus (CDV) is an enveloped negative-strand RNA virus that exhibits a high mutation rate and continuously expands the range of hosts. Notably, CDV has infected giant panda with spill over from viral reservoirs in canines. Giant pandas (Ailuropoda melanoleuca), especially captive pandas, are known to be susceptible to natural infection with CDV. The high fatality rate of CDV poses a serious threat to the safety of the giant panda population. However, vaccines or drugs for canine distemper in giant pandas have not been developed to date. Therefore, a rapid test that can achieve accurate onsite detection of CDV is important to enable the timely implementation of control measures. In this study, we established a nucleic acid visualization assay for targeting the CDV N gene by using combines reverse transcription recombinase polymerase amplification with a closed vertical flow visualization strip (RT-RPA-VF).

Results

The RT-RPA-VF assay does not require sophisticated equipment, and it was determined to provide rapid detection at 35 °C for 30 min, while the limit of detection was 5 × 10¹ copies/μl RNA transcripts and 10^0.5 TCID₅₀ ml^− 1 viruses. The results showed that the assay was high specific to CDV and had no cross-reactivity with other viruses infecting the giant panda. Compared with RT-qPCR, RT-RPA-VF assay had a sensitivity of 100% and a specificity of 100% in 29 clinical samples. The coincidence rate between RT-RPA-VF and RT-qPCR was 100% (kappa = 1), indicating that the RT-RPA-VF assay possessed good diagnostic performance on clinical samples.

Conclusions

The RT-RPA-VF provides a novel alternative for the simple, sensitive, and specific identification of CDV and showed great potential for point of care diagnostics for captive and wild giant panda.

A tail of two pandas- whole genome k-mer signature analysis of the red panda (Ailurus fulgens) and the Giant panda (Ailuropoda melanoleuca)

Background

The red panda (Ailurus fulgens) is a riddle of morphology, making it hard to tell whether it is an ursid, a procyonid, a mustelid, or a member of its own family. Previous genetic studies have given quite contradictory results as to its phylogenetic placement.

Results

A recently developed whole genome-based algorithm, the Whole Genome K-mer Signature algorithm was used to analyze the genomes of 28 species of Carnivora, including A. fulgens and several felid, ursid, mustelid, one mephitid species. This algorithm has the advantage of holistically using all the information in the genomes of these species. Being a genomics-based algorithm, it also reduces stochastic error to a minimum. Besides the whole genome, the mitochondrial DNA from 52 mustelids, mephitids, ursids, procyonids and A. fulgens were aligned to draw further phylogenetic inferences.

The results from the whole genome study suggested that A. fulgens is a member of the mustelid clade (p = 9·10^− 97). A. fulgens also separates from the mephitid Spilogala gracilis. The giant panda, Ailuropoda melanoleuca also clusters away from A. fulgens, together with other ursids (p = 1.2·10^− 62). This could be due to the geographic isolation of A. fulgens from other mustelid species. However, results from the mitochondrial study as well as neighbor-joining methods based on the sequence identity matrix suggests that A. fulgens forms a monophyletic group. A Maximum Likelihood tree suggests that A. fulgens and Ursidae form a monophyletic group, although the bootstrap value is weak.

Conclusions

The main conclusion that we can draw from this study is that on a whole genome level A. fulgens possibly belongs to the mustelid clade, and not an ursid or a mephitid. This despite the fact that previously some researchers classified A. fulgens and A. melanoleuca as relatives. Since the genotype determines the phenotype, molecular-based classification takes precedence over morphological classifications. This affirms the results of some previous studies, which studied smaller portions of the genome. However, mitochondrial analyses based on neighbor-joining and maximum likelihood methods suggest otherwise.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07531-3.

Isolation and characterization of mesenchymal stem cells from umbilical cord of giant panda

Umbilical cord-derived mesenchymal stem cells (UC-MSCs) constitute a class of cells with significant self-renewal and multilineage differentiation properties and have great potential for therapeutic applications and the genetic conservation of endangered animals. In this study, we successfully isolated and cultured UC-MSCs from the blood vessels of giant panda umbilical cord (UC). The cells were arranged in a vortex or cluster pattern and exhibited a normal karyotype, showing the morphological characteristics of fibroblasts. In addition, we found that basic fibroblast growth factor (bFGF) and epidermal growth factor (EGF) promoted cell proliferation, whereas stem cell factor (SCF) did not promote cell proliferation. Cultured UC-MSCs were negative for CD34 (hematopoietic stem cell marker) and CD31 (endothelial cell marker), but positive for MSC markers (CD44, CD49f, CD105, and CD73) and stem cell markers (KLF4, SOX2, and THY1). Similar to other MSCs, giant panda UC-MSCs have multiple differentiation ability and can differentiate into adipocytes, osteoblasts and chondrocytes. Giant panda UC-MSCs are new resources for basic research as cell models following their differentiation into different cell types and for future clinical treatments of giant panda diseases.

A novel missense mutation in the gene encoding major intrinsic protein (MIP) in a Giant panda with unilateral cataract formation

Background

Cataracts are defects of the lens that cause progressive visual impairment and ultimately blindness in many vertebrate species. Most cataracts are age-related, but up to one third have an underlying genetic cause. Cataracts are common in captive zoo animals, but it is often unclear whether these are congenital or acquired (age-related) lesions.

Results

Here we used a functional candidate gene screening approach to identify mutations associated with cataracts in a captive giant panda (Ailuropoda melanoleuca). We screened 11 genes often associated with human cataracts and identified a novel missense mutation (c.686G > A) in the MIP gene encoding major intrinsic protein. This is expressed in the lens and normally accumulates in the plasma membrane of lens fiber cells, where it plays an important role in fluid transport and cell adhesion. The mutation causes the replacement of serine with asparagine (p.S229N) in the C-terminal tail of the protein, and modeling predicts that the mutation induces conformational changes that may interfere with lens permeability and cell–cell interactions.

Conclusion

The c.686G > A mutation was found in a captive giant panda with a unilateral cataract but not in 18 controls from diverse regions in China, suggesting it is most likely a genuine disease-associated mutation rather than a single-nucleotide polymorphism. The mutation could therefore serve as a new genetic marker to predict the risk of congenital cataracts in captive giant pandas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07386-8.

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.

Lactobacillus plantarum BSGP201683 Improves the Intestinal Barrier of Giant Panda Microbiota-Associated Mouse Infected by Enterotoxigenic Escherichia coli K88

Giant pandas often suffered from gastrointestinal disease, especially the captive sub-adult one. Our study aims to investigate whether L. plantarum G83, a good panda-derived probiotic, can improve the intestinal barrier against the enterotoxigenic Escherichia coli K88 (E. coli K88) infection in giant panda microbiota-associated mice (GPAM). We treated SPF mice with antibiotics cocktail and transplanted the giant panda intestinal microbiota to set up a GPAM. Our results demonstrated that the microbiota of GPAM changed over time and was relatively stable in the short-term experiment (2-4 weeks). Whereafter, the GPAM pretreated with L. plantarum G83 for 15 days and infected with enterotoxigenic E. coli K88. The result indicated that the number of Bifidobacteria spp. increased in GPAM-G and GPAM-GE groups; the Lactobacillus spp. only increased in the GPAM-G group. Although the abundance of Enterobacteriaceae spp. only decreased in the GPAM-G group, the copy number of Escherichia coli in the GPAM-E group was significantly lower than that in the other groups. Meanwhile, the L. plantarum G83-induced alteration of microbiota could increase the mRNA expression of Claudin-1, Zo-1, and Occludin-1 in the GPAM-G group in the ileum; only Occludin-1 was increased in the GPAM-GE group. The sIgA in the ileum showed a positive response, also the result of body weight and histology in both the GPAM-G and GPAM-GE group. These results indicated that the L. plantarum G83 could improve the intestinal barrier to defense the enterotoxigenic E. coli K88 invasion.

from the giant panda, Ailuropoda melanoleuca, in China

BACKGROUND

Parasitic infections are among the important causes of death of giant pandas (Ailuropoda melanoleuca) that hamper their survival in the wild. There are about 35 species of parasites which have been identified in giant pandas, but no information is currently available regarding the infection of Babesia in giant pandas. Babesia spp. are common intraerythrocytic parasite in wildlife, transmitted by ixodid ticks, which cause babesiosis. Clinical signs of babesiosis include fever, hemolysis, anemia, jaundice and death.

METHODS

A species of Babesia was detected in the blood of a giant panda based on morphology and PCR amplification of the 18S rRNA gene. The phylogenetic relationship of Babesia sp. infecting giant panda was assessed by gene sequence alignment and phylogenetic analysis.

RESULTS

Our analysis revealed that the Babesia isolate detected was most similar to an unidentified species of Babesia identified in black bears (Ursus thibetanus japonicus) from Japan (Babesia sp. Iwate, AB586027.1) with a 99.56% sequence similarity, followed by Babesia sp. EBB (AB566229.1, 99.50%) and Babesia sp. Akita (AB566229.1, 99.07%).

CONCLUSIONS

To our knowledge, this is the first report of Babesia detected in the giant panda. The results indicate that this Babesia sp. may be a novel species, currently named Babesia sp. strain EBP01.

Identification and characterization of miRNA expression profiles across five tissues in giant panda

microRNA (miRNA) is a small endogenous noncoding RNA molecule that plays multiple roles in regulating most biological processes. However, for China's national treasure giant panda, a world-famous rare and protected species, reports of its miRNA have been found only in blood and breast milk. To explore the miRNA expression differences between different giant panda tissues, here, we generated the miRNA profiles of five tissues (heart, liver, spleen, lung, and kidney) from four giant pandas with Illumina Hiseq 2500 platform, and filtered the differentially expressed miRNAs (DEmiRs) in each tissue, predicted the target genes of miRNA from each tissue based on the DEmiRs. Then, the GO and KEGG enrichment analysis were conducted using the target genes predicted from DEmiRs in each tissue. The RNA-seq generated an average of 0.718 GB base per sample. A total of 1,256 known miRNAs and 12 novel miRNAs were identified, and there were 215, 131, 185, 83, and 126 tissue-specific DEmiRs filtered in the heart, liver, spleen, lung, and kidney, respectively, including miR-1b-5p, miR-122-5p, miR-143, miR-126-5p, and miR-10b-5p, respectively. The predicted target genes, including MYL2, LRP5, MIF, CFD, and PEBP1 in the heart, liver, spleen, lung, and kidney, respectively, were closely associated with tissue-specific biological functions. The enrichment analysis results of target genes showed tissue-specific characteristics, such as the significantly enriched GO terms extracellular matrix in the heart and insulin-like growth factor binding in the liver. The miRNA profiles of the heart, liver, spleen, lung, and kidney of giant panda have been reported in this study, it reveals the miRNA expression differences between different tissues of the giant panda, and provides valuable genetic resources for the further related molecular genetic research of the rare and protected species giant panda and other mammals.

Dynamics of bacteriophages in gut of giant pandas reveal a potential regulation of dietary intake on bacteriophage composition

Fecal samples of cubs and adults of giant pandas were examined to determine the effects of diets on the diversity and dynamics of gut bacteriophages. Enterobacteria phage, Salmonella phage, Escherichia phage, Shigella phage, Klebsiella phage, and Lactococcus phage were found to be dominant in both cub and adult samples. Citrobacter phage, Cronobacter phage, Pectobacterium phage, Erwinia phage, Dickeya phage, Erwinia phage, Enterococcus phage, and Pseudomonas phage were more abundant in adults, while Lactococcus phage, Streptococcus phage, Lactobacillus phage, and Leuconostoc phage were more abundant in cubs. The abundance and diversity of the majority of phage species were increased in pandas with age. There was an increase in the abundance of Pectobacterium phage and a decrease in the abundance of Lactobacillus phage, Leuconostoc phage, Bacillus phage, and Streptococcus phage in adults. As cubs and adults of giant pandas have different dietary habits, these observations suggest a significant effect of diets on the composition and abundance of gut bacteriophages in giant pandas.

Review on parasites of wild and captive giant pandas (Ailuropoda melanoleuca): Diversity, disease and conservation impact

The giant panda (Ailuropoda melanoleuca) is a rare species with a small global population size, and lives in the wild in only a few fragmented mountain ranges of Southwest China. Parasitic infections are among the important causes of death of giant pandas that hamper their group development. We reviewed the parasitic infections prevailing in giant pandas, and the parasitic diversity, diseases and their impact on conservation of this animal. A total of 35 parasitic species were documented in giant pandas, belonging to nematode (n = 6), trematode (n = 1), cestode (n = 2), protozoa (n = 9), and ectozoa (n = 17 (tick = 13, mite = 2, and flea = 2)). Among them, Baylisascaris schroederi had the highest prevalence and was the leading cause of death for giant pandas. Some parasites caused asymptomatic infections in giant pandas, and their health implications for the pandas remain unknown. As a whole, parasites are reported to be an important threat to the conservation of the giant pandas. Regular deworming and environmental disinfection appear to be effective ways to prevent captive giant pandas from parasitoses. In wild panda populations, parasitic control measures are suggested to include detailed examination of the ecology of the host-parasite assembly, with particular attention to density-dependent transmission. The parasitic pathogenesis and detection methods together with their biology, epidemiology, treatment, prevention and control need to be further studied for better protection of giant pandas from parasitoses.

Site-directed mutagenesis of odorant-binding proteins

Modifying the affinity of odorant-binding proteins (OBPs) to small ligands by replacement of specific residues in the binding pocket may lead to several technological applications. Thanks to their compact and stable structures, OBPs are currently regarded as the best candidates to be used in biosensing elements for odorants and volatiles detection. The wide and rich information on the structure of these proteins both in their apo-forms and in complexes with specific ligands provides guidelines to design reliable mutants to monitor specific targets. The same engineered proteins may also find applications in the slow release of pheromones and other chemicals in the environment, as well as in the fine purification of drugs, including the resolution of racemates. Apart from such useful applications, site-directed mutagenesis represents an interesting approach to dissect the specific interactions between small chemicals and amino acid residues in the binding pocket. These studies can lead to design of better ligands, such as pheromone analogues with desired physico-chemical characteristics. In this chapter we examine the different uses of mutagenesis applied to OBPs and report a couple of protocols that have been successful in our hands.

Effects of extender type on the quality of post-thaw giant panda (Ailuropoda melanoleuca) semen

Sperm cryopreservation is an essential approach for assisted reproduction and genetic resources conservation in captive giant pandas. Cryopreservation, however, leads to a significant decrease in sperm quality and, consequently, a low fertilization rate. Therefore, it is mandatory to disclose more suitable and efficient freezing strategies for sperm cryopreservation. In the present study, we compared for the first time the performance of two commercial freeze extender (INRA96 versus TEST) freezing methods on post-thawed semen quality. Semen cryopreserved with the INRA96 showed better total motility (73.00 ± 4.84% vs 57.56 ± 3.60%, P < 0.001), membrane integrity (60.92 ± 2.27% vs 40.53 ± 2.97%, P < 0.001) and acrosome integrity (90.39 ± 2.74% vs 84.26 ± 4.27%, P < 0.05) than stored with TEST. There was no significant difference in DNA integrity after thawing between the two extenders (95.69 ± 3.60% vs 94.26 ± 4.84%). In conclusion, the INRA96 method showed to be better for giant panda sperm cryopreservation and should therefore be recommended for use in order to increase success of artificial insemination.

Isolation and Preliminary Screening of a Weissella confusa Strain from Giant Panda (Ailuropoda melanoleuca)

Weissella confusa has recently received attention for its probiotic potential. Some W. confusa and Weissella cibaria strains isolated from fermented foods show favorable probiotic effects. However, the probiotic properties of W. confusa isolated from giant panda remain unreported to date. Thus, this study isolated a W. confusa strain from giant panda feces and then investigated its characteristics and probiotic properties. A lactic acid bacteria strain was isolated from giant panda fecal samples. The isolated strain was screened by in vitro probiotic property tests, including in vitro antimicrobial test, antioxidant test, surface hydrophobicity, and stress resistance. On the basis of biochemical identification and 16S rDNA sequencing, the W. confusa strain was identified as BSP201703. This Weissella confusa strain can survive at pH 2 and 0.3% (w/v) concentration of bile salt environment and inhibit common intestinal pathogens. It also possesses an in vitro antioxidant capacity, a high auto-aggregation ability, and a high surface hydrophobicity. BSP201703 might serve as a probiotic to giant pandas.

Data on diet and growth by giant panda in zoo Negara, Malaysia

In this data article we present the determinations of the diet preference and growth of a pair of the giant panda, Ailuropoda melanoleuca (David, 1869) from Zoo Negara Malaysia. Once considered as endangered, the captive giant pandas were given with nine species of local bamboo in separate indoor enclosures. We recorded data between May 25, 2014 and December 31, 2016 and analysed it based on food preference, the pattern toward food consumption and body weights using SPSS v25.0 (IBM, USA). Data on the bamboo preference, daily average bamboo provided and consumed, and factors predicting of body weight per individual are reported in this article. The data highlight correlation between panda growth (kg) to the part of bamboo consumed (kg) and exhibit the pattern of preferred part of food (i.e.: either the leaf, culm or shoots of bamboo variety) for panda consumptions. The food consumption toward the body weight was modelled using logistic regression analysis to help determine the pattern of food consumption and body weight of giant panda in the future and based on regression model 1, only consumed variable is significance to the model.

First demonstration of giant panda's immune response to canine distemper vaccine

The Canine Distemper Virus (CDV) is a high fatal virus to the giant panda (Ailuropoda melanoleuca), where CDV vaccination is a key preventative measure in captive giant pandas. However, the immune response of giant pandas to CDV vaccination has been little studied. In this study, we investigated the blood transcriptome expression profiles of five giant panda cubs after three inoculations, 21 days apart. Blood samples were collected before vaccination (0 Day), and 24 h after each of the three inoculations; defined here as 1 Day, 21  Day, and 42  Day. Compared to 0 Day, we obtained 1262 differentially expressed genes (DEGs) during inoculations. GO and KEGG pathways enrichment analysis of these DEGs found 222 GO terms and 40 pathways. The maximum immune-related terms were enriched by DEGs from comparisons of 21  Day and 0 Day. In the PPI analysis, we identified RSAD2, IL18, ISG15 immune-related hub genes from 1 Day and 21 Day comparison. Compared to 0 Day, innate immune-related genes, TLR4 and TLR8, were up-regulated at 1 Day, and the expressions of IRF1, RSAD2, MX1, and OAS2 were highest at 21  Day. Of the adaptive immune-related genes, IL15, promoting T cell differentiation into CD8⁺T cells, was up-regulated after the first two inoculations, IL12β, promoting T cell differentiation into memory cells, and IL10, promoting B cell proliferation and differentiation, were down-regulated during three inoculations. Our results indicated that the immune response of five giant panda cubs was strongest after the second inoculation, most likely protected against CDV infection through innate immunity and T cells, but did not produce enough memory cells to maintain long-term immunity after CDV vaccination.