METHODS: Oral swab samples were collected from smokers, smokeless tobacco users, and healthy controls (n = 44). Microbial DNA was extracted and the 16S rRNA gene profiled using the Illumina MiSeq platform. Sequencing reads were processed using DADA2, and taxonomical classification was performed using the phylogenetic placement method. Differentially abundant taxa were identified using DESeq2, while functional metagenomes based on KEGG orthology abundance were inferred using LIMMA.