16S rRNA V4 region amplified using 515F-806R primers with Illumina adapters and Golay barcodes. PCR amplification followed by quantification, pooling, and cleanup. Sequencing performed on Illumina MiSeq platform (2x300bp paired-end) using v3 reagent kit. Raw reads processed with QIIME2 v2019.10, quality filtered and denoised using DADA2 to generate ASVs. Taxonomic classification performed using Nave Bayes classifier trained on Greengenes database. Statistical analysis conducted in R using DESeq2 package.