Microbiome analysis was performed according to the previous reports.(22 (link),23 (link)) Briefly, the QIAamp Fast DNA Stool Mini Kit (Qiagen, Valencia, CA) was used to extract genomic DNA from the collected fecal pellets. Microbiome analyses by 16S rRNA gene sequencing were conducted by Takara Bio Inc. DNA specimen from feces was amplified using a 16S (V3–V4) metagenomic library construction kit for NGS (Takara Bio Inc., Kusatsu, Japan) with primer pairs 341F (5'-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3') and 806R (5'-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGGACTACHVGGGTWTCTAAT-3') corresponding to the V3–V4 region of the 16S rRNA gene. The amplicons were purified and prepared for the sequencing library by using AMPure XP beads (Beckman Coulter, Brea, CA). The purified amplicon library was sequenced on an Illumina Miseq platform (Illumina Inc., San Diego, CA) at the Biomedical Center at Takara Bio. The processing of sequence data, including chimera check, operational taxonomic unit (OTU) definition, and taxonomy assignment, was performed using QIIME 2,(24 (link)) DATA2,(25 (link)) and VSEARCH.(26 (link))