Diversity, functions, and antibiotic resistance genes of bacteria and fungi are examined in the bamboo plant phyllosphere that serve as food for the giant pandas

The phyllosphere of bamboo is rich in microorganisms that can disrupt the intestinal microbiota of the giant pandas that consume them, potentially leading to their death. In the present study, the abundance, diversity, biological functions (e.g., KEGG and CAZyme), and antibiotic resistance genes (ARGs) of bacteria and fungi in two bamboo species phyllosphere (Chimonobambusa szechuanensis, CS; Bashania fangiana, BF) in Daxiangling Nature Reserve (an important part of the Giant Panda National Park) were investigated respectively by amplicon sequencing of the whole 16S rRNA and ITS1-ITS2 genes on PacBio Sequel and whole-metagenome shotgun sequencing on Illumina NovaSeq 6000 platform. The results suggested that there were respectively 18 bacterial and 34 fungi biomarkers between the phyllosphere of the two species of bamboo. Beta diversity of bacteria and fungi communities exited between the two bamboos according to the (un)weighted UniFrac distance matrix. Moreover, the functional analysis showed that the largest relative abundance was found in the genes related to metabolism and global and overview maps. Glycoside hydrolases (GHs) and glycosyl transferases (GTs) have a higher abundance in two bamboo phyllospheres. Co-occurrence network modeling suggested that bacteria and fungi communities in CS phyllosphere employed a much more complex metabolic network than that in BF, and the abundance of multidrug, tetracycline, and glycopeptide resistance genes was higher and closely correlated with other ARGs. This study references the basis for protecting bamboo resources foraged by wild giant pandas and predicts the risk of antibiotic resistance in bamboo phyllosphere bacterial and fungal microbiota in the Giant Panda National Park, China.