design: Total DNA was extracted from honey bee gut samples using the DNeasy blood and tissue Kit (Qiagen). The V3V4 region of the 16S rRNA gene was amplified. Libraries were prepared using dual-indexed Illumina adapters and sequenced on the Illumina MiSeq platform with 2  300 bp paired-end reads.