design: The amplification of 16S rRNA gene V4 region was performed using 515F and 806R primers (Caporaso et al., 2012), along with Illumina adapters on the MiSeq platform (Illumina, San Diego, USA). The first PCR reaction followed Caporaso et al. (2011).