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.

Sexual Maturity and Life Stage Influences Toxic Metal Accumulation in Croatian Brown Bears.

The influence of reproductive and (early) life stages on toxic metal levels was investigated in the brown bear (Ursus arctos), the largest mammalian predator species in Croatia. The purpose was to examine critical clusters in a population that might be at a higher risk of adverse health effects caused by metals as environmental contaminants. Levels of cadmium, mercury and lead in muscle, liver and kidney cortex of 325 male and 139 female bears, quantified by inductively coupled plasma mass spectrometry, were analysed according to distinct bear life stages (young: cub, yearling, subadult; mature: adult). Metal levels did not differ among sexes in young animals (< 4 years), except for mercury in muscles (higher in females), and adult females had higher cadmium and mercury. A trend of renal cadmium accumulation with age in immature male animals disappeared once they reached maturity, whereas for females this trend has only slowly declined in mature compared to immature bears. In early life stage (< 1 year), bear cubs had lower cadmium, comparable mercury, and higher lead in the kidneys than the bears of the following age category (yearlings). Due to a higher proportion of renal lead transfer from the mother to the cub compared with cadmium, it may be that the high burden of cadmium found in kidneys of older females has lower toxicological concern for their cubs than the lead content. Sex, reproductive, and life stages of bears were confirmed as important in assessing toxic metal burden.

Growth trajectories in the cave bear and its extant relatives: an examination of ontogenetic patterns in phylogeny.

BACKGROUND: The study of postnatal ontogeny can provide insights into evolution by offering an understanding of how growth trajectories have evolved resulting in adult morphological disparity. The Ursus lineage is a good subject for studying cranial and mandibular shape and size variation in relation to postnatal ontogeny and phylogeny because it is at the same time not diverse but the species exhibit different feeding ecologies. Cranial and mandibular shapes of Ursus arctos (brown bear), U. maritimus (polar bear), U. americanus (American black bear), and the extinct U. spelaeus (cave bear) were examined, using a three-dimensional geometric morphometric approach. Additionally, ontogenetic series of crania and mandibles of U. arctos and U. spelaeus ranging from newborns to senile age were sampled. RESULTS: The distribution of specimens in morphospace allowed to distinguish species and age classes and the ontogenetic trajectories U. arctos and U. spelaeus were found to be more similar than expected by chance. Cranial shape changes during ontogeny are largely size related whereas the evolution of cranial shape disparity in this clade appears to be more influenced by dietary adaptation than by size and phylogeny. The different feeding ecologies are reflected in different cranial and mandibular shapes among species. CONCLUSIONS: The cranial and mandibular shape disparity in the Ursus lineage appears to be more influenced by adaptation to diet than by size or phylogeny. In contrast, the cranial and mandibular shape changes during postnatal ontogeny in U. arctos and U. spelaeus are probably largely size related. The patterns of morphospace occupation of the cranium and the mandible in adults and through ontogeny are different.

Estimation of the growth curve and heritability of the growth rate for giant panda (Ailuropoda melanoleuca) cubs.

Giant panda cubs have a low survival rate during the newborn and early growth stages. However, the growth and developmental parameters of giant panda cubs during the early lactation stage (from birth to 6 months) are not well known. We examined the growth and development of giant panda cubs by the Chapman growth curve model and estimated the heritability of the maximum growth rate at the early lactation stage. We found that 83 giant panda cubs reached their maximum growth rate at approximately 75-120 days after birth. The body weight of cubs at 75 days was 4285.99 g. Furthermore, we estimated that the heritability of the maximum growth rate was moderate (h(2) = 0.38). Our study describes the growth and development of giant panda cubs at the early lactation stage and provides valuable growth benchmarks. We anticipate that our results will be a starting point for more detailed research on increasing the survival rate of giant panda cubs. Feeding programs for giant panda cubs need further improvement.

[Intestinal fungal diversity of sub-adult giant panda].

OBJECTIVE: The fungi diversity in the guts of five sub-adult giant pandas was analyzed. METHOD: We analyzed the fungal internal transcribed spacer sequences (ITS) using restriction fragment length polymorphism (RFLP). ITS regions were amplified with fungal universal primers to construct ITS clone libraries. The fingerprints were analyzed by restriction fragment length polymorphism using the Hha I and Hae III enzymes. The cloned PCR products were analyzed by sequencing and diversities were demonstrated by phylogenetic tree. RESULTS: The gut fungi of 5 sub-adult giant pandas were mainly composed of Ascomycota (average of 46.24%), Basidiomycota ( average of 15.79%), unclassified (average of 29.14%), uncultured fungus (average of 8.83% ). Ascomycota was mainly composed of Saccharomycetes (average of 63.74%) and Dothideomycetes ( average of 35.91%); Basidiomycota was mainly composed of Tremellomycetes (average of 65.80%) and Microbotryomycetes (average of 33.15%). Four classes were mainly composed of Candida and Debaryomyces; Pleosporales and Myriangium; Cystofilobasidium and Trichosporon; Leucosporidium, and Leucosporidiella, whereas the proportions were different for each sample. CONCLUSION: Fungal flora existing in the intestines of sub-adult giant pandas expand our knowledge on the structure of the giant panda gut microbes and also help us to further study whether fungal flora can help giant pandas digest high-fiber foods.

Environmental, biological and anthropogenic effects on grizzly bear body size: temporal and spatial considerations.

BACKGROUND: Individual body growth is controlled in large part by the spatial and temporal heterogeneity of, and competition for, resources. Grizzly bears (Ursus arctos L.) are an excellent species for studying the effects of resource heterogeneity and maternal effects (i.e. silver spoon) on life history traits such as body size because their habitats are highly variable in space and time. Here, we evaluated influences on body size of grizzly bears in Alberta, Canada by testing six factors that accounted for spatial and temporal heterogeneity in environments during maternal, natal and 'capture' (recent) environments. After accounting for intrinsic biological factors (age, sex), we examined how body size, measured in mass, length and body condition, was influenced by: (a) population density; (b) regional habitat productivity; (c) inter-annual variability in productivity (including silver spoon effects); (d) local habitat quality; (e) human footprint (disturbances); and (f) landscape change. RESULTS: We found sex and age explained the most variance in body mass, condition and length (R(2) from 0.48-0.64). Inter-annual variability in climate the year before and of birth (silver spoon effects) had detectable effects on the three-body size metrics (R(2) from 0.04-0.07); both maternal (year before birth) and natal (year of birth) effects of precipitation and temperature were related with body size. Local heterogeneity in habitat quality also explained variance in body mass and condition (R(2) from 0.01-0.08), while annual rate of landscape change explained additional variance in body length (R(2) of 0.03). Human footprint and population density had no observed effect on body size. CONCLUSIONS: These results illustrated that body size patterns of grizzly bears, while largely affected by basic biological characteristics (age and sex), were also influenced by regional environmental gradients the year before, and of, the individual's birth thus illustrating silver spoon effects. The magnitude of the silver spoon effects was on par with the influence of contemporary regional habitat productivity, which showed that both temporal and spatial influences explain in part body size patterns in grizzly bears. Because smaller bears were found in colder and less-productive environments, we hypothesize that warming global temperatures may positively affect body mass of interior bears.

Age-structured Jolly-Seber model expands inference and improves parameter estimation from capture-recapture data.

Understanding the influence of individual attributes on demographic processes is a key objective of wildlife population studies. Capture-recapture and age data are commonly collected to investigate hypotheses about survival, reproduction, and viability. We present a novel age-structured Jolly-Seber model that incorporates age and capture-recapture data to provide comprehensive information on population dynamics, including abundance, age-dependent survival, recruitment, age structure, and population growth rates. We applied our model to a multi-year capture-recapture study of polar bears (Ursus maritimus) in western Hudson Bay, Canada (2012-2018), where management and conservation require a detailed understanding of how polar bears respond to climate change and other factors. In simulation studies, the age-structured Jolly-Seber model improved precision of survival, recruitment, and annual abundance estimates relative to standard Jolly-Seber models that omit age information. Furthermore, incorporating age information improved precision of population growth rates, increased power to detect trends in abundance, and allowed direct estimation of age-dependent survival and changes in annual age structure. Our case study provided detailed evidence for senescence in polar bear survival. Median survival estimates were lower (<0.95) for individuals aged <5 years, remained high (>0.95) for individuals aged 7-22 years, and subsequently declined to near zero for individuals >30 years. We also detected cascading effects of large recruitment classes on population age structure, which created major shifts in age structure when these classes entered the population and then again when they reached prime breeding ages (10-15 years old). Overall, age-structured Jolly-Seber models provide a flexible means to investigate ecological and evolutionary processes that shape populations (e.g., via senescence, life expectancy, and lifetime reproductive success) while improving our ability to investigate population dynamics and forecast population changes from capture-recapture data.

Reduced body size and cub recruitment in polar bears associated with sea ice decline.

Rates of reproduction and survival are dependent upon adequate body size and condition of individuals. Declines in size and condition have provided early indicators of population decline in polar bears (Ursus maritimus) near the southern extreme of their range. We tested whether patterns in body size, condition, and cub recruitment of polar bears in the southern Beaufort Sea of Alaska were related to the availability of preferred sea ice habitats and whether these measures and habitat availability exhibited trends over time, between 1982 and 2006. The mean skull size and body length of all polar bears over three years of age declined over time, corresponding with long-term declines in the spatial and temporal availability of sea ice habitat. Body size of young, growing bears declined over time and was smaller after years when sea ice availability was reduced. Reduced litter mass and numbers of yearlings per female following years with lower availability of optimal sea ice habitat, suggest reduced reproductive output and juvenile survival. These results, based on analysis of a long-term data set, suggest that declining sea ice is associated with nutritional limitations that reduced body size and reproduction in this population.

Free-roaming dogs limit habitat use of giant pandas in nature reserves.

Giant pandas (Ailuropoda melanoleuca) were historically hunted using dogs and are currently threatened by free-roaming dogs and their associated diseases. To better understand the spatial magnitude of this threat, we used a GIS approach to investigate edge effects of dogs on giant panda habitat. We first examined two nature reserves with contrasting free-roaming dog populations: Liziping, with many dogs (~0.44/km2), and Daxiangling, with few dogs (~0.14/km2). Spatial analysis indicated that giant pandas at Liziping (but not Daxiangling) showed a shift in habitat use away from populated areas consistent with a risk response to the foray distance of free-roaming dogs (10.9 km path-distance). Most giant panda locations (86%) from the 2014 census in Liziping were clustered around remote "dog-free zones." Expanding this analysis across the entire giant panda range revealed that 40% of panda habitat is within the foray distance of dogs. Our assessment will inform dog control programs including monitoring, education, veterinary care, and other measures. We recommend that reserves designated for the release of translocated pandas receive priority consideration for dog control efforts. Only by understanding and managing complex interactions between humans, domestic animals, and wild animals can we sustain natural systems in a world increasingly dominated by humans.

Identifying reliable indicators of fitness in polar bears.

Animal structural body size and condition are often measured to evaluate individual health, identify responses to environmental change and food availability, and relate food availability to effects on reproduction and survival. A variety of condition metrics have been developed but relationships between these metrics and vital rates are rarely validated. Identifying an optimal approach to estimate the body condition of polar bears is needed to improve monitoring of their response to decline in sea ice habitat. Therefore, we examined relationships between several commonly used condition indices (CI), body mass, and size with female reproductive success and cub survival among polar bears (Ursus maritimus) measured in two subpopulations over three decades. To improve measurement and application of morphometrics and CIs, we also examined whether CIs are independent of age and structural size-an important assumption for monitoring temporal trends-and factors affecting measurement precision and accuracy. Maternal CIs and mass measured the fall prior to denning were related to cub production. Similarly, maternal CIs, mass, and length were related to the mass of cubs or yearlings that accompanied her. However, maternal body mass, but not CIs, measured in the spring was related to cub production and only maternal mass and length were related to the probability of cub survival. These results suggest that CIs may not be better indicators of fitness than body mass in part because CIs remove variation associated with body size that is important in affecting fitness. Further, CIs exhibited variable relationships with age for growing bears and were lower for longer bears despite body length being related to cub survival and female reproductive success. These results are consistent with findings from other species indicating that body mass is a useful metric to link environmental conditions and population dynamics.

Gene expression profiles during postnatal development of the liver and pancreas in giant pandas.

Giant pandas are unique Carnivora with a strict bamboo diet. To investigate the molecular mechanism of giant panda nutrient metabolism from newborn to adult, the gene expression profiles of giant panda liver and pancreas tissues collected from three important feeding stages were investigated using RNA-seq. We found a total of 3,211 hepatic and 3,343 pancreatic differentially expressed genes (DEGs) from three comparisons between suckling and no feeding, adult and no feeding, and adult and suckling groups. Few differences in gene-expression profiles were exhibited between no feeding and suckling groups in both tissues. GO and KEGG analyses were performed to further understand the biological functions of the DEGs. In both the liver and pancreas, genes related mainly to cell cycle processes were highly up-regulated in newborn samples whereas genes related to metabolism and immunity were up-regulated in adult giant pandas. The high expression of metabolism-related genes in adult samples probably helps to fulfill the metabolic function requirements of the liver and pancreas. In contrast, several vital genes involved in cholesterol metabolism and protein digestion and absorption were over-expressed in newborn samples. This may indicate the importance of cholesterol metabolism and protein digestion and absorption processes in giant panda infancy.

Splitting hairs: dietary niche breadth modelling using stable isotope analysis of a sequentially grown tissue.

Stable isotope data from durable, sequentially grown tissues (e.g. hair, claw, and baleen) is commonly used for modelling dietary niche breadth. The use of tissues grown over multiple months to years, however, has the potential to complicate isotopic niche breadth modelling, as time-averaged stable isotope signals from whole tissues may obscure information available from chronologically resolved stable isotope signals in serially sectioned tissues. We determined if whole samples of brown bear guard hair produced different isotopic niche breadth estimates than those produced from subsampled, serially sectioned samples of the same tissue from the same set of individuals. We sampled guard hair from brown bears (Ursus arctos) in four regions of Alaska with disparate biogeographies and dietary resource availability. Whole hair and serially sectioned hair samples were used to produce paired isotopic dietary niche breadth estimates for each region in the SIBER Bayesian model framework in R. Isotopic data from serially sectioned hair consistently produced larger estimates of isotopic dietary niche breadth than isotope data from whole hair samples. Serial sampling captures finer-scale changes in diet and when cumulatively used to estimate isotopic niche breadth, the serially sampled isotope data more fully captures dietary variability and true isotopic niche breadth.

Change in body weight of mothers and neonates and in milk composition during denning period in captive Japanese black bears (Ursus thibetanus japonicus).

Japanese black bears, Ursus thibetanus japonicus, have been classified as a vulnerable species so that data on reproduction are needed to maintain and/or extend their population. They are known to have a peculiar style of reproduction, giving birth to their neonates and raising them during denning, a period of complete fasting. In this study, we investigated the metabolic rate and milk composition of mother bears raising neonates, and the changes in body weight of the neonates under captive conditions. Seven female bears kept in dens were weighed once a month, and the amount of energy they used was calculated. From birth, cubs were also weighed and their growth rate was determined. In addition, the milk composition was analyzed to investigate its characteristics. As a result, it was found that mother bears used 34% more energy than did solitary females. There was no significant difference in the energy used for nursing whether they had single or twin cubs. On the other hand, the body weight gain of single cubs was significantly higher than that of twin cubs, suggesting that the growth of the cubs was highly affected by the suppression of mother's energy consumption during the fasting period. The milk had high fat and low sugar concentrations. This indicates that mother bears used the fat accumulated prior to denning for their main energy source when raising cubs. Considering all results together, Japanese black bears showed remarkable efficiency in the use of energy for reproduction during the fasting period.

Evaluating giant panda as a surrogate species for conservation co-occurring species in the Baishuijiang National Nature Reserve.

The establishment of nature reserves is a key approach for biodiversity conservation worldwide. However, the effectiveness of nature reserves established by protecting the habitat needs of surrogate species is questioned. In this study, the Baishuijiang National Nature Reserve (Baishuijiang NNR), located in the Minshan Mountains, China, which is established mainly for the conservation of giant panda (a surrogate for the conservation of other endangered species) was selected. We quantitatively evaluated the conservation effectiveness of the reserve for giant panda and co-occurring species (here, seven protected species) using a maximum entropy model (Maxent), and analyzed spatial congruence between giant panda and other seven species. Results shown that the habitat of giant panda generally included the habitat of other seven protected species, suggesting that conservation of giant panda habitat also allows the conservation for the habitat of almost co-occurring species. Hence, the natural reserve established for giant panda as a surrogate species has a relatively high effectiveness. A high proportion of the suitable habitat for six species is inside the core zone, but a high proportion of the suitable habitat for two species is located in the experimental and buffer zones. Thus, the two species are affected by human activities. To improve the conservation effectiveness of the nature reserve, the management zones need to be amended. The result of the study will be beneficial for future conservation and management of the reserve. This study provides an effective method for evaluating the conservation effectiveness of nature reserves in other area of the worldwide.

Does hard mast production affect patterns of cementum annuli formation in premolar teeth of Asian black bears (Ursus thibetanus)?

Cementum annuli widths in mammals are is influenced by the nutrition of mammals. Reproductive stress has been is suggested to reduce the width of lead to narrower cementum annuli widths in female Asian black bears (Ursus thibetanus); however, food availability in autumn strongly impacts bear nutrition and likely impacts cementum widths as well. This study aimed to test how cementum annuli widths and the formation of false annuli were influenced by hard mast production. We established two hypotheses: (1) cementum annuli widths become narrower in poor mast years owing to inadequate nutritional conditions and (2) false annuli occur more frequently in poor mast years. We used teeth samples from male bears to avoid reproductive influences and separated width data into "adult" and "subadult" groups. We calculated the proportional width index (PWI) and used linear mixed models to estimate the masting effects on PWI. Generalized linear mixed models estimated the masting effects on false annuli frequency. True annuli widths and false annuli formation showed no significant relationship with mast production in adults. In subadults, poor mast production weak negative influence on false annuli formation. These new data resolve previous questions, allowing us to deduce that cementum annuli widths are a reliable index of reproductive success in female bears.

Using bear rub data and spatial capture-recapture models to estimate trend in a brown bear population.

Trends in population abundance can be challenging to quantify during range expansion and contraction, when there is spatial variation in trend, or the conservation area is large. We used genetic detection data from natural bear rubbing sites and spatial capture-recapture (SCR) modeling to estimate local density and population growth rates in a grizzly bear population in northwestern Montana, USA. We visited bear rubs to collect hair in 2004, 2009-2012 (3,579-4,802 rubs) and detected 249-355 individual bears each year. We estimated the finite annual population rate of change 2004-2012 was 1.043 (95% CI = 1.017-1.069). Population density shifted from being concentrated in the north in 2004 to a more even distribution across the ecosystem by 2012. Our genetic detection sampling approach coupled with SCR modeling allowed us to estimate spatially variable growth rates of an expanding grizzly bear population and provided insight into how those patterns developed. The ability of SCR to utilize unstructured data and produce spatially explicit maps that indicate where population change is occurring promises to facilitate the monitoring of difficult-to-study species across large spatial areas.

Giant panda protection: challenges and hopes.

In this paper, we discussed some major issues that hinder giant panda protection, such as diverse and mixed threats, habitat fragmentation, as well as the survey method to be improved, and some new protection actions appeared, such as the pilot program for the giant panda national park system and the administration of the giant panda national park was established. These information could provide important information for giant panda protection.

Role of nature reserves in giant panda protection.

Giant panda (Ailuropoda melanoleuca) is a flagship species in nature conservation of the world; to protect this species, 67 nature reserves have been established in China. To evaluate the protection effect of giant panda nature reserves, we analyzed the variation of giant panda number and habitat area of 23 giant panda nature reserves of Sichuan province based on the national survey data released by State Forestry Administration and Sichuan Forestry Department. Results showed that from the third national survey to the fourth, giant panda number and habitat area of 23 giant panda nature reserves of Sichuan province failed to realize the significant increase. Furthermore, we found that the total population growth rate of 23 nature reserves in the last 12 years was lower than those of the province total of Sichuan and the national total of China, and the total habitat area of the 23 nature reserves was decreasing in the last 12 years, but the province total and national total were all increasing. We propose that giant panda protection should pay more attention to how to improve the protective effects of nature reserves.

A Bayesian state-space model using age-at-harvest data for estimating the population of black bears (Ursus americanus) in Wisconsin.

Population estimation is essential for the conservation and management of fish and wildlife, but accurate estimates are often difficult or expensive to obtain for cryptic species across large geographical scales. Accurate statistical models with manageable financial costs and field efforts are needed for hunted populations and using age-at-harvest data may be the most practical foundation for these models. Several rigorous statistical approaches that use age-at-harvest and other data to accurately estimate populations have recently been developed, but these are often dependent on (a) accurate prior knowledge about demographic parameters of the population, (b) auxiliary data, and (c) initial population size. We developed a two-stage state-space Bayesian model for a black bear (Ursus americanus) population with age-at-harvest data, but little demographic data and no auxiliary data available, to create a statewide population estimate and test the sensitivity of the model to bias in the prior distributions of parameters and initial population size. The posterior abundance estimate from our model was similar to an independent capture-recapture estimate from tetracycline sampling and the population trend was similar to the catch-per-unit-effort for the state. Our model was also robust to bias in the prior distributions for all parameters, including initial population size, except for reporting rate. Our state-space model created a precise estimate of the black bear population in Wisconsin based on age-at-harvest data and potentially improves on previous models by using little demographic data, no auxiliary data, and not being sensitive to initial population size.

Age-associated microbiome shows the giant panda lives on hemicelluloses, not on cellulose.

The giant panda feeds almost exclusively on bamboo, a diet highly enriched in lignin and cellulose, but is characterized by a digestive tract similar to carnivores. It is still large unknown if and how the giant panda gut microbiota contributes to lignin and cellulose degradation. Here we show the giant pandas' gut microbiota does not significantly contribute to cellulose and lignin degradation. We found that no operational taxonomic unit had a nearest neighbor identified as a cellulolytic species or strain with a significant higher abundance in juvenile than cubs, a very low abundance of putative lignin and cellulose genes existed in part of analyzing samples but a significant higher abundance of genes involved in starch and hemicellulose degradation in juveniles than cubs. Moreover, a significant lower abundance of putative cellulolytic genes and a significant higher abundance of putative α-amylase and hemicellulase gene families were present in giant pandas than in omnivores or herbivores.