Raw sequences were sorted using assigned barcodes and cleaned up before analysis (barcodes removed and sequences above a quality score, Q ≥ 30 taken forward for analyses). For assembly and annotation of sequences, MetAMOS [28 (link)] pipeline or Partek Flow software (Partek® Flow®, Partek Inc., St. Louis, MO) were used. These softwares provide powerful tools to filter unique hits between human and mouse-specific genes versus microbial signatures. Alpha and Beta diversity calculations were done using embedded programs within the metagenomic pipeline, or using Stata15 (StataCorp LLC, College Station, TX) or EXPLICET software [29 (link)].
Functional profiling was performed using HUMAnN2-0.11.1 [30 (link)] with Uniref50 database to implement KEGG orthologies.