Extracted DNA was sent to the UCONN MARS facility for amplicon sequencing of the 16S rRNA gene targeting the V4 region (primers 515F and 806R) using the Illumina MiSeq paired-end sequencing platform via the v2 2x250 base pair kit (Illumina) following their standard protocol (Parada et al., 2016).