Metagenomic profiling of cecal microbiota and antibiotic resistome in rodents
Kai-Meng Shang, Hany M. Elsheikha, He Ma, Yongjie Wei, Jixin Zhao, Ya Qin, Jianming Li, Ziyu Zhao, Xiao‐Xuan Zhang
Abstract
The rodent gut microbiota is a known reservoir of antimicrobial resistance, yet the distribution of antibiotic resistance genes (ARGs) within rodent cecal microbial communities and the specific bacterial species harboring these ARGs remain largely underexplored. This study employed high-throughput sequencing of 122 samples from five distinct rodent species to comprehensively profile the diversity and distribution of ARGs and to identify the bacterial hosts of these genes. A gene catalog of the rodent cecal microbiome was constructed, comprising 22,757,369 non-redundant genes. Analysis of the microbial composition and diversity revealed that Bacillota and Bacteroidota were the dominant bacterial phyla across different rodent species, with significant variations in species composition among the rodents. In total, 3703 putative antimicrobial resistance protein-coding genes were identified, corresponding to 392 unique ARG types classified into 32 resistance classes. The most enriched ARGs in the rodent cecal microbiome were associated with multidrug resistance, followed by glycopeptide and elfamycin antibiotics. Procrustes analysis demonstrated a correlation between the structure of the microbial community and the resistome. Metagenomic assembly-based host tracking indicated that most ARG-carrying contigs originated from the bacterial family Oscillospiraceae . Additionally, 130 ARGs showed significant correlations with mobile genetic elements. These findings provide new insights into the cecal microbiota and the prevalence of ARGs across five rodent species. Future research on a wider range of wild rodent species carrying ARGs will further elucidate the mechanisms underlying the transmission of antimicrobial resistance. • High- throughput sequencing of 122 samples from cecum of 5 rodent species. • Gene catalog includes ∼22.8 M non-redundant genes: 21.3 % bacterial, 0.1 % archaea, 8.9 % viruses, 3.0 % fungi. • Bacillota and Bacteroidota are dominant phyla; significant variation in species composition observed. • Wild rodents showed higher Bacteroidota levels; significant differences in microbiota compared to lab mice. • 3703 antibiotic resistance genes (ARGs) across 392 unique types, and 130 ARGs correlated with mobile genetic elements.