Metagenomic analysis of gut microbiome and resistome of diarrheal fecal samples from Kolkata, India, reveals the core and variable microbiota including signatures of microbial dark matter

BACKGROUND: Metagenomic analysis of the gut microbiome and resistome is instrumental for understanding the dynamics of diarrheal pathogenesis and antimicrobial resistance transmission (AMR). Metagenomic sequencing of 20 diarrheal fecal samples from Kolkata was conducted to understand the core and variable gut microbiota. Five of these samples were used for resistome analysis. The pilot study was conducted to determine a microbiota signature and the source of antimicrobial resistance genes (ARGs) in the diarrheal gut. RESULTS: 16S rRNA amplicon sequencing was performed using Illumina MiSeq platform and analysed using the MGnify pipeline. The Genome Taxonomy Database (GTDB-Tk) was used for bacterial taxonomic identification. Diarrheal etiology was determined by culture method. Phylum Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria were consistently present in 20 samples. Firmicutes was the most abundant phylum in 11 samples. The Bacteroidetes/Firmicutes ratio was less than 1 in 18 samples. 584 genera were observed. 18 of these were present in all the 20 samples. Proteobacteria was the dominant phylum in 6 samples associated with Vibrio cholerae infection. Conservation of operational taxonomic units (OTUs) among all the samples indicated the existence of a core microbiome. Asymptomatic carriage of pathogens like Vibrio cholerae and Helicobacter pylori was found. Signature of Candidate phyla or “microbial dark matter” occurred. Significant correlation of relative abundance of bacterial families of commensals and pathogens were found. Whole-genome sequencing (WGS) on Illumina MiSeq system and assembly of raw reads using metaSPAdes v3.9.1 was performed to study the resistome of 5 samples. ABRicate was used to assign ARG function. 491 resistance determinants were identified. In 80% of the samples tetracycline resistance was the most abundant resistance determinant. High abundance of ARGs against β-lactams, aminoglycosides, quinolones and macrolides was found. Eschericia sp. was the major contributor of ARGs. CONCLUSIONS: This is the first comparative study of the gut microbiome associated with different diarrheal pathogens. It presents the first catalogue of different bacterial taxa representing the core and variable microbiome in acute diarrheal patients. The study helped to define a trend in the gut microbiota signature associated with diarrhea and revealed which ARGs are abundantly present and the metagenome-assembled genomes (MAGs) contributing to AMR.