DNA was extracted from the fecal samples using the Mo Bio PowerFecal DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA, USA) according to the manufacturer’s instructions. Variable region 4 (V4) forward primer: GTGYCAGCMGCCGCGGTAA, reverse primer: GGACTACHVGGGTWTCTAAT. The PCR reaction (50 µl total volume) contained 2 µL DNA (20 ng), 19 µL PCR grade water, 25 µL 2× Es Taq Master Mix (CW BIO, Beijing, China) (including 2× Es Taq Polymerase, 0.075 µM Mg2+, and 10 µM dNTP mix), 2 µL (0.4 µM) forward primer, and 2 µL (0.4 µM) reverse primer. PCR was performed at 94 °C for 1 min, followed by 30 cycles of 94 °C for 20 s, 59 °C for 25 s, and 68 °C for 45 s, followed by a final extension at 68 °C for 10 min. Variable region 4 was sequenced at the Beijing Genomics Institute (Beijing, China).
Shotgun metagenomic sequencing was performed at Novogene (Beijing, China). DNA libraries were constructed according to Illumina’s instructions. Briefly, DNA was sheared to 300–400 bp fragments, and the DNA from each individual sample was barcoded uniquely. Since most of the demographic data for the wild pandas are unknown, we intended to choose a wider range of samples from the captive samples to see if the environment is still the major driver of the gut metagenome. Therefore, we randomly chose three cubs (≤2-year-olds), one sub-adult (>2-year-olds and <5-year-olds) and three adults (≥5-year-olds and ≤20-year-olds) (including three males and four females) for metagenome analysis (
The datasets used in this study are accessible from the National Centre for Biotechnology Information Sequence Read Archive (SRA;